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Abstract 



Recent developments in the field of high precision calculations in the Standard Model are illustrated, 
with particular emphasis on the evidence for radiative corrections and on the estimate of the theoretical 
error in perturbative calculations. The most important high energy observables and the basic features 
of the renormalization program are discussed in the on-shell and MS frameworks, and a calculation 
of the relation between the electroweak mixing angle sin^^vp(Mz) and the effective coupling sva^ 9^^J^ 
measured at LEP and SLD is considered. I present strong indirect evidence for the contribution 
of bosonic electroweak corrections in the Standard Model, and argue that at the present level of 
experimental accuracy the full one-loop electroweak corrections are needed to describe the data. A 
complete calculation oiO{aas) effects in bosonic vacuum polarization functions is discussed in the most 
general case of arbitrary quark masses and momentum transfer. Compact exact formulas are derived for 
the self-energies of electroweak vector bosons and Higgs scalars, including the case of some extensions of 
the Standard Model; the connection between the calculations based on dimensional regularization and 
on dispersion relations is also investigated. I illustrate some applications to electroweak observables, 
giving an estimate of the error due to the truncation of the perturbative series, which appears to 
be well under control. The first complete calculation of the O^G^m^ M^) contributions to the p 
parameter in the Standard Model is then described for the case of neutrino-lepton scattering, and the 
result used to estimate the theoretical error in the prediction of and sva^ 9^^J^ . In view of the 
experimental precision expected at future colliders, this suggests that a complete two-loop calculation 
of the electroweak corrections is needed. Finally, in the appendices I list useful two-loop integrals, 
asymptotic formulas for the QCD corrections to two-point functions, and Ward identities connecting 
electroweak self-energies. 
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Chapter 1 

Introduction 



High energy experiments have reached in the last few years an impressive accuracy, which 
has called for an analogous precision in the predictions provided by theorists. The mass of 
the neutral gauge boson, for example, is now measured within few parts in 10^, and the 
electroweak mixing angle at the Z mass scale is known at the level of two parts in a thousand, 
with the prospect to reach a relative uncertainty of 5 X 10~^ in the very near future. 

As in the electroweak sector of the Standard Model (SM) the coupling constants are 
small, perturbation theory is generally valid: perturbative calculations beyond the lowest 
order in quantum field theory involve a considerable amount of conceptual and calculational 
complexity, but at present they are the most valuable tool to test the existing theories and to 
explore possible alternatives, in particular in a future of increasing precision. 

The SM has indeed been very successful in describing a large amount of data over a wide 
range of energies, from zero momentum transfer up to the Z mass scale. There are a few 
discrepancies, that need to be carefully considered, but they are at the level of 2-3(t, and 
certainly not compelling. At present there are no glaring contradictions with SM predictions. 
On the other hand, there are strong theoretical arguments indicating that there must be 
physics beyond the SM; the current situation suggests that we need to investigate the very 
fine details of the SM structure before being able to open a window on the physics beyond. 
Precision calculations are therefore of crucial importance: they serve the twofold purpose of 
verifying the SM as a fuU-fiedged Quantum Field Theory, uncovering the evidence for quantum 
corrections, and of allowing the detection of possible deviations from its predictions. 

In this framework, the recent reports of evidence for a top quark by the CDF and D0 
Collaborations at Fermilab [1,2] hit remarkably on the mark: the average of the values for the 
top quark mass given by CDF and D0 (m^ = 180 ± 12 GeV) agrees surprisingly well with a 
recent global fit of experimental data (m^ = 178 ± ll^ig GeV), obtained including all known 
radiative corrections of the SM [3]; these are often sensitive functions of the top mass, and the 
comparison of diff"erent observables allows us to set limits on this unknown parameter. As the 
precision in the measurement of improves, the focus is inevitably going to shift to the other 
standing enigma of the SM: the Higgs boson and the symmetry breaking mechanism, and the 
challenge to experimentalists and theorists would become more difficult, as the dependence of 
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most observables on this sector of the theory is much weaker. 

In this dissertation I will try to illustrate some of the recent developments in the field of 
high precision calculations in the SM of the electroweak and strong interactions. In particular, 
I will address two complementary questions: what is the present experimental evidence for 
quantum corrections beyond the leading fermionic contributions; what is the possible role of 
higher order (beyond one-loop) contributions to electroweak observables. Specifically, I will 
consider the cases of M^i^ and sin^ ^e//' ^^^^ appear to be most promising, among high-energy 
observables, for a significant improvement of the experimental accuracy. 

In Ch. 2 I introduce to two of the renormalization procedures most frequently employed in 
precision physics, and present a detailed calculation of the eff"ective mixing angle measured on 
the Z-resonance by experiments at LEP and SLC. In Ch. 3 I will discuss the phenomenological 
evidence for radiative corrections in the SM. A simple example will illustrate that there is 
very strong indirect evidence from high-energy experiments for the bosonic corrections of the 
theory. Ch. 4 deals with the universal QCD corrections, i.e. the perturbative QCD corrections 
to vector boson self-energies that enter the determination of most precision observables. I will 
present a complete calculation of these eff"ects in bosonic vacuum polarization functions at first 
order in a^, for arbitrary values of the momentum transfer and of the quark masses involved, 
and discuss some of the Ward identities that connect them. In Ch. 5 I will consider potentially 
large two-loop contribution of electroweak origin, focussing on a calculation of the leading and 
next-to-leading contribution to the p parameter in a heavy top expansion. The result is used 
to estimate the theoretical error in the determination of and sin^^^^^. Finally, Ch. 6 
contains my conclusions. 



Chapter 2 



Basic Observables and 
Renormalization 

In this chapter I will briefly introduce the basic features of the renormalization of the Standard 
Model (SM), the most important precision observables, and two renormalization procedures 
most frequently employed in the analysis of high-energy physics, the on-shell and MS schemes^. 
Finally, in Sec. 2.4, 1 will present a detailed calculation of the efl"ective mixing angle measured 
at LEP and SLC in the MS scheme. 

2.1 Input parameters and tests of the theory 
2.1.1 Input parameters 

The electroweak lagrangian is characterized by a number of free parameters: the fermion 
masses, the Higgs boson mass M^, the mixing angles and phase of the Kobayashi-Maskawa 
quark matrix^, and three independent parameters that describe the gauge sector of the SM. 
Two possible choices for the latter are 

g, g', V, or g, M^, M^, (2.1) 

where g and g' are the SU(2)l and U(1)y couplings, v the vacuum expectation value of the 
Higgs fleld, and My^r and Mz the vector boson masses. As the theory is renormalizable, this 
minimal set of parameters can in principle describe the physics at any order in perturbation 
theory. For what concerns the gauge sector of the SM, we therefore need only three input 
quantities, to be chosen among the most precise experimental observables; any additional 
measurement will test the theory. As can be seen from Table 2.1, based on Refs. [1-3,7], at 
present the four most precise high-energy observables are the masses of the Z and W vector 

^ Among the many reviews on the SM, the procedures to renormalize it, and on precision tests, one can refer 
for instance to Refs. [4-6]. 

^In the event that neutrinos are massive, additional parameters have to be considered, masses and mixing 
angles, but neutrino masses are not expected to contribute significantly to precision physics, and I will disregard 
them in the following. 
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observable 


experimental value 


SM fit Ref. [3] 


M^(GeV) 


91.1888(44) 


91.1887 


r'/'(GeV) 


2.4974(38) 


2.4973 


Mw{GeY) 


80.23(18) 


80.32 ± 0.06 ± 0.01 




zU.795(4Uj 


zO.786 


^FB 


0.0170(16) 


0.0153 




U.i4o(^iU j 


U.i4o 


A 


n 1 1 "> 
u. ioo^^ i i j 


n 1 

U. 140 


7?, 
It}) 


u.zzuzi zu 1 


n 91 c;a 

U.ZiOo 


J) 


U. 10001 JO 1 


n 1 79 
u. 1 1 z 


jO,6 
^FB 


0.0967(38) 


0.1002 


^FB 


0.0760(91) 


0.0714 


1 - MllMlivN) 


0.2253(47) 


0.2242 


^^^'e[%{ALR) 


0.2294(10) 


0.2320 ± 0.0003t°:°°°° 


mt(GeV) 


180(12) 


178± llti^ 


M„(GeV) 


^ 65.1 


^ 730 




0.118(6) 


0.125 ± 0.005 ± 0.002 



Table 2.1: High-energy observables compared to their SM expectations. The SM predictions are 
based on a fit including all LEP, SLD, vN , and pp collider data, but not the new rrit measurements; 
the first error and the central value refer to Mjj = 300GeV, and the second error corresponds to 
the variation of the central value when varies between 60GeV and ITeV. 

bosons, the total decay width of the Z boson, and the leptonic asymmetries. In addition 
to these high-energy quantities, and to the masses of the fermions except the top, two other 
electroweak observables are measured at low-energy with great accuracy: the fine structure 
constant a, which is most precisely determined from the anomalous magnetic moment of the 
electron, and the Fermi constant Gfj,, which is extracted from the muon lifetime. In Table 2.2 
I show this set of six observables with their experimental errors. 
The most obvious choice in the present situation is to use 

a (2.2) 

together with the fermion masses and mixing parameters as inputs to derive precise predictions 
for the other observables sin^ ^e//' ^'^^ ™ principle for any electroweak observable. I 

will try to estimate the theoretical errors involved in this process, including the effects of the 
uncertainties of the inputs, which can be however expected to be small. In addition to these 
specifically electroweak input parameters, I will use the value 

a,(M^) = a,(M^)g= 0.118 ± 0.006 (2.3) 

for the strong interaction coupling constant, corresponding to the 1994 world average [8] for 
the MS definition, with five active quark flavors. 
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observable 


experimental value [Ref .] 


uncertainty 




137.0359895(61) [14] 


4.5 X 10"^ 


G, 


1.16639(1) X lO-^GeV"^ [14] 


8.5 X 10-6 




91.1888(44) GeV [3] 


4.8 X 10-5 




0.2316(4) [3] 


2 X 10-^ 




80.23(18)GeV [15] 


2 X 10-^ 


"ptot 


2.4974(38)Gev [3] 


1.5 X 10-^ 



Table 2.2: The six most precise electroweak observables. The last column shows the relative 
uncertainties. The values for and are the 1994 LEP average, sin^ ^^^^ is given as the 
average of LEP and SLC asymmetry values, and Myy is the 1994 collider world average. It does 
not include the value Myy = 80.24(25)GeV from measurements in i/N low-energy scattering. 

A number of natural relations link the couplings and the masses of the gauge bosons 
occurring in the bare SM Lagrangian. Most notably, we see that, although initially defined in 
terms of the gauge couplings, after spontaneous symmetry breaking the electroweak mixing 
angle 9^ relates both bare masses and couplings among themselves: 

tanC = -, eo = ffosinC, = cos0°^M°. (2.4) 
9o 

At the tree level, using \/2Gjti = q"^ l{^^w) — ^^/(^ti"), the above relations can be 
re- written in terms of the input parameters (2.2) 

sin^ 9-„ cos^ 9-„ = —f= , (2.5a) 

Ml = cos^ 0y^Ml. (2.5b) 

Radiative corrections modify the natural relations of the bare Lagrangian, and it is convenient 
to summarize their effect in a compact way [9-11] 

''^ c^l- Ar)= fc^l- Ar) 



V2G^Ml 

Ml= Ml =c?f>Ml (2.6) 

The two different possibilities for the right-hand side refer to the fact that, when quantum 
corrections are incorporated, the relations of Eqs.(2.5) cannot hold true simultaneously, and 
a precise definition of the parameters involved is needed. The connection between the renor- 
malized parameters of the Lagrangian and the experimental inputs is precisely the content of 
a renormalization scheme. It seems reasonable to define the masses as the physical masses 
of the gauge bosons, as they are the quantities that are presently being measured at experi- 
ments. In particular, one considers the complex- valued position of the pole of the propagator. 
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s = ttIq + A[s), where mo is the bare mass and A[s) the vector-boson self-energy. This is 
certainly a gauge invariant quantity, as it is based on the fundamental properties of the S- 
matrix. Writing s = m\ — im2T2, two frequently used definitions are rn^f^y^ = Res = m| 
and rripf^yg = m\ = m\ + T\. It turns out that, in a large class of renormalizable gauges, m\ 
diff"ers from the on-shell definition rn? = ttiq + ReA(m^) by gauge- dependent terms of 0{a^) 
which lie beyond current and expected accuracies (in the complementary class of gauges, a 
gauge- dependent diff"erence arises in O(a^) but it is bounded and very small). Furthermore, 
mi can be identified with the mass parameter measured at LEP [16], and can therefore be 
used in precision calculations. A similar pattern is expected for the mass of the W boson. 

Concerning the mixing angle, diff"erent choices are possible. Unlike the vector boson masses, 
the definition of a renormalized sin^ 9^ is purely conventional, and one can choose either of 
Eq.(2.5) to relate a renormalized mixing angle to observables^, or a diff"erent path, using the 
MS (Modified Minimal Subtraction) scheme [21]. All definitions diff"er by finite calculable 
radiative corrections. 

In Eq.(2.6) s, c, and s, c, are therefore abbreviations for, respectively, the "on-shell" sin^^i^, 
cos 9yfr, and the MS quantities sin^^i^, cos^^i^. The first mixing angle, 9^^, is defined by the 
requirement that the third relation of Eq.(2.4) be valid at any order in perturbation theory [9], 

while the second one, 9yfr, is defined by [20] 

sin^^;.(/x)^^^^. (2.8) 

Here a[fi)^ and q:2(m)ms" &re the electromagnetic and SU(2)l couplings defined at a mass scale 
/X in the MS scheme, a renormalization procedure that I will introduce in the next sections. 
The MS angle 9yfr[fi) is therefore a scale dependent quantity; for Z-physics the natural choice 
is to evaluate it at /x = M^. 

2.1.2 The electromagnetic coupling 

Although the fine structure constant has practically no experimental uncertainty, the use of 
a to describe physics on the Z-resonance is somewhat anomalous: among the three input 
parameters of (2.2), a, defined at = 0, is the only one that involves long-distance dynamics. 
This is witnessed by the appearance of mass singularities in the conventional renormalization of 
the electromagnetic coupling [9]. In fact, physics at the Z-scale is naturally described by short- 
distance parameters defined at the Mz scale, like a[Mz), the conventional electromagnetic 
running coupling evaluated at M| 

a(Ml) = 77^-^ = — (2.9) 

^ l + Ren(^^)(M|)-n(^^)(0) 1-Aa' ^ ' 



^The first solution (or improvements thereof) is used in the so-called scheme [17-19], and the second in 
the on-shell scheme, that I will discuss in detail. 
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where Wy-y'^q"^) is the fermionic photon two-point function^, defined at squared transfer mo- 
mentum g^. Numerically, a(M|) ~ 1/128. Because of the importance of this input parameter, 
the precise knowledge of the denominator of Eq.(2.9) is crucial. Note that any other defini- 
tion of a at the Mz scale would be appropriate; another possibility is the MS definition 
that I will discuss in the next section. Unfortunately, the evaluation of Ren-y-y(M|) — n-y-y(O) 
involves long-distance QCD dynamics in the light quark contribution for which perturba- 
tion theory cannot be used. The standard procedure [9, 12, 13] is to use experimental in- 
put from e+e" hadrons and numerically integrate the spectral function obtained in this 
way (see Sec. 4.5). This method introduces a substantial uncertainty in the determination 
of a(M|) stemming from low-energy data, a fact that has important consequences in the 
analysis of the precision electroweak data. In the following, I will use a very recent state- 
of-the-art determination [13] for the five light quark contribution to the running of a, which 
gives n^7^(0) - Rel6y^^{Ml) = 0.0280 ± 0.0007^ The leptonic contribution to the denomi- 
nator of Eq.(2.9), Aai, keeping all leptonic masses [14] in the one-loop QED expressions, is 
Aai = 0.031419. This implies Aa = 0.0594 ± 0.0007, and a-i(M|) = 128.90 ± 0.10. The un- 
certainty is quite substantial, ~ 7 X 10~^, and aff"ects any determination of the e.m. coupling at 
a high mass scale based on the fine structure constant (i.e. it is scheme independent). In fact, 
we will see that it is a real bottleneck to achieve greater theoretical accuracy. The increas- 
ing experimental precision indeed requires new accurate measurements of the cross-section 
e~^e~ hadrons at low energies. The problem could, in principle, be circumvented by using 
high-energy experiments to determine a(M|), but this entails a severe loss of information and 
is quite premature at present. 

2.1.3 M;^ and sin^ ^^7_f* 

Once the input parameters are specified, as we have done in the preceding sections, predictions 
for SM observables can be obtained and serve as tests of the theory. In the following, I will 
focus on two pseudo-observables^: the mass of the W boson and sin^ ^e//' ^^^^ among 
the most precise quantities measured in high-energy experiments, and their experimental 
precision is expected to improve significantly in the next decade, as can be seen in Table 2.3 

^The convention adopted for the vacuum polarization functions throughout this work is the following: 
Aabiq^ ) is the transverse component of the vacuum polarization tensor for the transition between two vector 
bosons a and b, defined as —i times the standard Feynman amplitude. For the photon, Ajj(q^) = —q^Iijj(q^), 
and in general IiT''{q'^) = —Aab{q'^), in the notation introduced in Ch. 4. 

^Recently, two other re-evaluations of the hadronic contribution to the running of alpha have appeared 
[22,23], and diff'er from Ref. [13] by 1-1. 5<7. However, one of them [23] relies on three-loop perturbative QCD 
up to very low energies, instead of fully exploiting the experimental data, and it is not yet clear whether the 
procedure of Ref. [22] is completely consistent. The discrepancy is potentially significant, but its clarification 
lies beyond the scope of this work. 

^ Mw and sin^ ^'eff' well as a number of other high-energy quantities, are sometimes called pseudo- 
observables, in contrast to real observables such as cross-sections and asymmetries. The former are related to 
the latter by some well-defined set of specific assumptions that constitutes a sort of deconvolution procedure 
and involves a certain amount of theoretical input. The distinction is not unimportant, as this procedure 
introduces a large theoretical error, mainly connected to the initial state QED corrections, that is generally 
taken into account by the experimental groups. 
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observable 


present uncertainty [exp.] 


^ year 2000 


far future 




180[w.a.] 
230[CDF] 


50 ± 20[Tevatron] 
40[LEP200] 


20[NLC] 




4 X 10-4[SLD + LEP] 
10 X 10-4[SLD] 


1 X 10-4[SLD + LEP] 
2.5 X 10-4[SLD] 


6 X 10-5 
[LHC] 


6mt{GeY) 


12[CDF + D0] 


4[Tevatron] 


2[LHC] 



Table 2.3: Precision goals at future colliders. NLC stands for the Next Linear Collider, a 500GeV 
e+e" collider under study. "Tevatron" implies extensive upgrades of the collider and detectors, 
and the start of Main Injector. 

(see Ref. [24]; the estimate of the error on sin^ ^^^^ at LHC is taken from Ref. [25]). They are 
less affected by QCD corrections than r^°'; the QCD corrections to and sin^ ^^^^ will be 
studied in Ch. 4. 

The effective weak interaction angle employed by the LEP groups is defined by 

l-4sin2^^7; = Re (2.10) 

9 A 

where gy and are the effective vector and axial couplings in the —>■ II amplitude at res- 
onance (with the exclusion of photon exchange effects), and t denotes a charged lepton [3,26]. 
sin^ ^^^^j is extracted from the on-resonance left-right and forward-backward asymmetries. 
In the determination of the effective angle, it is common practice to use the h and c quark 
asymmetry data together with the leptonic ones, as this procedure is practically model inde- 
pendent. In order to establish the connection with the MS parameter sin^ dw{Mz), or with s, 
the on-shell sine, we note that this amplitude is proportional to [27] 

< u\j^\o > = -ua^ 

where the electroweak form factor k{q^) and can be replaced by ki^q^) and in the on-shell 
scheme, and vi and ui are the lepton spinors^. The actual radiative corrections that enter into 
the determination of these two pseudo-observables therefore depend on the renormalization 
scheme. In the following I will present the case of the on-shell and MS schemes. 

2.2 On-shell scheme 

In the on-shell scheme [9,28], in analogy with the traditional renormalization of QED, the 
fundamental parameters are chosen to be a and the physical masses of all the particles of the 

'^Alternatively, one could define eff'ective couplings directly through the so-called hare asymmetries A'^ , 
derived from the experimental asymmetries after subtraction of soft photon and photon exchange effects. For 
example, ^^'^ = 2x/[l + e^), with x = Qv/Qa- The difference comes from the imaginary part of fc(g^), but we 
will see in Sec. 2.4 that the it is numerically very small. 



75 



(2.11) 
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SM and, in particular for the gauge sector, and Mw From the experimental values of 
the vector boson masses, we find (^5^/5^)exp ~ — 2(c^/s^)(^Mpi^/Mpi^)exp ~ 1-6%, a very poor 
precision for a fundamental coupling. As the experimental precision on Myy is not very good 
at present, and even after LEP200 will not be better than the one on Gfj,, one has to resort to 
the relation, Eq.(2.6), 

Ml\ Ml) V2G^M|[1- Ar(M;^,m4,M«)]' ^" ^ 

where /S.ri^My^r, rrit, M^) is a finite radiative correction that includes the term Aa defined above 
and the non-QED electroweak corrections to the muon decay [9]*. The inclusion of (1 — Ar) 
in the denominator of Eq.(2.12) eff"ectively resums all leading logarithms of 0(a" log" (in 
this way replacing a with a(M|)), and to good approximation the terms of 0(a^ log [29]^. 
Eq.(2.12) can be solved iteratively for and yields a prediction for or sin^ 9^ for any 
given value of rrit and M^- The values for and derived in this way are then used in 
the calculation of any other observable in the on-shell scheme. It is interesting to show the 
asymptotic behavior of the function Ar[mt, M^), as the comparison between the Myy^mt, M^) 
obtained from Eq.(2.12) and allows us to set limits on the possible range of and M^- 

The leading asymptotic behaviors for large rrit and are 

3a m? , rrit 

Ar ~ - + ... f — ^ > 1) 

11a , f Mi\ Mh 

Note that the dependence is quadratic in rrit and logarithmic in Mjj, and that the eff"ects have 
opposite signs. This partially accounts for the fact that, in global analyses of the electroweak 
data, small Mjj values favor relatively small values of rrit] it also explains why it is much 
more difficult to set bounds on Mjj than on rrit. More generally, the whole set of precision 
observables can be calculated as a function of and Mjj, and the result fitted to the existing 
data, yielding the bounds for rrit and Mjj shown in Table 2.1. 

In the on-shell framework it is natural to define the renormalized mixing angle according 
to Eq.(2.7). That definition is manifestly based on short-distance parameters, and it therefore 
absorbs large photon vacuum polarization eff"ects that would be present adopting, for instance, 
the first of the relations in Eq.(2.5). The use of the on-shell mixing angle in Born amplitudes 
will not induce large loop eff"ects of that kind. However, Eq.(2.7) does not absorb m sm u-^ 
large isospin violation eff"ects that occur in the on-resonance amplitudes, and are due to the 
mass splitting in the top-bottom isodoublet. These manifest themselves in the counterterm of 



'As Mz, Mw , and a are physical observables, it follows that the same is true of Ar. 
^In the on-shell scheme the incorporation of leading higher order reducible contributions of 0(G^m() is 
subtle. It can be implemented following the procedure outlined in Ref. [30]. 
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sin^ 9v 



6M^ _ SMl 



Ss' = -c^Y, Y = f^, (2.14) 



where SM^ = ReAww{M^) and SM^ = ReAzz(M|) + {■J,Z) mixing terms, are the mass 
counterterms of the vector bosons. At one-loop level we find Y = Xf + . . ., where Xf = 3 ^^'^ 
is the leading top correction to the p parameter [71], and the ellipses stand for the rest of the 
one-loop contribution. This problem is frequently circumvented by introducing an eff"ective 
parameter = [1 + [c"^ / s'^)xt]s'^ or, equivalently, = c^(l — a;^). 

The determination of Myy through Eq.(2.12) is aff"ected by the parametric uncertainty on 
the inputs Gfj,, a, and M^, and by the theoretical uncertainty on Ar, due to the truncation 
of the perturbative series. For instance, a change SAr induces a shift in the determination of 
Myy through the following relation 

= f w/"^: y (2.15) 

2 [c' — s'^ — zc^Xt) (1 — Aa) 

where some higher order eff"ects have been retained [36]. Using the result of Ref. [13], we 
see that the uncertainty on the hadronic contribution to the running of a implies ^M^jf <^''- ^ 
13MeV. This is the ultimate limit on the precision of a theoretical determination of My^r in 
the present situation. In the next chapters, however, we will see that the theoretical error 
due to electroweak higher order eff"ects in this determination can be of the same size. We 
can compare this uncertainty with the one coming from the experimental error for the input 
parameters. For instance, in the case of M^, 8Mz — 4.4MeV. From Eq.(2.12) it is easy to see 
that this induces on My^r an uncertainty 8My^r ~ 5.4MeV, smaller than the one induced by 



2.3 MS scheme 

The on-shell scheme has the privilege of simplicity, but as we have seen above, it induces 
potentially large higher order corrections in the on-resonance amplitudes. The on-shell cou- 
plings do not seem therefore adequate expansion parameters unless one employs the eff"ective 
couplings a(M|) and . A very good alternative is provided by the MS scheme, which is the 
prevalent framework for QCD studies. 

The on-shell definition of the mixing angle is based on physical quantities; the renormalized 
parameter is set to be equal to an observable quantity, the ratio of the physical masses of the 
vector bosons. On the other hand, it is possible to define a renormalized parameter by choosing 
the counterterm associated to it. The MS definition of a renormalized parameter [21] is based 
on a theoretical prescription for the associated counterterm and requires that the calculation 
be performed in one particular regularization framework, i.e. dimensional regularization. The 
idea is to subtract only the divergent part (poles in n — 4) and the associated constants 7 and 
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ln4x from the bare parameter^", as these are fixed exclusively by the bare Lagrangian, unlike 
the finite parts of the counterterm, that depend on the renormalization condition. The result 
is a scale dependent quantity, and a first example is a[fi), the MS e.m. coupling constant, 
defined by 

where A^(/x) = ll^'y(l{0)^+ ^(7/2 ln(Mp^//x) - 1/6), with the subscript MS to indicate that 
l/(n — 4) poles have been subtracted with their related constants. In addition to the fermionic 
self-energies, A-y contains the contribution of W bosons. The hadronic part is evaluated 
adding and subtracting II^^'' (M|), which can be calculated perturbatively, and using the 
value n^^^O) - ReIlir^{Ml) = 0.0280 ± 0.0007 [13]. For Mz and as given in Table 
2.2, rrit = 175GeV, and including two-loop QCD eff"ects following the results of Ch. 4, we 
obtain A^{Mz) = 0.0653 ± 0.0007 and d-'^{Mz) = 128.09 ± 0.10. However, rut > Mz, and 
according to the common practice in the MS scheme, we may subtract the top eff"ects from 
Eq.(2.16), absorbing them in the definition of a, as it would be the case for another heavy 
unknown charged particle. Decoupling the top quark from the definition of d, we therefore 
have A^{Mz) = 0.0666 ± 0.0007, and ^-^(M^) = 127.91 ± 0.10, which diff"ers from a(M2), the 
conventional running coupling constant of Sec. 2.1.2 by 0.73%, and has a similar uncertainty 
of ^ 7 X 10"^ 

An analogous definition can be employed for the other couplings and masses. The bare 
relations Eq.(2.4) hold for all MS couplings and masses at a given scale /x. In particular, we 
can relate the bare sine to the MS definition of Eq.(2.8) through (the ellipses indicate higher 
order contributions) 



sin^ 9w{iJ')ms = sin^ Ow 



oo /n , 19 



4x I 3 6 sin^ 0° / e 



(2.17) 



Here I have subtracted only the poles and the related constants (not indicated), but one 
can think of implementing the decoupling of heavy particles as done above for a. How- 
ever, in this case a subtlety arises, as Ci2{n)^ can be defined by charged or neutral current 
couplings, and the top aff"ects them diff"erently. According to the Marciano-Rosner [M-R] con- 
vention [38], adopted also in Ref. [36], the logarithmic term is subtracted in the evaluation of 
Re a''^°^\mI)/MI, which enters the form factor ki{Ml) of Eq.(2.11) as 



cRe 47^(M|) 



a 



Ml 6x52 



1 - -5 



as\ 15 

— log6 - 

X / 4 X 



D{t)- (2.18) 



Here = ^?/-^| and D[l/^t) represents small terms that decouple in the limit oo- 
We see that at two-loop there is also an m^-independent term that should be subtracted as 
well. This additional finite counterterm keeps the 7 — Z mixing independent of rrit in neutral 



^"This is consistently done at any order in perturbation theory by rescaling the 't Hooft mass scale n 
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current amplitudes and allows sin^ ^w(m)ms" to be continuous at /x = rrit, up to small 0{aas) 
terms. The aim of the prescription is to make the value of sin^ dw{Mz), as extracted from the 
on-resonance asymmetries, very insensitive to heavy particles of mass M > M^. As can be 
seen from Eq.(2.18), the numerical difference between sin^ 6w(Mz) evaluated with or without 
the M-R decoupling convention is very small, but not negligible, ~ 0.0002 in the present rrit 
range. 

The use of MS couplings would suggest the use of MS masses Mz{fi) and Myy^fi) for the 
vector bosons^^. However, as we have seen, the experiments provide the values of the physical 
or pole masses of the vector bosons, and the relation between M and M(/x) involves large 
radiative corrections of 0{Gfj,m^\nmt/ fi). At high orders in perturbation theory, it can be 
seen that a whole class of potentially large 0(G^mj"~^M^) contributions, corresponding to 
fermionic loop insertions on vector boson propagators, can be absorbed in the definition of 
the on-shell masses. The on-shell mass counterterms cancel completely the leading quadratic 
dependence of the heavy fermion loop. Thus, the use of MS couplings together with on-shell 
masses in the propagators seems preferable in order to keep higher order radiative corrections 
small [11]. An explicit example of this kind will be presented in Ch. 5, in the case of the 
two-loop 0(G^mjM|) contribution to the p parameter, which turns out to be smaller when 
employing this prescription. 

If we choose this route, the auxiliary parameter sin"^ 9 w{Mz) is derived from the inputs 
setting the mass scale equal to Mz, and using 

s^c^ = , (2.19a) 

V2G^Ml[l-Ar{M^,muM^,s)y 

and Myy is determined through Eq.(2.6) 

Ml = p{My„mt,MH,s) Ml (2.19b) 

Eqs.(2.19) can be solved simultaneously by iteration. The radiative correction 
p{mt, Mjj, Myy, s), which now governs the Myy-Mz interdependence, incorporates the custodial 
symmetry violation eff"ect, quadratic in m^. p can be readily obtained from the diff"erence 
between the two definition of sin^ ^Pi^. Using Sq = s^ — Ss^ = — ^i^, we have = (1 — i^Js")? 
where MS indicates that ultraviolet poles and related constants have been subtracted from 
Y in Eq.(2.14), and p, has been set equal to Mz- It then follows that y5 = (1 — Y^)~^. The 
function Ar is obtained calculating the electroweak corrections to the muon decay in the 
MS framework [10,11]. The asymptotic behavior of these two functions is similar to that of 
Ar, but the top dependence of sin^ 6w(Mz) is milder, as part of it has been absorbed in the 
definition [10] 

Ar ~ -2 S + {mt>Mz). (2.20) 

Explicit formulae including some higher order eff"ect can be found in Refs. [11,36]. The use 
of a instead of a in the numerator of Eq.(2.19a) implies the appearance in Ar of the mass 



^This approach is illustrated in the first article of Ref. [6]. 
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singularities contained in A-y, and the corresponding hadronic uncertainty. The latter affects 
the determination of P according to 

which leads to an uncertainty ~ 2.5 X 10~^ due to the hadronic contribution to the running 
of a; this is one order of magnitude larger than the one induced by the uncertainty in M^. 
The determination of M^i^ through Eq.(2.19b) differs by the one through Eq.(2.12) by higher 
order subleading contributions, when the leading rrit terms are consistently implemented. 
A numerical analysis [36] of this scheme dependence has shown a maximum discrepancy of 
^ 4MeV. 

An important part in the choice of the renormalization framework is played by the sim- 
plicity of the formulation. For instance, the use of sin^ ^^^^ as definition for a renormalized 
coupling would complicate beyond reason the analytic computation of electroweak radiative 
corrections, as this parameter is based on a very specific physical amplitude. If this may ap- 
pear as an aesthetic point of view in one-loop calculations, it becomes a vital matter in higher 
order calculations. The MS framework does provide a clear and simple environment for such 
involved computations^^. 

Finally, it is worth mentioning that a precise determination of the MS parameters at the Z 
scale is a very convenient basis for applications to Grand Unified Theories (GUT), and for new 
physics in general. In particular, using the renormalization group, it is possible to relate the 
values of these parameters at Mz to the ones at a very high energy scale /x, and in this way to 
test different scenarios of Grand Unification. For example, the minimal non-supersymmetric 
SU(5) GUT [31] predicts = 3/8 at the unification scale, which extrapolates to P = 0.2100 ± 
0.0025 ± 0.0007 at Mz [32] using a and as inputs. Together with proton decay data, this 
prediction rules out the simplest GUTs. In contrast, in the minimal supersymmetric SM 
(MSSM) Grand Unification yields sin^ Ow{Mz) = 0.2334 ± 0.0025 ± 0.0025, which is clearly 
compatible with present data. In practice, as the precision on sin^ 6w(Mz) is now much 
better than the one on as{Mz), it is standard procedure to use the MSSM determination 
of sin 9yfr(^Mz) together with ^(Af^)? and to predict ^^(Af^)? assuming Grand Unification. 
Present data are compatible with this hypothesis [32]. This is just one illustration of the 
importance of precision calculation for uncovering new physics. 



2.4 The leptonic effective sine in MS 

As we have seen above, the MS definition of sin^ 9^ seems particularly suitable for the de- 
scription of physics on the Z-resonance. It has been known for some time, for example, that 
the use of MS parameters in the Born amplitudes provides a very good approximation of the 
neutral current complete amplitude^^, and in particular sin^ dw{Mz) is very close numerically 

A comprehensive comparison of different renormalization schemes for the electroweak sector of the SM can 
be found in the second article of Ref. [6]. 

^^An exception is the Zbb vertex [35], where large non-universal top terms are also present. 
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to the parameter sin^ 6^J'^ extracted from the asymmetries. However, the reason and extent 
of this coincidence and the precise numerical relation between the two has not been spelled 
out in the literature. In view of the very accurate experimental determination of sin^ ^e//' 
it is important to know the precise relation between the two mixing angle. In this section, I 
present a detailed calculation [34] of the form factor fc(M|) which controls this relation. 
Up to terms of order 0(d) in the MS framework we have [27] 

" 77MSV" / 

where A^zi^"^) is the 7 — Z mixing self-energy, the subscript MS means that the MS renormal- 
ization has been carried out (i.e. the pole terms have been subtracted and the 't Hooft scale has 
been set equal to Mz), the superscript / stands for fermionic part, a{Mz) = [127.91 ± 0.10]~^, 
= M^/MI and Vi{q'^) is a finite vertex correction. Explicitly, 

where f{x) and g[x) are defined in Eqs. (6d) and (6e) of Ref. [27]. I have included the photon 

77MS * 



self-energy A^^^-^il"^) in the second term of Eq.(2.22) because, as it will be explained later, it 



gives rise to relatively large O(a^) terms. 

It is clear from Eq.(2.11) that the ratio of the vector and axial vector couplings at resonance 
is given by 1 — 4 ^^(Ml) P. We can now discuss the various contributions to k^^Mz). 

To 0(d) the fermionic contribution to the self-energy of Eq. (2.22) can be written in the 
form 



^ Qi ^75 ' (2-24) 



s Ml X52 ^ V 4 Ml 



where Qi, Ci, and rrii are the charge, third component of weak isospin (with eigenvalues ±1), 
and mass of the i-th fermion, the summation includes the color degree of freedom, 11^ is 
the vacuum polarization function involving vector currents, according to the normalization 
of Ch. 4, and henceforth the MS renormalization is not indicated explicitly. For the leptons 
we can safely neglect the masses and using the result in App. B find that the contribution 
to Eq.(2.24) is (d/xs2)(5/12 + ix/4)(l - 45^) = (3.3 + 6.2i)10"^. In this calculation and 
henceforth I employ = 0.2317, which corresponds to the central values rrit = 178 GeV and 
ruH = 300 GeV in the global fit of Ref. [3], and d = (127.9)"^ 

For the first five quark fiavors I set again = (it can be verified that this is an excellent 
approximation) and, including 0{das) corrections available from the expressions in App. B, 
obtain a contribution (d/xs^) (7/12 - 11 P/9) [5/3 + (a^/x) (55/12 - 4C(3) + iir)] = 
(5.35 + il0.51) X 10"^, where I have used = a^iMz) = 0.118 and C(3) = 1.20206.... 

The top quark contribution to Eq.(2.24) can be derived from Eq.(2.18), and the complete 
expression can be again gleaned from App. B. For the current range 150GeV^m4^210 GeV [3], 
D[l/^t) varies from 6 X 10~^ to 3 X 10~^ and is of the same order of magnitude as neglected 
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two-loop contributions ~ (d/xs^)^ ~ 10 ^ to Eq.(2.22). We reach the conclusion that when 
the M-R prescription is applied, the top quark contribution to Eq.(2.24) is very small. 

The other contributions to ^^(Ml) in Eq.(2.22) can be readily obtained from the literature. 
This form factor is gauge invariant, but several individual components are not. I evaluate them 
in the 't Hooft-Feynman gauge, using Myy = 80.23 GeV [3]: (i) the bosonic contributions 
-(c/i) - ^^^(O)]/^! ^an be extracted from Ref. [39] and amount to -5.94 X 10"^; 

(m) -(d/4xs2) Vi{Ml) can be obtained from Eq.(2.23) above and Eqs.(6d,e) of Ref. [27], and 
gives +(3.33 + 2. 78i) 10"^; (m) (d/xs^)^^ log = -2.11x10"^; {iv) although two-loop effects 
have not been fully calculated, I include the O(d^) contribution arising from the product of 
ImA^z{M^) and ImA^^(M|) in the second term of Eq.(2.22). The ratio / {M^ - Ai/-y\M^)) 
is given by 0.992 — 0.0175i for the current range of rrit, and the interference between the two 
imaginary parts gives a contribution to Eq.(2.24) of ~ 1.9 X 10~^. It is quite sizable, relative 
to typical O(d^) contributions, because these imaginary parts involve several additive terms. 
On the other hand, the large logarithmic O(a^) corrections associated with the running of a 
are already taken into account, in the MS scheme, by employing d in the evaluation of A^z- 

Combining all the above results we have ^^(Ml) = 1 + (0.33 + 0.62i + 5.32 + 105.1i-5.94) X 
10"^(0.9921 - 0.0178i) + £'(1/6) + (3.33 + 2.78i - 2.11) x 10"^, which to good approximation 
becomes 

h{Ml) = 1.0012 + i 0.0139, (2.25) 

where the real part slightly decreases for a heavy top. It is clear, on the basis of Eq.(2.25), 
that at the one-loop level the ratio of effective vector and axial vector couplings in the Z ^ ££ 
amplitude is a complex number. This is also expected from general principles. On the other 
hand, the LEP groups interpret both sides of Eq.(2.10) as real quantities. This can be justified 
on the grounds that the imaginary component of k^^M^) gives a negligible contribution to the 
leptonic T-even bare asymmetries and partial widths For instance, the bare left-right 
asymmetry is given by A^£^ = 2 Ke(gv / gA)/[^ + Iffy/ffjlP] and one readily verifies that the 
inclusion of Im ki(Ml) decreases its value by only —0.02%, and S sin^ ^^^^ ~ —2%S A^£^/ A^£^ 
for the current value of sin^ ^e//' Similarly, is modified by ~ —0.03%. Therefore we 

identify 

sin^ 0[J} = Re ki{Ml). (2.26) 

Using Eq.(2.25) we have 

sin^ 0[J^^ - sin^ 0w{Mz) = 2.8 X 10"^ ^ 3 X 10-^ (2.27) 

which tends to approximate to 2.7 X 10~^ for a heavy top (m4^185GeV).The following ob- 
servations are appropriate at this stage: (a) because the Higgs boson does not contribute at 
the one-loop level to Eq.(2.22), the results of Eqs.(2.25, 2.27) are independent of ttih; (6) it 

^^The imaginary part of ki(m%), however, gives important contributions to T-odd leptonic asymmetries. See 
Ref. [40]. 
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is clear that the closeness of Re k^^M^) to unity and, correspondingly, the small difference 
in Eq.(2.27) are due to the cancellation of significantly larger terms. For instance, the light 
quark and bosonic contributions to the f-Z mixing self-energy are of the roughly expected 
order of magnitude ~ a/(2xs^) ~ 5.4 X 10~^, but they largely cancel against each other. On 
the other hand, the 0{a) contributions to the Im ki{Ml) are ~ 1%, an order of magnitude 
larger; (c) as the relation between sin^ dw{Mz) and sin^ dw = 1 — M^/M^ is well-known (see 
Sec. 2.3), Eq.(2.27) determines the connection between the three parameters. 

If the Marciano-Rosner decoupling convention is not applied, so that in the MS renormal- 
ization one only subtracts poles and sets the 't Hooft scale equal to M^, there is a further 
contribution to Re k^^M^) arising from the first term in Eq.(2.18). Using as{mt) ~ 0.11, this 
amounts to -6.1 X 10"^, -8.7 X 10"^, -1.09 X 10"^, for rut = 150, 180 and 210 GeV, re- 
spectively. Correspondingly, Re k^^M^) becomes 1.0006, 1.0003, 1.0001, even closer to unity. 
As a consequence, although the diff"erence between sin^ ^^^^ and sin^ 9w{Mz) depends more 
on rrit when the M-R prescription is not applied, it is actually smaller for the current range 
150^mt^210GeV. In fact, we find that it is 1 X 10"^ for 150 < rut < 191 GeV, and there is 
no diff"erence in the fourth decimal for 192 < rrit < 210 GeV. 

One can obtain a rough consistency check of the order of magnitude of Eq.(2.27) by 
comparing the fits of Ref. [3] (Table 2.1) with the calculations of Ref. [36]. Using the LEP, 
SLD, collider and v data, Ref. [3] finds rut = 178±11+J9 and sin^ e[j*f = 0.2320±0.0003+° °°°2- 
Their central values assume ttih = 300 GeV, the first error represents experimental and 
theoretical uncertainties, while the second refiects changes corresponding to the assumptions 
rriH = 60 GeV and Mh = 1 TeV. According to Eq.(2.27), the corresponding central value for 
P should be 0.2317. On the other hand, from Ref. [36] one finds P = 0.2317 for rut = 178 GeV 
and rriH = 300 GeV. Thus, the comparison of the conclusions of Ref. [3] with the calculations 
of Ref. [36] is roughly consistent with Eq.(2.27). Of course, such consistency checks are not a 
substitute for precise, ab initio calculations, like the one leading to Eq.(2.27). 

Two additional comments can be useful at this stage: (a) Consistently with Eq.(2.10), 
sin^ ^^^^ can be defined in terms of a bare forward-backward asymmetry A^^, which is ob- 
tained from the physical asymmetry Ap^ after extracting the eff"ect of photon-mediated con- 
tributions and other radiative correction eff"ects [26]. Therefore, we should not attempt to find 
the numerical relation between sin^ ^^^^ and sin^ 6w{Mz) by comparing detailed MS calcula- 
tions of the physical asymmetry Ap^, as those in Ref. [27], with theoretical expressions for 
A^pg expressed in terms of sin^ ^e//' "^^^ point is that Ap^ contains electroweak eff"ects not 
contained in A^^. (6) Global analyses often cite the value of sin^ ^^^^ as extracted only from 
the on-resonance asymmetries, while they give the prediction for rrit derived from the complete 
data base. Current asymmetry results lead to determinations of sin^ ^^^^ that are somewhat 
smaller than the sin^ dw{Mz) numbers corresponding to the central rrit. This, however, is not 
a contradiction with Eq.(2.27), because the on-resonance asymmetries represent only a part of 
the experimental information. This is quite visible in the detailed report of Ref. [3], in which 
one finds sin^ ^^^^ = 0.2317 ± 0.0004 from the on-resonance asymmetries and, as mentioned 
before, larger values from the global fits. 
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In summary, I have attempted to clarify the connection between sin^ 9^^^ and sin^ dw{Mz) 
and obtained the value of their difference by means of a detailed calculation, both with and 
without the M-R decoupling convention [34]. I recall that sin^ ^^^^ is defined in terms of a 
physical amplitude, while this is not the case for sin"^ 9 w{Mz), the running MS parameter 
evaluated at the M^-scale. Thus, it is clear that the two parameters are conceptually very 
diff"erent. On the other hand, we have found that their numerical values are very close (cf. 
Eq.(2.27)). As illustrated in the discussion after Eq.(2.27), this fact is probably due to a 
fortuitous cancellation of radiative corrections. 

2.5 Higher order corrections 

In the preceding sections I have outlined the main features of the renormalization program for 
electroweak high-energy physics. Although the bulk of the radiative corrections is doubtless 
represented by the one-loop contributions, it is necessary to consider all possible sources of 
large higher order contributions. In some case, as for the electromagnetic coupling, renormal- 
ization group techniques help to resum large contributions at any order. Whenever they are 
possible, however, explicit multi-loop calculations are the best way to investigate the conver- 
gence of the electroweak perturbative series. In fact, what we really need is an estimate of 
the theoretical error due to the truncation of the perturbative series. There exist heuristic 
methods for estimating the higher order contribution in QCD, but in general no such method 
applies to the more complicated situation of electroweak interactions. In the last few years, 
multi-loop calculations in electroweak physics have become increasingly important, as the 
lower bound on the top mass increased, and the experimental accuracy improved. One of 
the eff"ects of the activity in this field has been the realization that the convergence of the 
electroweak series is not as smooth as it was believed: large higher order eff"ects do exist, and 
must be taken into account. In the following I will restrict the analysis to the two observables 
considered in this chapter, and sin^ ^e//' distinguish two kinds of possibly large 

higher order eff"ects that can aff"ect the predictions for these observables. 

• Large mass effects and rrit are heavy, of the order of the spontaneous symmetry 
breaking scale. Heavy masses do not decouple in the SM and therefore we can expect 
loop contributions that are not bounded as rrit and go to infinity. In particular, the 
so-called decoupling theorem of Appelquist and Carazzone [43] does not apply in this 
context. If we exclude the case of amplitudes involving external Higgs bosons, we can see 
that at the one-loop level the radiative corrections scale as the square of the top mass, 
and as the logarithm of M^, as in Eq.(2.13)^^. At two-loop level, the leading heavy 
mass eff"ects are of 0{G'^mf) and 0{aGfj,M^). Since complete two-loop calculation do 
not exist at present, a good starting point would be to have all power-like heavy mass 

^^The logarithmic dependence on Mh of the one-loop radiative corrections in the gauge sector of the SM has 
been first shown in Ref. [41]. In the limit of large Mh the SM can be viewed as a gauged non-linear <7-model 
and the mass of the Higgs boson can be interpreted as an ultraviolet regulator. It is then easy to see that at 
one-loop the non- linear <7-model contains only logarithmic divergences [42]. 
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corrections under control. This kind of effects will be considered in Ch. 5. 

• QCD effects Because of the large coupling constant, perturbative QCD corrections 
can be sizable, even at a high energy scale. In particular, the determination of Myy 
and sin^ ^^^^ proceeds through physical amplitudes (muon decay and Zll vertex) that 
involve only oblique QCD corrections, i.e. QCD corrections to quark loop insertions on 
vector boson propagators. We have seen that large quadratic top effects can arise as a 
consequence of the SU(2)i^ symmetry breaking in the (t, b) isodoublet; we can therefore 
expect relevant QCD corrections of 0{Gfj,asm^). This kind of effects will be studied in 
Ch.4. 

As we have seen, in the MS framework, the determination of sin^ ^^^^ proceeds through 
Eq.(2.27) and Eq.(2.19a). The relation between sin^ O^J^ and sin^ 0w{Mz), however, has 
already been evaluated keeping leading QCD effects in the hadronic self-energies, and, if 
the Marciano-Rosner prescription is employed, should be free of large power-like heavy mass 
effects of 0(G^mjM|). Therefore, the precise determination of M^i^ and sin^ dw{Mz) through 
Eq.(2.19) becomes the most important task in precision physics. Some aspects of this effort 
will be illustrated in the following. 

But first, in the next chapter, I will present a simple method to attain indirect evidence 
that a full one-loop calculation is needed to describe current data. 



Chapter 3 

Evidence for radiative corrections 



We know that a renormalizable quantum field theory includes quantum corrections of 0{h), 
and we know that the SM is very successful in describing current experimental data. One 
might then wonder what do we know phenomenologically about the quantum field structure 
of the SM. In this chapter I will discuss the current phenomenological evidence for radiative 
corrections. A simple example [33] will illustrate that there is very strong indirect evidence 
from high-energy experiments for the leading contributions due to fermionic loops, as well as 
for the numerically subleading contributions that involve virtual boson exchange. 

3.1 Evidence from low-energy 

The fundamental importance of radiative corrections in the electroweak theory has been known 
for a long time, and has been established first by low-energy experiments [47]. For example, 
in the absence of radiative corrections, a large violation of the unitarity of the Cabibbo- 
Kobayashi-Maskawa matrix (which is mainly determined from low-energy experiments) would 
appear. When all the one-loop electroweak corrections are included and the nuclear overlap- 
ping is appropriately taken into account, an analysis [44,45] of the ^^O super-allowed Fermi 
transitions leads to Vud = 0.9745 ± 0.0005 ± 0.0004. and Vub are much smaller, and are 
extracted from AS" = 1 and B decays, yielding [45] 

\Vud\^ + \Vus\^ + \Vub\^ = 0.9983 ± 0.0015, (3.1) 

very close to unity. The radiative corrections aff"ecting the derivation of Eq.(3.1) are quite 
sizable, about 4.1%, and involve vector boson exchanges. As the momentum scale is small, no 
large electromagnetic corrections due to the running of a are present, and the sensitivity to 
rrit is cancelled in the ratio between the amplitudes for /3-decay and muon decay, so that we 
can ascribe this large contribution to bosonic eff"ects. When we analyze the data employing 
only tree level expressions and the Fermi-Coulomb function, we obtain 1.0386 ± 0.0013 [47], 
which is about 30a away from unity. In fact, the appearance of a large radiative correction in 
the determination of Vud can be traced back to the observation that the photonic corrections 
to the /3-decay in the local V-A electroweak theory are not convergent, unlike the ones to the 
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muon decay^. They involve a logarithmic divergence ~ \n—, where nip is the mass of the 
nucleon. In the one-loop SM calculation [44], the cutoff is simply replaced by M^, yielding a 
sizable contribution. 

On the other hand, one could argue that the effect of SM loop corrections that go beyond 
the QED corrections can be mimicked by simply setting the cutoff equal to M^. However, 
in order to match the experimental precision, an accurate determination of the cutoff A in 
the Fermi effective field theory would be necessary. Although the order of magnitude of 

— 1/2 

the cutoff can be easily guessed (the symmetry breaking scale v ~ Gfj, ), without a precise 
determination the local V-A theory result for Eq.(3.1) can still be very different from unity [50]. 
Such a precise determination can only be provided in the framework of the full renormalizable 
SM. In conclusion, it appears that low-energy experiments provide a very strong evidence for 
electroweak radiative corrections of bosonic origin. 

3.2 Evidence from high-energy 

While the low-energy corrections in the problem of universality involve virtual fermions and 
W^, 7, and Z bosons, at high-energy the dominant corrections come from fermionic loops 
(an exception is the Z ^ bb vertex, where large corrections are induced by loops involving 
virtual bosons). As we have seen in the previous chapter, one of the main sources of large 
corrections is the running of a, with the associated large logarithms. Evidence for radiative 
corrections beyond the running of a has been analyzed in Ref. [47]. Direct evidence of such 
contributions comes from the comparison of the experimental data for specific observables with 
the predictions of a Born approximation in which a has been replaced by a[Mz) or d. This 
kind of direct evidence is not particularly compelling, and presently reaches ~ 1.5 — 2a [48]. 
However, it has been argued [47] that very strong (at the level of 5-6a) indirect or inferred 
evidence can be uncovered when the comparison is made between the Born approximation 
and the results of a global fit in the full SM to the whole body of electroweak data. The global 
fit provides more precise values than experiments for single observables, and incorporates all 
the interlocking relations that highly constrain the SM. 

It is also interesting to note that no Born approximation can accurately describe all the 
available experimental electroweak data. Indeed, we have seen in the previous chapter that two 
very different values for sin^ 6^ can be extracted from experiment: sin^ ^^^^ and 1 — M^/M|, 
respectively 0.2317 ± 0.0004 and 0.2259 ± 0.0035; they differ by 1.7a. No Born approximation 
for the electroweak mixing angle can accommodate both values. Again, the discrepancy is 
much stronger if we compare the values for these observables in the SM fit of Table 2.1: we 
find sin^ e[j*f = 0.2320 ± 0.0004 and I -Ml I Ml = 0.2242 ± 0.0012, which differ by 5.6a. Note 
also that the fit of Table 2.1 does not include the latest CDF and D0 values for m^. 

beautiful and simple explanation of this difference, based on Fierz transformations, can be found in 
Ref. [46]. 
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Figure 3.1: Determination of sin"^ 9 w{Mz) from asymmetries (horizontal lines), and from Mz 
(slanted curves), in the full SM with = 300GeV. The la errors are indicated. 

3.2.1 Bosonic radiative correction 

Besides the running of a, another source of large fermionic contributions is the sharp breaking 
of isospin symmetry in the t — b quark doublet, which introduces terms that are quadratic in 
the top quark mass, and allows the determination of rrit through radiative corrections. For 
instance, the radiative correction Ar, introduced in the previous chapter (Eq.(2.12)), can be 
decomposed in 

Ar = Aa - —xt + Ar^em (3.2) 

where Aa and Xt = 3 ^ are the large fermionic contributions coming from the running 
of the e.m. coupling constant and from the t — b doublet, while Ar^em collects all other 
contributions, mostly bosonic in origin. The typical size of Ar^em is 1%, for example in the 
case mt = ISOGeV, Ar ^ 3.8%, Xt ^ 1.0%, and Ar^em ~ 1-3%. 

Especially after the experimental observation of the top quark, we can safely argue that 
the fermionic sector of the SM is phenomenologically very well established. Fermionic loops 
are unambiguous theoretically, and the couplings have been tested in many high-energy pre- 
cision experiments. On the other hand, we know that, in addition to these fermionic con- 
tributions, there are a number of conceptually very important but numerically "subleading" 
corrections involving virtual bosons (Z, W^,f, and H). They involve the plethora of bosonic 
couplings of the SM, including the well-known tri-linear vertices, which have never been tested 
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Figure 3.2: Determination of sin^ Ow{Mz) in the full SM with Mh = lOOGeV. 

independently^, and affect self-energies, vertices, and box diagrams (in four fermion processes, 
bosonic vertex and box diagrams also include virtual fermions). Since the separation into 
fermionic and bosonic contribution is gauge invariant and finite at the one-loop level, and 
therefore unambiguous, it is possible to regard these two subsets as completely independent. 

With the increasing experimental accuracy the question arises of what is the sensitivity 
of present high-energy experiments to subleading contributions as the ones mentioned before, 
and to what extent present data provide evidence for the virtual structure of the theory. In this 
respect, the idea is to analyze the effect of the bosonic part of the theory on basic observables 
and their relations, and show that the current accuracy already allows us to recognize the 
presence of the associated virtual corrections. As an illustration, I will consider two different 
precise determinations of the MS parameter siii^ 6w(Mz) = i^, which conveniently describes 
physics at the Z-mass scale, as discussed in the previous chapter. Alternatively, this procedure 
can be viewed as a determination of sin^ ^^^^ in the MS framework. 

The first determination of sin^ dw{Mz) is obtained from the effective weak mixing angle 
sin^ ^^^^ which is now known with a very good precision. As mentioned above, the average of 
LEP and SLC asymmetry measurements gives sin^ ^^^^ = 0.2317±0.0004. The two parameters 
P and sin^ ^^^^ are related by Eq.(2.27) which, to very good accuracy, is independent of both 
rrit and M^- The parameter can also be obtained using Mz, G^, a through Eq.(2.19a), 



"2 "2 

s c 



V2G^MI{1 - Ar) 



(3.3) 



^The present experimental resolution on tri-boson couplings is very poor, but is expected to improve signif- 
icantly at LEP200. 
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Figure 3.3: Determination of sin^ Ow{Mz) in the full SM with Mh = 600GeV. 

and it has been shown that the radiative correction Ar is a sensitive function of rrit and M^- 
The results of the two determinations of P in the SM are shown in Figs. 3.1, 3.2, and 3.3, 
as a function of m^, in the case of respectively equal to 300, 100, and 600GeV, with the 
corresponding la errors. The error on the determination from Eq.(3.3) is almost completely 
due to the uncertainty on the light hadron contribution to the running of a. The comparison 
between the two determinations of siii^ 6w(Mz) forces a bound on the top quark mass. In 
the case of the SM with = 300 GeV this bound is in good agreement with the quoted fit 
of Ref. [3], and from Figs. 3.2 and 3.3 we can see that there is a slight preference for a light 
Higgs boson. In the evaluation of Ar I have neglected small O(d^) and 0{aas) corrections, 
except for the leading two-loop effects of 0{asGij,mf), which I have incorporated using a simple 
method recently discussed by Sirlin [51] and illustrated in Chap. 4. The values of the inputs 
are as reported in Chap. 2. 

In Fig. 3.4 the situation is analyzed in a truncated version of the theory, where all contribu- 
tions involving virtual bosons (vertices, boxes, and bosonic self-energies) have been removed. 
From Sec. 2.4, removing the bosonic contributions, one obtains Refc£(M|) = 1.0060, so that 

(i^)^^ = 0.2303 ± .0004 {Asymmetries). (3.4) 

Henceforth the subscript tr reminds us that this value corresponds to a "truncated" version 
of theory, with bosonic contributions removed in the electroweak corrections. In evaluating 
Refc£(M|) we have neglected small 0{das) corrections that were retained in Sec. 2.4. On the 
other hand, removing all the bosonic contributions from Ar leads to 
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Figure 3.4: Determination of sin^ dw{Mz) as a function of rrit from asymmetries and from Mz 
without bosonic contributions. 

Here I have also removed the irreducible two-loop effects of 0{a'^mf) [115] as they involve 
virtual bosons. It is apparent from Fig. 3.4 that the removal of the bosonic component of the 
electroweak corrections leads to a sharp disagreement. The two determinations do not overlap 
for any value of rrit in the present experimental range rrit = 180 ± 12GeV [1,2]. At the very 
conservative lower bound of = 150 GeV, the value obtained from Eq.(3.5) is 

(i^)t^ = 0.2275 ± 0.0003 {G^, a, Mz,mt = IbOGeY) (3.6) 

and we see that the discrepancy with Eq.(3.4) amounts to 5.6a (if the SLC value were not 
included, the top curve in Fig. 3. 4 would be shifted upwards by ~ 0.0004, the error would be 
slightly increased, and the difference would be 4.5a). As shown in Fig. 3.4, the discrepancy 
rapidly increases with m^. For instance, it is more than 7a for rrit = 180 GeV. Comparing 
Figs. 3.4 and 3.1, we see that the removal of the bosonic corrections leads to lower values of 
P. The discrepancy arises because the effect is much more pronounced in the [Gfj,, a, Mz) 
determination. An analogous investigation of the effects of bosonic contributions on the de- 
termination of Myy gives much less interesting results; the prediction derived using Eq.(2.12) 
is in acceptable agreement with the experimental value over most of the accepted rrit range, 
and no glaring contradiction is uncovered. Somewhat similar conclusions have been drawn 
in Ref. [52], where the effect of bosonic contributions was studied in s'm^ O^'^^p M^, and the 
decay width F^. 

We conclude that the accuracy currently reached in high-energy experiments is such that 
there is strong indirect evidence for the presence in the SM of these important subleading 
corrections. In the next few years precision experiments are likely to become sensitive to the 
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Figure 3.5: Determination of sin^ dw{Mz) as a function of rrit from asymmetries and Mz without 
Higgs contributions. 

structure of the gauge and spontaneous symmetry breaking sectors of the SM. The whole set 
of one-loop radiative corrections of the SM is definitely needed to describe current high-energy 
data. 

3.2.2 The Higgs contribution 

Given the sharpness of the signal, it is natural to ask whether one can use the same approach to 
probe specific components of the bosonic corrections. For instance, can one search for signals 
of Higgs boson contribution by removing them from the electroweak corrections, retaining the 
rest? As H does not contribute, at one-loop level, to ki^M^), the value of P derived in this 
case from sin^ ^^^^ is that of the full SM: 

= 0.2314 ± 0.0004 {Asymmetries), (3.7) 

instead of Eq.(3.4). In order to be physically meaningful, the removal of the Higgs contribution 
from Ar must be done in a gauge-invariant and finite manner. Fortunately, in the SM the 
diagrams involving H in the self-energies contributing to Ar form a gauge-invariant subset. 
On the other hand, they are divergent. Therefore, one must specify the renormalization 
prescription and the scale at which these partial contributions are evaluated. As P is the 
MS parameter and the electroweak data are dominated by information at the Z°-peak, it is 
natural to subtract the MS-renormalized Higgs boson contribution evaluated at the Mz scale. 
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Neglecting non- leading 0{a ) terms, the latter is given by 

« ( 1 S^n^-c^lnc^ 19 , 

' MS 



where ^ = M^/M^, H[^) is a function studied in Ref. [9,37], and the subscript MS reminds 
us that the MS renormalization has been carried out and the scale /x = Mz chosen. The need 
to specify the scale in defining the Higgs boson contribution can be most easily understood 
in the on-shell method of renormalization, where one employs sin^ 6w = 1 — M^/M^ instead 
of i^. In that case, the relevant radiative correction is Ar, rather than Ar. Although Ar 
is a physical observable and is therefore /x- independent, the Higgs-boson contribution is /x- 
dependent. Thus, a specification of the scale is necessary in its definition. 

Subtracting then the Higgs boson contribution one obtains a new truncated version of Ar, 
independent of M^, from which we can compute the corresponding {P)tr via Eq.(3.5). The 
comparison with Eq.(3.3) is given in Fig. 3.5. In contrast with Fig. 3. 4, where the complete 
block of bosonic contributions was removed, there are no signals of inconsistency. This is easily 
understood by noting that {Ar)HB. vanishes for ~ 113 GeV. Thus, the subtraction of 
Eq.(3.8) is equivalent to a SM model calculation with a relatively light Higgs scalar, ~ 113 
GeV, and this is consistent with current electroweak data. Similar results, from a diff"erent 
viewpoint, have been found in Ref. [53]. 

In summary, I have presented strong indirect evidence for the presence in the SM of 
bosonic electroweak corrections (Fig. 3. 4). If one probes just the Higgs component of these 
corrections, no evidence has been uncovered in our very simple analysis. However, it is likely 
that the signals will become sharper as the precision increases and rrit is measured. Finally, 
the approach I have illustrated, namely the removal of the bosonic electroweak corrections 
and their components from the relevant corrections, could be also extended systematically to 
global analyses. 



Chapter 4 

QCD corrections 



What is common to all men of goodwill is this: that in the final analysis our work disappoint us, 
that we must always start again at the beginning, that the sacrifice must always be made afresh. 

H. Hesse, Narcissus and Goldmund 

In the last few years many interesting results have been obtained in the field of perturbative 
QCD corrections to the electroweak phenomena, and the inclusion of leading two and even 
three-loop eff"ects has become routine in most phenomenological applications. Indeed, in most 
cases relevant for precision physics the eff"ective scale of the processes involved is high enough 
for perturbation theory to be reliable. A notable exception is the hadronic contribution to 
Aa discussed in Ch. 2. 

In this chapter I will first illustrate two examples of two-loop calculations, concerning the 
QCD corrections to electroweak bosonic two-point functions. The idea behind these calcula- 
tions is to provide a complete set of exact and compact 0{aas) expressions for electroweak 
bosonic vacuum polarization functions (VPF), in the most general possible case, including 
some extensions of the SM. From this set of general expressions it is possible to obtain asymp- 
totic expansions for some special case of physical interest (very large and vanishing external 
momentum and/or masses), as well as the imaginary parts, which in turn are related to the 
gluonic corrections to the hadronic decay widths of the corresponding boson. These calcula- 
tions generalize and extend previous work on the subject [54-57], which was limited to the 
transverse parts of the vector boson self-energies, and relied on some approximation on the 
masses involved, destined to cover only the present most important phenomenological cases. 
The final goal of this kind of eff"orts should be to place analytical calculations in the Stan- 
dard Model at a level comparable to what is known in QED where, in a pioneering work, the 
authors of Ref. [58] have derived exactly the VPF of the photon at two-loop order. 

As we will see in Sec. 4.6, the universal part of the QCD corrections to four-fermion 
electroweak processes (and the whole QCD correction in the case of no external hadrons) 
comes from the transverse parts of the electroweak vector boson self-energies. In particular 
this is the case of the most important precision observables, based on the muon decay and the 
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leptonic asymmetries. The longitudinal parts of the vector boson VPF's can however prove 
to be important, for example, in problems involving the production of heavy particles, and 
the scalar VPF is necessary in the case of Higgs decay, and to study radiative corrections to 
the masses of the Higgs particles involved in the symmetry breaking, both in the SM and in 
many of its extensions. Furthermore, an exact calculation without any assumption on the 
quark masses might turn out to be mandatory in the case of a fourth generation of fermions, 
the existence of which is still allowed by present experimental data if the associated neutrino 
is heavy enough [3]. 

Several checks of the cumbersome calculations are possible and will be illustrated. A nice 
example is provided by the use of Ward identities, described in Sec. 4.3, while the comparison 
of the dimensional regularization method with a dispersive approach will be discussed in Sec. 
4.4. 

QCD calculations, because of the strong scale dependence of the coupling constant as(/x^), 
are very sensitive on the choice of mass scale /x to be considered. As this scale dependence 
is unavoidable in perturbative calculations, where the perturbative series is truncated, this 
choice is crucial; a brief discussion of the subject will be given in Sec. 4.5. 

Finally, in the last section of this chapter, I will briefly review the incorporation of these 
QCD perturbative efl"ects in the most important electroweak corrections, and illustrate the 
flrst and most natural applications of the results of the preceding sections. In App. A and 
B the integrals used in the calculation, some asymptotic formulae for special situations of 
physical interest, and the imaginary parts are reported for completeness. 

4.1 QCD corrections to electroweak vector boson self-energies: 
complete calculation 

In this section, I present exact and compact analytical expressions for the 0{aas) contributions 
of quark doublets to the vacuum polarization functions of the electroweak gauge bosons in 
the most general case: real and imaginary parts of both the transverse and longitudinal 
components, for difl"erent and non-zero quark masses and for arbitrary momentum transfer [59]. 

The contribution of a fermionic loop to the vacuum polarization tensor of a vector boson 
i, or to the mixing amplitude of two bosons i and j, denoted II^j^, can be written as 

nifU?') = -i/ d'xe''^'' < 0|T* [Ji{x)4\0)] |0 > (4.1) 

where T* is the covariant time ordered product and q the four-momentum transfer; J^^,Jl 
are fermionic currents coupled to the vector bosons i,j and constructed with spinor flelds 
whose corresponding masses are ma,m}). The deflnition of 11^^, corresponds to times the 
standard Feynman amplitude. The vacuum polarization tensor can then be decomposed into 
a transverse and a longitudinal part, 

= {9,. - '-f-) n'/(g^) + ^ n-(g^) (4.2) 
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and the two components can be directly extracted by contracting n^j^(g^) by the two projectors 
Qfj.i' ~ ipiivli'^ ^iid ipiiv I i"^ '■> in n = 4 — 2e dimensions one has 



= ^n-U?^) (4.3) 

To set the notation it seems convenient to discuss briefly the well-known one-loop expres- 
sions of the fermionic vector boson self-energies: the one-loop transverse and longitudinal 
components can be written as (with s = q^) 

^Il{s) = ^N, [ {v'v^ + a'a^) s 0+^(5) + {vV - a'a^) m„m6 li-^s) ] (4.4) 

with Nc the number of colors of the fermionic doublet and and a' the vector and axial- 
vector couplings of the gauge boson i to the fermions expressed in units of the proton charge 
e = -v/ixa. The vector and axial- vector components of the vacuum polarization function, with 
coupling constants factored out, are then simply given by 

n^,i(5) = S n+^(s) ± maTTlb n~^(s) (4.5) 

which exhibits the fact that n^'^(s) can be obtained from n^'^(s) by making the substitution 
nia —rria (or m;, —m,})) as expected from 75 reflection symmetry. 

At the one-loop level the vacuum polarization amplitude receives contributions only by the 
graph of Fig. (4.1). The regularized amplitude must be multiplied by /x^^ for each momentum 
integration, where /x is the 't Hooft renormalizationmass scale introduced to make the coupling 
constant dimensionless in n = 4 — 2e dimensions; I have also introduced an extra term (e'''/4x)^, 
where 7 is the Euler constant, to prevent uninteresting combinations of ln4x, 7..., in the 
flnal result. After contracting the amplitude with the projectors [g^'^ — q^q^ I q^) and <l^<f ■, 
the resulting integrals can be expressed in terms of the scalar one-loop integrals that are given 
in App. A.l for completeness. 

The expressions of 11^^(5) in the general case rUa 7^ rri}) ^ are then 

, , , 2 + 3a + 3/3 1/2 \, , 5 a + /3 (a-jif 

- ^(a - /3)3ln I + [3(1 + a + /3) - A] (lna;„ + Ina;;,), 

n+(^) = ^ - i(« + PtPa + />6) + a + + ^(a - 

+ ^(A-l)(a-/3)ln|+ ^aI(A - 1 - a - /3)(lna;„ + Ina;;,), 

n-(.) = n-(.) = i - P^±P^ + 2 + ^ln| + \\H^^^a + Inx,), (4.6) 
where we use the variables 
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Figure 4.1: Feynman diagrams for the contribution of a quark isodoublet to the vacuum polar- 
ization function of a gauge boson at the one-loop (a) and two-loop (b) levels. 
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X = l + 2a + 2/3 + {a- /3)\ (4.8) 

and (3, p^, and X}, are defined accordingly, exchanging nia^b- 

In all the previous expressions the momentum transfer has been defined to be in the 
space-like region, < 0. When continued to the physical region above the threshold for the 
production of two fermions, s > (rrZa + mbY, the vector boson self-energies acquire imaginary 
parts, which are related to the decay widths of the vector bosons into fermions. Adding an 
infinitesimal imaginary part —ie' to the squared fermion masses, the analytical continuation is 
consistently defined. The imaginary parts of the vector boson VPF's are discussed in Appendix 
B.2. 

The knowledge of the imaginary part of the polarization function, which can be calculated 
directly using Cutkosky rules [64], allows an alternative way for obtaining the real part TZeI[[s), 
which can be expressed as a dispersive integral of Iml[[s). The connection between the 
dispersive approach and the results which are derived using dimensional regularization, will 
be discussed in Sec. 4.4. 

The imaginary parts of Eq. (4.1) are proportional to the decay width of the corresponding 
vector boson into quarks. This also explains why the expressions of Eq. (4.6) are free of mass 
singularities in the limit of vanishing quark masses: the real parts can be derived by dispersive 
integration of their imaginary parts, which are guaranteed to be free of these singularities by 
a theorem by Kinoshita, and Lee and Nauenberg [61,62]. Clearly, this must hold true at any 
order in perturbation theory, and will provide a check of the two-loop expressions. 



4.1.1 Mass definition and scheme dependence 

The two-loop diagrams contributing to the VPF n'^t(?^) at ©(aaj are shown in Fig. 4. lb. In 
the 't Hooft-Feynman gauge of QCD, using the routing of momenta shown in the figure and 
following the notations introduced in the previous section, one can write the bare amplitude 
as 

nifU?^) =--aa.(16x^) f^l" kif. + el (4-9a) 

^""^^ ' bare 3 ^ ' \ ^TT ) J (2x)" J (2x)" i 1^^^ ^ 



with 



fJ,J/ 



(fci - (kl - ml) {kl - ml) [(fci - g)2 - ml] [{k^ - q)^ - ml] 
(^1 + ma)jf,{v^ - a^j5)ih- k + 'mb)j^{v^ - a^j^){^i + ma)jxiM + 'ma)j^ 



ih - k2)^ [kl - ml)^ {kl - ml) [(fci - g)2 - m^] 
+ [rua ^ mt) (4.9b) 

In addition to these amplitudes, the expansion of the bare masses and couplings in the one- 
loop amplitude provides the necessary counterterms. By virtue of a QED-like Ward identity, 
the vertex and fermion wave function counterterms cancel each other, so that the electroweak 
couplings do not get renormalized at ©(a^), and only quark mass renormalization has to be 
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included. The latter is obtained by considering the diagram shown in Fig. 4. 2, the amplitude 
of which reads in dimensional regularization 

where p is the four-momentum of the quark and m its bare mass. This expression can be 
decomposed into a piece proportional to ( ^ — m) which will enter the wave function renor- 
malization and another piece proportional to m which will give the mass counterterm. After 
integration over the loop momentum, the latter is given by 

„2 



X \ m"^ / e 1 — 2e 




(4.11) 



The mass counterterm depends now on the choice of the renormalization condition, i.e. on 
the definition of the quark mass. For instance, the fermion mass can be defined in a gauge 
invariant way as the real part of the complex pole position s = — iV^ of its propagator. 
In applications, however, it is usually more convenient to use the zero of the real part of the 
inverse propagator: m = mo + ReS(m). These two definitions diff"er by terms of 0(r|/m), 
which are completely negligible in view of the foreseen experimental precision on the mass of 
the top quark. This mass will be denoted simply by m, and it is generally referred to as the 
on-shell, pole, or physical mass; the mass counterterm is then given by 

Us m ( iJ?Y 1 - 2e/3 / , \ ^ ^ 

" = VTi^j TTiT (1+12' +■■■]• 

where the ellipses stand for higher order terms in e. Pole masses imply the possibility of 
defining asymptotic states of the corresponding fermion fields. Although this is perfectly 
reasonable for leptons, the confinement of quarks makes this concept much less suitable to 
quark masses. A clear indication that in QCD the pole mass is not a well-defined concept 
will be described in Sec. 4.5. However, in the case of heavy quarks, the definition of mass 
as the pole of the propagator maintains a validity. Indeed, the lifetime of heavy quarks, 
and in particular of the top quark, is much smaller than the time-scale involved in strong 
interaction (of the order of 1/Aqcd)- Equivalently, the top decay width Tt, of the order of 
l-2GeV, is much larger than Aqcd ~ 200 — 300MeV. Moreover, the pole mass appears to be 
the natural parameter in the experimental determination of the top mass through kinematic 
reconstruction currently carried out at Fermilab [1,2]. From a formal point of view, the pole 
mass is the mass parameter implicitly used in the evaluation of VPF's through dispersion 
relations. 

Another possibility is the use of the MS scheme in which the mass is defined by simply 
removing the divergent term in the expression of Sm(m^) (the related constants ln4x,7, . . . 
have already been absorbed in the normalization). The MS mass rn(/x) is related to the pole 
mass (considering only QCD eff"ects) by [65,66] 

m = m{m) 1 + 1^ + if ^ + 0(«?) , (4.13) 
ox X 
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p — k 



(a) 




k — q 




(b) 



Figure 4.2: Feynman diagrams for the one-loop quark self-energy diagram (a) and for the contri- 
bution of the mass counterterm to the VPF at two-loop level (b). 



where K = 10.96 for the top quark, 12.5 for the bottom, etc. A study of the relation between 
fermion pole mass and MS Yukawa coupling in the SM can be found in Ref. [63]. 

What is relevant here is that the definition of mass is completely arbitrary and there may 
be different legitimate choices, depending also on the problem at hand. Once the choice is 
made one inserts the counterterm in the one-loop self-energies, as depicted in the diagrams of 
Fig. 4. 2b, which is equivalent to calculate 
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-Sm, 
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m„ 



1— loop 



1— loop 



(4.14) 



where the one-loop vacuum polarization function is given by Eq. (4.6). The renormalized 
two-loop self-energies will then read 



nii(g^) = n^(g2) 



bare 



nii(g^) 



CT 



(4.15) 



Having at hand the result of the vacuum polarization function in the on-shell scheme, one 
can obtain the polarization function in any other scheme X: one simply has to add to the 
expression of the two-loop self-energy Eq. (4.15) the quantity^ 
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" d 



m„ 



nii(g^) 



1— loop 



d 



drrib 



nii(g^) 



1— loop 



(4.16) 



^Note that since at this stage we are only discussing the difference between two renormalization schemes, we 
do not need the 0(e) terms in the one-loop result of the vector boson self-energies. These terms are of course 
needed to evaluate the renormalized two-loop self-energies in a given scheme. 
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where Am^j^ = rUafi — ''^ab '^^ difference between the pole mass rUafi and the quark mass 
m-^, defined in the scheme X. 

a,o 

Decomposing the one-loop transverse and longitudinal components as in Eq. (4.4), one 
obtains for the derivative of these components 

d ct \ \ ct 

^^^::r nJ(5)|i_ioop = - - o«(/'a + Pb) + - a) - -a{a - pfln 



drua 



a(l - a + - 2l3\-^){lnxa + In Xb), 



n+(s)|,_i^^p = ^ - -a{pa + Pb) + a(2 - 3/3 + 3a) 



drUa 



^a[3(a-/3)2 +4a]ln^ 



a 1 

— A2 

2 



/3 

l + S{a-(3) + 2^{l + (3-a) 



(4.17) 



rua-z mamb\lj. j^{s) 

drUa 



11 1 r,\^ ^ 

mamb ---(/)„ + /)6) + 2 + -(3a - /3) In - 
i-ioop le z z p 

^A2 [l + 2aA-i(l + a - /?)] (lna;„ + In a;;,) 



and similarly for the piece involving rrib which can be obtained by making the substitution 
nia nib. Once the complete result of the VPF is known in the on-shell scheme, it is 
therefore straightforward to obtain the corresponding results in any renormalization scheme 
using Eqs. (4.16,4.17). In practice, the most useful cases are the ones with equal quark masses, 
corresponding to the Z boson, and with one nearly massless quark, corresponding to the [t — b) 
doublet contribution to the W boson VPF. The explicit expressions for these two cases are 
(up to a factor Ncu/tt t;^(a^)) 

A^n^ = 

'-'-T,ma=mi, 

A"'n^;:,=o = -mi =^ [--pa-a+{l-a')ln ^ ] (4.18) 

A few remarks are necessary at this stage. 

First, one notices in the previous expressions the occurrence of terms that are inversely 
proportional to the velocity factor A^/^. In principle, after an analytical continuation to the 
physical region beyond the threshold for quark pair production, s > [ma + rubY, these terms 
would diverge for energy values near the production threshold, A ~ 0. However, as we will 
see later, when evaluated in the on-shell scheme, which is the only scheme where the physical 
threshold is well defined, the renormalized two-loop polarization function is free of these A~^/^ 
factors. Near threshold, the dominant terms will be constant and would correspond, once the 
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vector boson self-energy is normalized to its one-loop value, to the well known A~^/^ Coulomb 
singularities which require a non-perturbative treatment [62]. Of course, the A~^/^ terms can 
be present in the vacuum polarization function when it is evaluated in a different scheme but 
in this case the threshold is not well defined since the masses are not "physical" masses. 

In principle one can even define the masses nia and m;, in two completely diff"erent schemes; 
this could be useful if, for instance, diff"erent scales are involved in the evaluation of the VPF's. 
In the top-bottom isodoublet, one could employ the on-shell mass for the heavy top quark, 
and use a running mass evaluated at the scale for the relatively light b quark, which in 
general avoids the appearance of large logarithms for ^ ml (see the example at the end of 
App. B). In the rest of this section, however, I will employ exclusively on-shell masses. 

4.1.2 Exact two-loop results 

In this section I give the expression of the vacuum polarization function at order 0{aas) in 
the general case nia ^ m,}) ^ and for arbitrary momentum transfer [59]. The result will be 
expressed in terms of on-shell quark masses. I will use the same notation previously introduced 
at the one-loop level, and nia^b will stand for the on-shell masses. 

At O(aas), the transverse and longitudinal components of the vacuum polarization func- 
tion n^j^(g^) can be written as (the color factor iVc = 3 is now included^) 

nt%(5) = [ {v'v^ + a'a^)sl[+^{s) + {v'v^ - a'a^) mam^ H" Js) 1 , (4.19) 

TT TT I ' 'J 

where 11^^ are the sum of the corresponding components in the bare two-loop amplitude 
Eqs. (4.9) and the mass counterterm which can be obtained from the preceding section. 

In analogy with the one-loop case, after contracting 11^^, in Eqs. (4.9) by the tensors 
{9^'^ — i^i^ li"^) and q^q^ Iq^ and expressing the scalar products of momenta appearing in 
the numerators in terms of combinations of the polynomials in the denominators, one is led 
to the calculation of a set of scalar two-loop integrals. Most of these integrals have been first 
calculated by Broadhurst in Ref. [133]; the remaining integrals reduce after straightforward 
computations to the previous ones [55,134] and are all listed in App. A.l. After a cumber- 
some computation, and taking advantage of the symmetry in the change a we can write 
11^^(5) in a relatively simple and compact form^ 

n+ = {-^ + ^ + 2ap^ - ^a) + iE^±P^ (11a _ 1 _ 9a/,„ + 2ap,) 

+ 24 " 24" " T"^ + 3«[G'(2^&) - G{xa)] 

+ ^ [(a - /3 + A2) (11 + 19a + 19/3 + 12a/3 - 5A) + (a - /3) 

X (42 - 5a - 5/3)] + | In In [3 + /3(4 + 2a)] 

^The case of SU(N) with N ^ Z can be obtained by multiplying our result by (iV^ - l)/8, while for QED 
it is sufficient to multiply by 1/4 for leptons, and by 3/4 QiQj for quarks, Qi being the charge of the i quark 
in unit of the positron charge. 

^Minor differences with the expressions in Ref. [59] are due to trivial algebraic simplifications. 
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2^ 
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3a 



4na;„ [a - /3 + (a + /3 - 3)(a - /3 + A^)] - (1 + 2a)2: - 21' 



In Xa 



[l + a + /3-X+{/3- a)A2] + ^ In In a;;, [1 + 2a - 2a/3] 



+ {a^(3}, 
with I and T' given by 



(4.20) 



I = F{1) + F{XaXb) - F{xa) - F{xb) 

I' = X^G{xaXb) - ^(/3 - a + X^)G{xa) - ^{a - /3 + X^)G{xb) (4.21) 

In terms of the polylogarithmic functions [135] Li2{x) = — Jq \n[l — xy)dy and Li3(a;) = 
Jq lni/ln(l — xy)dy, the functions F and G are given by 



F{x) = 6Li3{x) - 4Li2(a;)lna; - In^ a;ln(l - a;) 
G{x) = 2U2{x) + 21na;ln(l - a;) + In^ x 



(4.22) 
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or also 

F{x) = ^[(2 - nlnxf + 2] a;"/n^ 

n=l 

These two functions also admit a simple and useful integral representation [133] 

The imaginary parts of the VPF are derived along the same lines as in the one-loop case, 
and are given in App B.2. Incidentally, we note that the analytic structure of the VPF is 
the same as at one-loop level, with a unique threshold at (rrZa + rubY. This is obviously due 
to the masslessness of the gluon. In general, at two-loop level the presence of new virtual 
massive particles in the loops prompts new thresholds, and it is this feature that complicates 
the evaluation of two-loop massive integrals. For instance, it has been noted [67] that two-loop 
scalar integrals can be expressed in terms of polylogarithmic functions only in a very restricted 
set of cases including the one we are considering (see also App. A). 

The gluonic corrections to the hadronic widths of vector bosons in the case of arbitrary 
quark masses and for both transverse and longitudinal components have been known for some 
time [56,68,69] (in the first reference only the transverse part is given) and were obtained by 
directly calculating the QCD corrections to the flavor changing decay of a vector boson. These 
calculations were mostly motivated, in the context of QCD sum-rules, by the need to calculate 
at flrst order in perturbative QCD the Wilson coefficients of dimension two operators. The 
results presented here were obtained using a completely different method, and agree with those 
of Ref. [56,68], and with Ref. [69], if some misprints in the integrals of their Appendix (Ji 
and J2) are corrected (see also Ref. [70]). This is a very important check of our calculation. 

In App. B the asymptotic expressions for the most important special cases, corresponding 
to two equal mass quarks, one massless quark, and ^ ^ are presented in detail. They 
represent excellent approximations of the case in which we are going to be mostly interested: 
the [t — b) doublet, where the mass of the bottom can be safely set to zero as in most 
applications the neglected terms will be of 0(m^/g^) = 0(m^/m^). The case of vanishing 
external momentum is discussed in the next section. We also verify that in the limit of 
vanishing quark mass, the self-energies are free of mass singularities as required by the mass 
singularity theorem [61,62]. 

Even though the real parts of Eq. (4.20) are a new result, a number of partial results allow 
us several useful checks of the calculation. In the equal mass case, 11^(5) has been derived in 
the 1950s by the authors of Ref. [58] by means of a dispersive approach, and, more recently, 
IIt'^{s) have been calculated in Ref. [56] by the same method. Real parts of the transverse 
components in the cases m;, = and nia = m,}) had been previously derived using dimensional 
regularization in Ref. [54,55] ^. Finally, the expressions of Il^'^(s) in these special cases has 
also been obtained in Ref. [57], again by means of dispersive integration. 

^As explained in Ref. [59], the expressions of Ref. [54,55] differ from the ones reported in App. B by a 
constant term proportional to tt^ , due to an omission in the mass counterterm used in that calculation. This 
does not affect most physical observables. 



CHAPTER 4. QCD CORRECTIONS 



38 



4.1.3 Ap and a fourth generation of quarks 

An immediate application of the expressions of the previous section is the calculation of the 
effects of a heavy quark isodoublet with SM couplings on the p parameter, which measures 
the ratio of neutral and charged couplings at zero momentum transfer. 

The fermionic one-loop contribution to the p parameter can be expressed in terms of the 
difference between the transverse components of the Z and W bosons self-energies at zero- 
momentum transfer, = 



Ap 
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ng^(o) 

Ml 
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(4.24) 



The contribution of a quark pair with different masses, nia ^ rrib ^ 0, to 11^,'^ (0) in the 
one-loop approximation is 
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(4.25) 



where the upper (lower) sign refers to vector (axial) current. Note that, in this limit, the 
longitudinal and transverse components of the self-energies are equal. This is a trivial con- 
sequence of the analyticity of n^'^(g^) at = 0; in particular, the coefficient of q^q'^ has to 
be regular in the limit q^ —>■ 0, thus forcing the identity nT'^(O) = n^'^(O). Similarly, the 
conservation of equal mass vector current prompts a Ward identity which implies 11^(0) = 
for ma = mi). Ap can therefore be written in terms of either the longitudinal or the transverse 
components of the VPF: 



with 



fL,T{s,ma,mb) 



E E 
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(4.26) 



ni,T(s,mi,mi) 



III^t{s, ma, mb) - n^,T(s, rua, mb), 



(4.27) 



where I have kept some n^(s,m, m) terms even if they do not contribute to Ap for future 
convenience. Using the previous expressions for ll^'^^(0) one obtains the one-loop correction, 
which is quadratic in the fermion masses [71] 



f^'\0,ma,mb) 



ml 



2mlm? , m? 
— —In — - 

ml - ml ml 



(4.28) 



At O(aas), the last equation gets modified in the case of virtual quarks. In the limit of zero 
momentum transfer, the expressions of 11^'^ (-s) simplify considerably, and one obtains from 
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Eqs. (4.20) (n^'^(o) = nl'^(o) = n^'^(O)) 




4 ml -ml {ml - mlf ml 2 



(4.29) 



where the function G{x) has been introduced in the previous section, and again the upper 
(lower) sign refers to vector (axial) current. 

In the two special cases ma = m^ and m^ = the last equation simplifies further to 



As expected, there is no singularity in the self-energies in this limit. In addition, besides the 
manifest symmetry in the exchange ma ^ m^, the previous formulas exhibit the fact that 
n^'^(O) can be obtained from n^'^(O) by simply making the substitution ma —ma (or 
mi) —mb) as expected from 75 refiection symmetry. 

Using the previous expression, one readily obtains the QCD corrections to the contribution 
of a heavy quark isodoublet to the p parameter in the general case ma 7^ m^ ^ 0. Defining 
the contribution to the p parameter in analogy with Eq. (4.26) 



"6 1 

ll^{0,ma = mb) = ml ^ + -(11 
n^'^(0,m, = 0) = ^[^ + ^(ll 



12pa) + I2pl - 22pa + Y + 
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(4.31) 



the function f^'^\ defined as the two-loop analogue of f^^\ 



is given for iVc = 3 by 



/(^)(o, 




= - <ma + mb+ 2 — ^ In — 1 + 

[ ml -ml ml 




+ {ml-ml) Uu, f4U21n4lnfl 




(4.32) 



Notice that /(^^(O, rrZa, m;,) is free of ultraviolet divergences as it should be since Ap is an 
observable physical quantity, and that it does not depend on the 't Hooft mass scale /x. Note 



CHAPTER 4. QCD CORRECTIONS 



40 



0.05 - 



0.04 - 



0.01 - 



0.06 - 



0.03 - 



0.02 - 



- 








100 



200 



300 



400 



500 



600 



700 



Figure 4.3: Contribution of a fourth generation of quarks to Ap, in the case of one quark mass 
equal to SOOGeV and = 0.11, as a function of the other quark mass (GeV). The full line 
represents the one-loop contribution; the dashed line incorporates also two-loop 0{aas) terms. 

also that the cofactor of ml - ml in Eq.(4.32) is simply {G{ml/ml) - G{ml/ml))/2, which 
makes the symmetry in the interchange of ma and m^ explicit (we have used the identity 
G{l/x) + G{x) = 2x^/3). In the limit of large mass splitting between the two quarks, ma ^ mf,, 
the QCD corrections to the p parameter reduce to the known result for m^ = [54,55] 



where it is apparent that the two-loop QCD correction significantly "screens" the one-loop 
effect (for = 0.11, which roughly corresponds to at the top mass scale, the QCD 
correction amounts to about 10%). 

The contribution of a new quark isodoublet to the p parameter is exemplified in Fig. 4.3 
both in the one-loop approximation (full line) and including the QCD radiative corrections 
(dashed line). I have set the mass of one quark to 300GeV and varied the mass of the 
other quark from to 700 GeV; for the strong coupling constant I have used the numerical 
value Us — 0.11 (see Sect. 4.4 for a discussion of the choice of scale). As one can see, the 
contribution of the heavy quark isodoublet to Ap vanishes for degenerate quarks, ma = mi. 
This is expected, because the deviation of the p from 1 is due to a violation of the SU(2) 
custodial symmetry [72] which can be induced only by hypercharge or by isospin violation in 
the Yukawa sector. In our case, only the latter actually contributes, and we see that the QCD 
corrections are always negative, thus screening the one-loop contribution to Ap. In fact, in 
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Figure 4.4: QCD correction to the contribution of a fourth generation of quarks to Ap, normalized 
to the one-loop contribution. One quark mass is set to SOOGeV and the other one is varied from 
to 700GeV; = 0.11. The dashed line represents the same correction when the one-loop 
contribution is expressed in terms of MS masses. 



this mass range the QCD corrections are practically constant and amount to a decrease of Ap 
by approximately 7-8%, as can be seen from Fig. 4. 4. In the case of nearly degenerate quarks, 
Eqs. (4.31,4.32) reduce to 
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(4.34) 



It is interesting to see what happens if one uses a different definition of mass: for example, 
if we express the one-loop correction in terms of MS masses rhafi(rnafi), using Eq.(4.13) and 
neglecting higher orders in a^, we obtain an extra term in 0{aas) equal to |^ times the 
one-loop correction. In the case of one vanishing mass, Eq.(4.33) is modified into 
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(4.35) 



In this way the QCD correction to Ap is much smaller, about 0.7%; the comparison with the 
on-shell case is shown in Fig. 4. 4. 

For what concerns the top-bottom doublet, I illustrate the effect of keeping rrib ^ for 
rrit = 200GeV, and using = 0.11. From Eq. (4.33) we see that the one-loop contribution 
in the approximation of neglecting the mass of the bottom amounts to 1.2535%, while using 
Eqs. (4.26-4.28) with rub = 5GeV we obtain 1.2427%. The ©(a^) terms in Eq.(4.33) and 
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Figure 4.5: Contribution of the {t, b) doublet to Ap for 150 < rut < 220GeV and = 0.11. The 
one-loop contribution with m;, = 5GeV (full line), and the sum of one and two-loop contribution 
for m,}) = (dotted line) and for m;, = 5GeV (dashed line) are shown. 

(4.31-4.32) add to these results -0.1255% and -0.1225%, respectively. We conclude that the 
finite bottom mass effect on the one-loop Ap amounts to about 10""*, while the effect on the 
QCD corrections to Ap is a mere 3xl0~^. In other words, a 5GeV bottom mass decreases 
the QCD corrections to the one-loop formula by about 2.4% for rrit = 200GeV, roughly of the 

2 2 

order rrib/mt, as an effect of terms like ^In^ These contributions are well beyond the 
accuracy of current analyses. The finite bottom mass effect is illustrated in Fig. 4. 5 , in the 
range 150GeV< rrit < 220GeV. 
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4.2 QCD corrections to the Higgs vacuum polarization func- 
tions 

The existence of at least one scalar particle, the Higgs boson, is required in the SM to generate 
the masses of the other fundamental particles, leptons, quarks and weak gauge bosons [75]. 
The discovery of this particle and the study of its fundamental properties will be the most 
important mission of future high-energy colliders. 

The phenomenological properties of the unique Higgs particle of the SM have been stud- 
ied in great details in the literature [75,76]. In fact, because the precise knowledge of Higgs 
decay widths, branching fractions and production cross sections is mandatory, quantum cor- 
rections must be included and this subject has received much attention recently [77]. In 
particular, QCD corrections to Higgs decay and production processes are of the utmost im- 
portance. For instance, the Higgs decays into quark pairs and gluons, which together with 
the H T~^T~ decays are the most important decay modes in the intermediate mass range, 
M^^Mjf^ 140GeV, receive very large QCD corrections. In the case of H ^ qq, they are 
known exactly to 0{as) [78,79,68] and up to 0{a1) [80] in the approximation rUq <C Mh', 
in the case of the gluonic decay, H ^ gg, the QCD contributions are known up to next-to- 
leading-order [81], and the leading electroweak correction have also been calculated [82]. 

The electroweak radiative corrections to the SM Higgs decays [83] are of significance as 
well, as the leading contribution is proportional to the squared mass of the heavy top quark. 
In fact, a fourth generation of heavy fermions, the existence of which is still allowed by present 
experimental data with the proviso that the associated neutrino is heavy enough, would have 
a dramatic eff"ect on the Higgs decay widths. Its contribution is universal in the sense that it 
does not depend on the final state particle, and will also increase quadratically with the heavy 
fermion masses. The universal part of the two-loop mixed ©(asG^mg) corrections, which 
have been calculated very recently [84,85], will screen the leading one-loop contribution by a 
non negligible amount. 

Many extensions of the Standard Model predict the existence of a larger Higgs sector. 
For instance, supersymmetric theories (SUSY) require the existence of at least two isodoublet 
scalar fields $i and $2 to give masses separately to isospin up and down particles, thus 
extending the physical spectrum of scalar particles to five [75]. The physical Higgs bosons 
introduced by the minimal supersymmetric extension of the Standard Model (MSSM) are 
of the following type: two CP-even neutral bosons h and H (where h will be the lightest 
particle), a CP-odd neutral boson A (usually called pseudoscalar), and two charged Higgs 
bosons H^. Besides the four masses Mh, Mjj, Ma and Mjj±, two additional parameters 
define the properties of the scalar particles and their interactions with gauge bosons and 
fermions: the ratio of the two vacuum expectation values tan/? = v^jvi and a mixing angle 
a in the neutral CP-even sector. Note that, unlike a general two-Higgs doublet model where 
the six parameters are free, supersymmetry leads to several relations among these parameters 
and only two of them are independent. These natural relations also receive large radiative 
corrections [86]. 
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In this section I present a calculation [74] of the hadronic contributions to the Higgs boson 
self-energies at 0{aas). The most general case has been considered, leaving the momentum 
transfer arbitrary, with different quark flavors u ^ d and consequently difl"erent quark masses 
rriu 7^ md- The motivation for performing such a calculation is threefold: (i) As in the case 
of the SM Higgs boson, the strong [78,79,68,87] and some of the electroweak [88] radiative 
corrections to SUSY Higgs boson decays are known at the one-loop level. In some limiting 
cases, as for the gluonic corrections for nearly massless quarks, SM two-loop results can be 
adapted to the SUSY case. Here, we provide the necessary material which allows to derive 
the universal part of the mixed 0{asGF) radiative corrections to these Higgs decays. This 
is a generalization to a multi-Higgs doublet model of the recent SM calculation [84], which is 
just a special case (for a CP-even neutral Higgs boson) of the result presented here, (ii) The 
imaginary parts of the Higgs boson self-energies are related, through the optical theorem, to the 
partial decay widths of the Higgs bosons into quarks. Here I give exact analytical expressions 
of the QCD corrections to the Higgs decay widths in the most general case ^ m^, which 
is not available in the literature^. In the special cases = m^, the known results for the 
QCD corrections to CP-even and CP-odd [78,79,68] neutral Higgs bosons can be recovered. 
Since the expressions reported here have been obtained by a completely difl"erent method, this 
serves as an independent check, (iii) We can derive Ward identities that relate scalar VPF's 
to the longitudinal components of vector boson self-energies. In Sec. 4.3 I will show explicitly 
that it is possible to use these Ward identities to provide a powerful consistency check of the 
calculations described in the previous and in the present section. 

4.2.1 Notation and one-loop results 

In order to set the notation, I rederive here some one-loop results which will be relevant to 
our next discussion. 

The contribution of a quark loop to the self-energy of a scalar Higgs boson $ will be denoted 
by H*(s = q^) where q is the four-momentum transfer, and will correspond to —i times the 
standard Feynman amplitude. To treat the cases of neutral CP-even, neutral CP-odd and 
charged Higgs bosons on the same footing, it is convenient to work in the general situation 
where the internal quarks in the loop are of difl"erent flavor, and have difl"erent masses. This 
corresponds to the case of a charged Higgs boson which couples to an up-type and down-type 
quark, with arbitrary masses mu^d] the self-energies of neutral scalar and pseudoscalar Higgs 
bosons will be special cases of the previous one. 

The coupling of charged Higgs bosons to fermions is a P-violating mixture of scalar and 
pseudoscalar couplings 

g{H+ud) = i(GF/V2y^^[K{l-'f5) + Ml + 75)]. (4.36) 

In the so-called type II two-Higgs doublet model [75], in which one doublet couples to the up 
quarks and neutrinos while the second doublet couples to down quarks and charged leptons, 

^QCD corrections to the charged Higgs decay for ^ md have been calculated in Ref. [87], but not in a 
full analytical form as the integrals for the real corrections have been performed numerically. 
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hu/rriu 


hd/rrid 


HsM 


1 


1 


h 


cos a/ sin/3 


— sin a/ cos 


H 


sin a 1 sin (3 


cos a/ cos (3 


A 


1 /tan/3 


tan/3 



Table 4.1: Neutral Higgs couplings to up-type and down-type fermions in the SM and MSSM. 



we have 



hu = mu/tan/3, = tan/3. (4.37) 

The MSSM belongs to this type. It is often convenient to use the scalar and pseudoscalar 
components of this coupling 

V = hd + K, a= hd- K- (4.38) 

The couplings of scalar, that we will denote by S, and pseudoscalar A Higgs bosons take the 
general form 

g{Sqq) = -i (^^.72)'^' /i„ g{Aqq) = - (^^.72) 75 K (4-39) 

In the SM, the reduced couplings hq are just the quark masses; in the the MSSM these 
couplings for ^ = S,A are given in Table 4.2.1 for u and d quarks. 

In the one-loop approximation, the contribution of a quark loop to the vacuum polarization 
amplitude of a charged Higgs boson, H'^(g^), corresponds to the diagram of Fig. 4.2.1a. For 
arbitrary fermion masses 7^ rrid 7^ and momentum transfer g^, this amplitude can be 
written as 



n^(,) = 



(4.40) 



^2 

where, using a = —m^/q'^, (3 = — m^/g^, pafi = In-^^, and all the other variables introduced 
in the last section, H^(s) take the following form (H+ = H^ = H"*") 

n+ W = + + a + ffl - + 1 + |(a + ffl 

+ ^(l + a + /3) [(a-/3 + A2)lna;„ + (/3-a + A2)lna;6] , 
n-{s) = ---2 + _ 1 [(a _ /3 + Xl)hix^ + (/3 - a + X^)hiXb] . (4.41) 
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Figure 4.6: Feynman diagrams for the contribution of a quark isodoublet to the self-energy of a 
Higgs boson at the one-loop (a) and two-loop (b) levels. 
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As already done in the previous section, ln4x and 7 factors have been absorbed in the 
normalization. The expressions of the self-energies for neutral scalar and pseudoscalar Higgs 
bosons can be simply obtained by setting = = niq in Eqs. (4.40) and (4.41), and by 
using the relevant couplings which are given in Tab. 4.2.1. The result is 



n^'^{s) = ^^^shl [n+(.) ± {ml/s)Il-{s)] 



(4.42) 



leading to 



n-(.) 
n-(.) 



— ;= — s n„ 
2^2x2 ^ 

SGf ,2 
— ;= — s n„ 
2^2x2 ^ 



^ - Pa^ (1 + 6a) + 1 + 5a + ^(1 + 4a)2 Ina; , 
^ 0-/'a)(l + 2a) + l + a + ^(l + 4a)2lna; , 



(4.43) 



where x = 4a/(l + ^/TT~^Y^ 

The imaginary parts of the scalar VPF's can be obtained along the same lines as described 
already for vector boson self-energies. The analytic structure of the VPF's is also identical. 
Considering the normalization we have chosen, the imaginary parts of a charged scalar self- 
energy are related to the partial decay widths of a charged Higgs boson into ud quark 
pairs through 



T{H+ ud) 



N,GfMh+ 

2^2x2 

rriumd 
2huhd — 



[hllrnU.+ {s) + hllrnU+{s) 
IrriUri^s) , 



(4.44) 



and the case of neutral Higgs can be treated in a similar way. 



4.2.2 Exact two-loop results 

At O(aas), the two-loop diagrams contributing to the Higgs boson self-energies H*(g2) are 
shown in Fig. 4. 2.1. In the 't Hooft-Feynman gauge, using the routing of momenta shown in 
the figure, and using^ Nc = 3, one can write the bare amplitude as 



H*(g2) 



WttGf 
3V2 



Nra, 



4x 



j J (2x)" j 



d^k2 



[A + B] 



(4.45) 



where 



A=Tt 



(^1 + mu){v - a75)(^i- + md)7A(^2- ^ + md){v + a75)(^2 + "iu)7^ 



,B=Tr 



{h - k^Y {kl - ml) {kl - ml) [{k, - qf - ml] [{k^ - qf - ml] 
{jki + mu){v - a75)(^i- ^ + md){v + a75)(^i + m„)7A(^2 + "iu)7^ 



(fci - ^2)2 {kl - mlY (k\ - ml) [(fci - g)2 - ^2] 
md- 



(4.46) 



^The case of different color group can be obtained as described in Sec. 4.1.2. 
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This bare amplitude has to be supplemented by counterterms; we can proceed in analogy 
with the case of vector boson VPF, with the difference that the Yukawa couplings of the 
quarks with the scalars are proportional to the quark masses, and they must be renormalized 
according to the mass definition. The result is that again only quark mass renormalization has 
to be considered, and the diagrams with quark mass insertions into the one-loop amplitude 
are the analogues of the ones in Fig. 4. 2b, with the addition of appropriate mass insertions in 
the vertices^. This is of course equivalent to calculate 



CT 



orriu 



1— loop 



8mdd 
drud 



1— loop 



(4.47) 



where the one-loop VPF is given by Eqs. (4.40, 4.41) and h^^d have been expressed in terms 
rriu^d before diff"erentiation; 0(e) terms not present in Eqs. (4.40,4.41) have to be included. 
The two-loop self-energies expressed in terms of renormalized masses then read 



n*(g2) = n*(g2) 



n*(g^) 



CT 



(4.48) 



The calculation proceeds exactly in the same way as for the vector boson VPF's, and involves 
the same integrals (see App. A.l). 

Here I give the expression of the contribution of a (u, d) isodoublet to the charged Higgs 
boson two-point function at order O^aa,) in the case of arbitrary quark masses [74]. The 
result is expressed in terms of on-shell quark masses rriu^d- Using the same notations as in the 
one-loop case, the charged Higgs boson self-energy H'^(g^) at the two-loop level is given by 



H^(.) 



Gf as 
2^2x2 X 



hlJl+{s) + hlJl+,{s) + 2Khd 



ruumd 



H-(.) 



(4.49) 



with H=^(s) given by the relatively simple and compact expressions 



H^ 



^ + 14(a + 13)- 3pa - Uapa - 6/3(/>„ + pb) 



{Pa + Pb) 



4a + 4/3) - - 
^ + 14(a + /3)-3(/>„ + />,) Q 



a + /3 -3(a-/3)ln 



1 " 



5 + f (<» + /?) - f (1 + 20 + 2/3)(p. + p,) + |(a' - f')^J 

53 3 1 a 

y (a + /3) + -A2 (1 + a + /3)(ln a;„ + In a;;,) In - 



^A2 (1 + a + I3){\uxa + Inxb) + In | (40 + 9(a + /?)) 

In^ Xa + In^ Xb 



[(1 - 2a - 2/3) (1 - A + a + /3) - 3(a + /3) + 3(a - pf] 



'^Having at hand the expressions of the two-loop self-energies in the on-shell scheme, the procedure for 
obtaining the corresponding result in any other renormalization scheme is straightforward and completely 
analogous to the one described in Sec. 4.1.1. In the following I will express all results in terms of on-shell 



'Minor differences with Ref. [74] are due to trivial algebraic simplifications. 
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+ — In - (lna;„ + In Xb) {I + a + /3) - —(1 + 4a + 4/3) 

+ ^ In In Xb{\ + 2a + 2/3 + 6a/3) - (1 + a + l3fl - 2(1 + a + l3)l', 

n+ = n+ [ m„ ^ rrzd ], (4.50) 

6 1 

n- = — + -(14 - 6pa - 6pb) - 14(/9a + Pb) + ^Pa + Pbf + 20 + 

- 61na;„(a-/3 + A2)-ln2a;„ [a - 1 - a - /3 + (a - /3)A2] 

- 61na;6(/3 - a + A2 ) - In^ a;;, [a - 1 - a - /3 + (/3 - a)A2] 

- 61na;„lna;6(l + a + /3) + 2(1 + a + /3)2: + AT', 

where the functions F and G and the integrals I, I' are given in Sec. 4.1.2. One can derive 
the expressions of the self-energies for neutral scalar and pseudoscalar Higgs bosons from 
Eqs. (4.50) by setting = = rriq and using the proper couplings. This is done in 
App. B.l, where the corresponding formulae are shown, together with the limits for infinite 
and vanishing momentum transfer. The expression of 11^ recently derived in Ref. [84] in the 
case of the SM Higgs boson is in agreement with our result. Explicit 0{Gfj,asmf) formulae for 
ng'^(g^) are crucial in order to calculate the leading radiative corrections to the Higgs boson 
fermionic decay widths. In the case of leptonic decays, for instance, there is only a universal 
factor coming from the Higgs wave function renormalization and the renormalization of the 
Higgs vacuum expectation value^. 

Analytic results for TmH^(s) [68,69,78,79] and TmH^(s) [79,68] have been obtained in 
the past by a number of authors by directly calculating the QCD corrections to the decay 
of a scalar Higgs boson into quark pairs. The results that we obtain from Eq. (4.50) using 
a completely diff"erent method agree with the previous ones; this serves as a check of our 
calculation in the general case. The expressions of the imaginary parts in those cases and for 
arbitrary quark masses are reported in App. B.2. Note also that for the value tan/3 = 1, we 
recover the expression of the imaginary part for the longitudinal component of the electroweak 
vector bosons in the general case, which is given in Refs. [68,69]. Indeed, because of a Ward 
identity discussed in the next section, the imaginary part of the longitudinal component of 
the vector boson self-energy is the same as the one for the Higgs boson self-energy for this 
value of tan/3. This feature provides also a powerful check of the calculation presented here. 

4.3 Ward Identities: a consistency check 

It is well-known that there are Ward identities relating the longitudinal components of the 
electroweak vector bosons and the corresponding Goldstone bosons [89] (see also Sec 5.3). 

^This factor has been calculated independently in Ref. [74,84]. For a more detailed discussion see [74] and 
references therein. 
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In this section, I will use the current algebra of the Standard Model to derive these Ward 
identities; I will then briefly show how to relate the Higgs boson self-energies calculated in the 
preceding section to the longitudinal parts of the electroweak vector boson self-energies, given 
in Sec. 4.1.2. 

Deflning the fermionic contribution to the VPF of the W boson and of the corresponding 
Goldstone boson $ as [g is the SU(2)i coupling constant) 

nrw(g') = -iy / d-xe-''^-^0\T*jl^{x)rUO)\0) (4.51) 

n$$(g') = I d-xe-''i-^0\T*S\x)S{0)\0) (4.52) 

where Jw{x) and S{x) are the charged fermionic currents coupled to the W and to the $ bosons 
and T* denotes the covariant time ordering product; for the notation and normalization of 
the currents, I follow Ref. [44] (see also Sec. 5.3, where I show them explicitly). Contracting 
liww with the tensor q^q'^, one obtains 

/ {0\T*jl^{xmO)\0) = [ {0\T*S^{x)S{0)\0) - ^(0|5i(0)|0) (4.53) 

with = J d"'xe~'^'^'^ , Si the current coupled to the Standard Model Higgs boson, and where 
we have used 

d^J{!,{x) = iS{x) (4.54) 

and 

iCix), S\y)lo^^o 6\x - y)[Si{x) - iS2{x)] (4.55) 

where 52 is the current coupled to the neutral Goldstone boson. This, in turn, can be written 
as 

q'^wwiQ') = -M^n$$(g2) _ ^(o|5i(0)|0) (4.56) 

This equation relates the longitudinal part of the vacuum polarization function of the W 
boson to the self-energy of the corresponding unphysical charged boson. One can see that the 
subtraction term (0 151(0)10) (a tadpole) is needed to cancel a spurious quartic dependence on 
the mass of the fermions. 

Even though the previous derivation was at the one-loop level in the electroweak inter- 
actions, it is valid at any order in the strong interactions as the QCD generators commute 
with the ones of the electroweak group. To derive the self-energy of the charged Goldstone 
boson at 0{a) and 0{aas), we therefore need only the expressions of the electroweak vector 
boson self-energies given in Ref. [59] in the general case and the one of the tadpole diagrams of 
Fig. 4.7 where both the two quarks of the same weak isodoublet are running in the loop. Using 
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(a) 




(b) 




(c) 



Figure 4.7: Tadpole diagrams relating the Higgs boson self-energies to the vector boson self- 
energies at the one-loop (a) and two-loop (b) levels, and mass counterterm (c) contributing at 
two-loop. 



the same notation as in the previous sections, and for a single quark of mass niq (renormalized 
on-shell), we obtain for the tadpole amplitude up to order 



(0|5i(0)|0) 



3m^ 
4x2 



- + 1 



-30 



281n^ 



12 In 




(4.57) 



This equation, added to the one and two-loop expressions for the longitudinal part of the W 
boson self-energies in the general case m„ ^ ^ given in Sec. 4.1.2, leads to the one and 
two-loop expressions of the charged Higgs boson self-energy given in Eqs. (4.41,4.50). This 
is just because for tan/? = 1, the charged Higgs boson couples to fermions exactly like the 
charged Goldstone of the SM, up to a relative minus sign for up-type and down-type quarks. 

In a completely analogous manner, one can derive the Ward identity in the case of the 
neutral Goldstone boson $2 

,2 



g^H^^(g^) = -MjH$,$,(g^) 



9 



-(0|5i(0)|0) 



(4.58) 



4 cos^ 9w 

which allows to check the expressions of the self-energies of the pseudoscalar Higgs boson. For 
tan/3 = 1, the latter has exactly the same couplings as the neutral SM Goldstone, again up to 
a relative minus sign for isospin up and down quarks. 

This result provides a powerful consistency check of both the calculation of the electroweak 
vector boson self-energies described in Sec. 4.1, and the one of the neutral pseudoscalar and 
charged Higgs boson self-energies presented in Sec. 4.2. 
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4.4 Dispersive approach 

In the preceding sections, I have discussed the calculation of vector boson and scalar self- 
energies using dimensional regularization. An alternative method consists in the use of dis- 
persive integrals of the imaginary parts of the VPF's. The imaginary part of the self-energy 
of a vector boson or of a Higgs scalar is related to the partial decay width of the particle 
into quarks. This in turn can be easily calculated at one-loop by integrating in the phase 
space of the fermions the tree level decay amplitude. At two-loop, 0{aas), the calculation of 
the decay amplitude must be pushed up to include virtual gluon corrections as well as real 
gluon production. This procedure is equivalent to the use of Landau- Cutkosky rules [64] on 
the relevant one or two-loop Feynman amplitudes, and implicitly assumes the use of on-shell 
fermion masses.^ Once the imaginary part is known, the real part can be obtained solving a 
dispersive integral 



where the cutoff scale A ^ ^/s, rUa, m,}) regulates the ultraviolet divergence, and a general 
dependence on the masses nia^b is understood. However, n(s) ~ s for large s, and in order 
for the previous integral to be meaningful, the dispersion relation of Eq. (4.59) needs to be 
subtracted. For a conserved current, as in the case of QED, the presence of spurious quadratic 
divergence is removed by a subtraction of the integral at s = 0, in this way restoring the Ward 
identity which requires the photon VPF to be transverse and proportional to s. The situation 
becomes more complicated for a non-conserved current, for which some logarithmic divergences 
are expected to be proportional to the squared fermion masses, and should not be removed in 
the subtraction procedure [56]. 

Subtraction procedures for Eq. (4.59) have been proposed in the past by Chang, Gaemers, 
and Van Neerven [56] and, more recently, by Kniehl and Sirlin [90,91], The dispersion relation 
proposed by the authors of Ref. [56] reads 



and was based on the requirement that the subtraction remove all quadratic divergences 
and that at one-loop level the resulting amplitudes amount to finite renormalizations of the 
result in dimensional regularization, so that the subtraction terms cancel out in physical 
observables, like the shifts of the vector boson masses. The prescription proposed in Ref. [90,91] 
instead manifestly preserves the generic Ward identities constraining the VPF [90]. In the 
following section, as an application of the calculation of section 4.1, 1 generalize the approach 

^The on-shell mass is the parameter that governs the start of the cut on the complex s-plane in S-matrix 
theory, and is naturally associated to the application of Cutkosky rules. 




(4.59) 




(4.60) 
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of Kniehl and Sirlin to the case of arbitrary quark masses, and verify the consistency of their 
prescription up to two-loop level [60]: it is clear that for the prescription to be viable, the 
result of perturbative computations of physical observables in dimensional regularization and 
via dispersive integration must be exactly the same. 

In addition to provide an alternative method for the computation of amplitudes within 
perturbation theory, the dispersive approach can be used to implement additional contribu- 
tions to the amplitudes, whenever perturbation theory is either not applicable or not reliable. 
For instance, in the case of the light-quark contribution to the running of a, where the photon 
VPF has to be evaluated at zero momentum transfer and the effective scale is of the order of 
the masses of the quarks, perturbation theory cannot be used, and the standard procedure 
is to employ experimental data for the cross section a[e'^e~) hadrons and obtain through 
the optical theorem the imaginary part of the photon VPF as a function of s. The photon 
self-energy is then calculated by evaluating numerically the dispersion integral [9,12]. 

Another potentially relevant case where perturbation theory fails to give an acceptable 
description is associated with a heavy quark threshold. Perturbation theory up to 0{aas) 
predicts a discontinuous step-like threshold behavior, which is a rather unrealistic approxima- 
tion. On the other hand, one expects a dense spectrum of resonances just below the threshold, 
and only for a heavy top, above ~ 130GeV, the decay width becomes large enough to smear 
out the resonances, as the lifetime becomes shorter than the revolution period of a tt bound 
state [92]^. It is also well-known that the Coulombic effects near threshold cannot be accom- 
modated when the perturbation series in is truncated at a finite order, as they correspond 
to soft multi-gluon exchange.^ 

An acceptable description of a[e~^e~ tt) therefore requires contributions to the absorp- 
tive part of the photon and Z self-energies besides the usual perturbative ones. As the t — h 
doublet gives the leading electroweak corrections to most precision observables, the possibility 
arises that additional effects connected to the ti threshold be important in precision physics. 

The question of whether these additional "non-perturbative" contributions to the imagi- 
nary parts of VPF's affect the real parts, and consequently the radiative corrections to elec- 
troweak precision observables, has been the subject of considerable attention in the last three 
years (see Refs. [90,91,36,51,92] and references therein). Several methods have been proposed 
to take into account the ti threshold contributions, based on the use of phenomenological 
potentials like Richardson's in the context of a non-relativistic approximation: a thorough 
discussion can be found in Ref. [91,36]. The resulting effect on the electroweak parameters 
seems to be generally sizable, even though with large theoretical uncertainties: in Ref. [36] 
it has been estimated that the threshold effects would enhance the usual QCD perturbative 
corrections to the leading top contribution of Eq.(4.33) by 25-40%. 

^Compared to it resonances, the contribution of tfe bound states is small, as threshold effects are proportional 
to the squared reduced mass of the quarks, and can be safely neglected [93] 

^The resummation of these long-distance effects amounts to multiplying the bare spectral function around 
threshold by the factor |-7r-^ (1 — e"^"""''^"')"^ , for which the expansion parameter is as/v rather than as, 
with V = \/^ — 4m( / s the relative velocity of the quark pair in the center of mass system. 
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The suggestion has been made, however, that when one works at a specified order in per- 
turbation theory, the eff"ect of the potentially large threshold contributions on the real part of 
the VPF should cancel out in the dispersive integration over the real axis, if all perturbative 
and non-perturbative contributions are taken into account [94,96], as a result of the analyt- 
icity of the S-matrix. This is in contrast to what happens in many other cases where the 
absorptive part of the VPF is not integrated over the cut in the complex plane (for instance, 
top production, e~^e~ t~^t~ [95]). In other words, although individual hadronic bound 
states are governed by long-distance dynamics, and are therefore inherently non-perturbative, 
the dispersion integrals over these states are determined uniquely by the short-distance dy- 
namics, and can be safely calculated in the framework of perturbation theory [97]. This point 
of view is supported by arguments by Takeuchi and collaborators [99], and Sirlin [51], that I 
briefiy summarize in section 4.4.2. A significant progress in clarifying the problem of QCD 
corrections to electroweak observables has however come from a recent three-loop calculation 
of the QCD corrections to Ap [103]. I will discuss some of the consequences of this new result 
in Sec. 4.5. 

As the subject of the role of non-perturbative QCD corrections to precision observables is 
not settled to date, I will not dwell further on it. What is of relevance here is that dispersion 
relations might possibly be of use to implement non-perturbative eff"ects, and that it is impor- 
tant to know the precise relationship between them and dimensional regularization. Indeed, 
as will be illustrated in Sec. 4.4.2, the implementation of non-perturbative eff"ects depends 
sensitively on the subtraction prescription that is used. In order to test the consistency of 
the procedure, it is then crucial to verify whether the results for physical observables given by 
the dispersive method in perturbation theory coincide with the ones obtained in dimensional 
regularization. Finally, once the prescription has been chosen, it may be useful to implement 
threshold eff"ects in the MS scheme, and the correspondence between the two methods must 
be elucidated. These are the points that I will investigate in the next section, using the results 
obtained in the first part of this chapter. 

4.4.1 Connection with dimensional regularization 

First of all, let us recall the basic elements of the procedure proposed by Kniehl and Sirlin 
(for a more detailed discussion see Refs. [90,91,36]). Using Eq.(4.1), we can write the Ward 
identity 



where I dropped the indices i, j of the vacuum polarization tensor for simplicity. Contract- 
ing the previous equation with q'^, one obtains A(s) = 11^(5); then introducing the linear 
combination $(«) = (11^(5) — I[t[s))/s, the cofactor of q^q'^ in Eq.(4.1), we can rewrite the 
transverse part in terms of A(s) and $(«), 




(4.61) 



5$(5) + ni(5) 



5$(5) + A(5), 



(4.62) 
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The advantage of expressing 11^ in this way is that and A(s) are only logarithmically 

divergent for s oo, unlike 11^(5) which is quadratically divergent as discussed above. For 
this is obvious for dimensional reasons (two powers of the external momentum have been 
extracted), while for A(s) this is because the soft breaking of the current Jfj^ by mass terms 
implies that its divergence involves at most operators of canonical dimension three, and one 
power of the momentum has been extracted. It is therefore possible to write unsubtracted 
dispersion relations for and 11^(5) = A(s). The use of Eq.(4.62) will then produce a 

dispersion relation for the transverse part of the VPF that automatically enforces the relevant 
Ward identities. Re-expressing Tmlli in terms of Tm$ and ImllT, we obtain 



ds' 

I (mi+m2)^ 



— + Tm$(s ) 



1 

n^{s) = - / 

TT J{mi- 



where the quadratic divergence is now removed by the contribution of Tm$(s'), and the 
logarithmic divergence is regulated by the cutoff A. The integral over the large circle, where 
the s in the denominator has been neglected compared to s', comes from the fact that 11^(5) 
and $(5) do not vanish asymptotically but behave like constants modulo logarithms, so that 
they must be considered in the application of Cauchy theorem. Note that for conserved 
currents dfj,J^{x) = 0, and 11^(5) vanishes reducing then the Ward identity to the QED form 

We can now analyze the relation between the two calculational approaches at the one and 
two-loop level, and use the one and two-loop results in dimensional regularization of section 
4.1 to write the real parts of Il^'^(s) as 

~ - TtlP' - 

K'^{s) = sX + {ma T mi,fY + {ml - ml) In^ Z + F'''^{s), (4.64) 

where the divergent constants X, Y and Z involve only the poles in e and the logarithms 
of the scale /x of the expressions of section 4.1, the function F^'^ behaves like a constant as 
s ^ 00, and the — ( + ) sign refers to vector (axial- vector) current. The one and two-loop 
components of the functions F^'^ and of the divergent constants can separated by writing 

X = Xi + ^X2, f = fi + —Y2, Z = Zr + ^Z2. 

Their expressions at the one-loop level are given by (see Eq.(4.6))^ 

- 1 mgmi, - 3 3 mgmi, 

Xi = --ln-^, F, = -_ + -ln— ^, (4.66) 



^The definition of F^2^{s) used here diff'ers from the one used in Ref. [91], as we normalize the divergent 
constants to the dimensional regularization result Eq. (4.64) keeping only the ultraviolet divergences and the 
related logarithms; remember also that ln47r and 7 constants have been absorbed in the normalization of the 
VPF. 
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with Zi = 0, while at the two-loop level, as can be easily read off from the asymptotic 
expression for s ^ oo reported in App. B.l, they are given by 

- 1 ma,mb 
Xo = m — 

- 3 1 /II ma,mb\ 11 ma,mb 2 "^<i"^6 
12 = :r-:7 H — I ^ 3 in — ) ^ in h 3 in 



2e2 e V 4 fi^ J 2 fi^ fi^ 

^2 = -— + 31n —. (4.67) 

If one integrates the imaginary parts of the vector bosons VPF's at the one and two-loop 
level (given for example in App. B.2), the result of the dispersive integral Eq. (4.63) for 
arbitrary quark masses and for vector and axial-vector currents can be expressed as 

ml 



ml 



n^'^(s) = sX + {m^T rribfY + {ml - m^) In ^ Z + F^'^(s) 

ds' 

2xi /|s'|= 



where the upper sign again refers to vector current. Here the divergent constants X,Y,Z 
involve logarithms of the cutoff scale A, and their one and two-loop parts can be separated as 
in Eq.(4.65); their explicit form is given at one-loop by 

3 A2 

Xi=ln -, yi = --ln + 3, (4.69) 

mamb 3 2 mamb 

with Zi = 0. The finite function F^'^[s) must be the same as the one in Eq. (4.64). Since 
the result of Eq. (4.63) should be consistent with the one obtained using dimensional regular- 
ization, the integral over the large circle must be absorbed in a redefinition of the divergent 
constants X, Y, Z. In particular, if 11^(5) and $(5) are evaluated in dimensional regularization, 
the eff"ect of the second integral in Eq. (4.63) should be simply to replace Xi^2? 5^i,2? and Zi^2 
by Xi^2? 5^,2? and ^1,2- This is indeed true at one-loop level, and one has 



iT^l J\s'V 



2wi /|s'|=A2 

A_ I d^'Hl:^ = (^^ ^ rnbf{Y^ - Fi). (4.70) 

This point is crucial to ensure that the cancellation of divergences occurring in physical ob- 
servables like the p parameter and Ar yields a total independence of these observables from 
the calculational approach used to obtain the vacuum polarization functions. 

At the two-loop level the dispersive integral of Eq.(4.63) has been calculated in the cases 
nia = mb and mb = [57], for which it has therefore been verified that the functions F^'^ 
in Eq.(4.64,4.68) coincide up to 0{aas) once the integrals over the large circle are evaluated. 
In the case of arbitrary masses the integral Eq.(4.63) has not been calculated yet, but the 
divergent part (involving logarithms of A) can be easily obtained and is reported in Ref. [93,57]. 
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It is therefore possible to use our asymptotic expressions for the 0{aas) VPF's of App. B.l to 
evaluate the integrals of the left-hand side of Eq.(4.68), obtain the two-loop constants X2,Y2 
and Z2, and compare at least the divergent part. Here again, since the result of Eq. (4.68) 
should be consistent with the one obtained in dimensional regularization, the integral over 
the large circle must be absorbable in a redefinition of the divergent constants X2,Y2 and 
Z2. Thus, when we use dimensionally regularized expressions for $(5) and 11^(5), the integral 
involving $(5) in Eq. (4.68) will replace X2 by X2, and the one involving 11^(5) will replace 
Y2 and Z2 by respectively Y2 and Z2. One has therefore at two-loop 



1 / .j^l:^) 

s' 

,nr'^(.') 



^ ds' ^ = X2- X2, 

2^1 J\si\=A2 S' 



L I = (m„ T m,f{Y2 - Y2) 

2X2 J\s'\=A^ S' 



+ {ml-ml)ln^ (^2-^2). (4.71) 

5 



Assuming that the procedure is consistent, we solve for X2,Y2, and Z2, and obtain' 

^2=ln^+4C(3)-f| 
niamb 12 

Y2 = hn^^-hnJ^-hn^^+'-- 6C(3) - '-.^ 
2 mam,}, 2 mam}, 2 m}, 8 4 

3 A^ 

Z2 = --lu + 5 (4.72) 

2 mam}, 

The divergent part of these constants agree with Ref. [57,93]. Inserting Eq.(4.71) into Eq.(4.63) 
finally leads to 



0^(5) = - / ds' + Tm^{s') 



s — le 



(4.73) 



ml 



s{X -X) + {ma T m},y{Y - Y) + {ml - mf) \n^{Z - Z) 

m^ 



Using Eq.(4.73), it is now possible in principle to implement in the general case of arbitrary 
quark masses non-perturbative eff"ects [91,36] also in the MS scheme for which a dimensional 
regularization calculation is required. This result [60] generalizes part of Ref. [91]. 

4.4.2 An application 

Let us now briefiy discuss the application of the dispersive method to physical observables, 
taking again the example of the fermionic contribution to the p parameter. Using the fact 
that at zero-momentum transfer the transverse and longitudinal components of the W and 



^Note that the sign of the C(3) term in Y2 in Ref. [60] is incorrect due to a misprint, and that the expressions 
for the special cases considered in Ref. [91] were obtained using Ref. [54], and are therefore affected by the 
omission of the tt^ term in the mass renormalization mentioned in section 4.1.2. 
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Z self-energies are equal, and recalling that the vector component vanishes in this limit for 
equal mass quarks, we have written Ap in terms of the the value at s = of a function 
fh^T^s, TTia, m,})) defined in section 4.1.3. Since Ap is an observable and finite physical quantity, 
fL^T^^iiT^ain^b) can be evaluated by means of an unsubtracted dispersion relation 

fL,T{^^"^a,mb) = - ds j : (4.74) 

TT J(mi+m2y s' - S - 16 



The condition for Eq.(4.74) to be valid is that 

ds' 

A2^oo y|s'|=A2 



/ ds 

lim <p — fL,T{s',ma,mb) = (4-75) 

^U'I=A2 S' 



as a consequence of Cauchy theorem. 

Using Eq.(B5), it is easy to verify at the two-loop level that /i,,t(-5, rrZa, m;,) is free of 
ultraviolet divergences and tends to zero as s ^ oo and therefore satisfies the condition of 
Eq.(4.74) in the case of arbitrary masses. Note that this happens for both transverse and 
longitudinal components of the VPF. This means that, as long as perturbative eff"ects (up to 
O(aas)) are concerned, a dispersion relation based on the longitudinal components as the one 
proposed in Ref. [90,91] ^ gives the same result as a dispersion relation based on the transverse 
components, as the one proposed in Ref. [56], despite the very diff"erent asymptotic behavior 
of the individual components. In fact, the two subtraction prescriptions of Refs. [56, 90] 
coincide at one-loop for arbitrary quark masses, but they diff"er by finite renormalizations to 
leading order in QCD (with the obvious exception of conserved vector current). The extra 
pieces cancel however in physical observables like Ap and Ar^, giving the same results as 
in dimensional regularization. This is by no means a general feature: it has been shown in 
Ref. [98] that in the case of the leading order QCD corrections to the decay amplitude of the 
SM Higgs boson into leptons, only the prescription of Ref. [90] reproduces correctly the result 
obtained in dimensional regularization. 

The observation, Eq.(4.75), that the two dispersive procedures are equivalent up to 0{aas) 
motivated the work of Ref. [99], where the Operator Product Expansion (OPE) was used 
to show that non-perturbative QCD corrections to Ap can be described equivalently using 
dispersion relations that involve longitudinal or transverse components of the VPF's. The 
basic idea of Takeuchi and collaborators is that the non-perturbative asymptotic behavior of 
the VPF's nT^i(s) as |s| oo can be extracted from their OPE. The latter can be found in 
the appendix of Ref. [100], together with the relevant Wilson coefficients Ci{q), which can also 
be extracted from the high-momentum expansion given here in App. B, 

nr'^(g') = (m„(g) T M<l)YC^,i{q) + (m„(/i) T mb{f^)YC^,2{q) + • • • , (4.76) 
$-(.^) = + C.,(,) (^M^l^^^ + ^^^^(^) (^.(.)Trn.(,))^ ^ 



q^ 



^Note that for s = the integrand in the first line of Eq.(4.63) reduces to ni(s')/s'. 
'^I have verified this explicitly in both cases for arbitrary masses. 
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where the ellipses stand for terms suppressed by inverse powers of g^. The long distance 
dynamics is embodied in the running masses rn(/x) in terms of which the OPE is built (products 
of masses are the only dimension two operators available), and in the vacuum expectation 
values of operators of dimension higher than two which are suppressed by inverse powers of 
g^. The conclusion is that, since the Wilson coefficients depend only logarithmically on q, 
Eq.(4.75) is valid at any order in QCD and even including non-perturbative QCD corrections. 

The immediate consequence of this result is that the incorporation of non-perturbative 
threshold effects through either of the two dispersion relations should give the same result.* 
The calculations of Refs. [57, 99] in the context of a non-relativistic approximation instead 
indicate that the two approaches give opposite results! 

Eq. (4.75) and the argument by Takeuchi et al. can be recast in a different form, which is 
particularly simple. Sirlin has observed [51] that the difference of the two forms of Eq.(4.75) 
can be expressed as the integral over the large circle of a combination of functions $^'^(s) 
(the denominator cancels). Applying Cauchy theorem again, the result can be alternatively 
phrased as the dispersive integral of a combination of spectral functions Tm$(s), a sum rule 
similar to the well-known Weinberg's [102] 



With this sum rule in mind, Sirlin argues that, in the case of the top-bottom doublet, poten- 
tially large threshold effects can arise only from Tm$^'^(m4, m^), because the other channels 
are suppressed by reduced mass considerations. Consequently possible large threshold contri- 
butions should cancel in the integral /°° dslm[^^[s, rrit, rrit) + $^(-s, rrit, rrit)]. The possibility 
that they cancel against each other seems inconsistent with the non-relativistic approxima- 
tion developed in Ref. [91], where it is found that Tm$^(s, m^, m^) ~ Tm$^(s, m^, m^) in 
the threshold region. The only reasonable conclusion is that large threshold contributions 
to Tm$^'^ cancel against other contributions when the integrals over all s are considered. 
These integrals can be related to integrals in the asymptotic region, and can be expanded per- 
turbatively, so that the cancellation must occur order by order in perturbation theory. This 
argument, without proving rigorously the cancellation of the threshold contributions, seems 
to point in the same direction as Ref. [94]. 

4.5 The choice of the scale of 

Before concluding this chapter on QCD corrections with a overview of the incorporation 
of QCD corrections in the calculation of precision electroweak observables, let us turn our 
attention to a central issue for any QCD calculation: the choice of the scale at which the 

'in Ref. [101] the same authors of Ref. [99] apply the OPE method to the H ^ II case, and confirm to all 
orders in perturbation theory the result [98] that only the dispersion relation based on the Ward identity is 
consistent with dimensional regularization in that case. 




i=a,b j=V,A 



(4.77) 
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running coupling constant must be evaluated. In the preceding sections I have purposely 
avoided this problem, but as soon as we consider a practical application of the calculations 
illustrated above, we are bound to specify the scale, and the question cannot be further 
postponed. 

In general, QCD predictions for a measurable quantity A have the form 



The coefficients Ci{Q) depend both on the exact definition of (i.e. the scheme) and the 
choice of the scale Q. If a calculation could be performed at all orders in as{Q), the choice of 
scheme and scale would be immaterial: A must be the same in any scheme, and independent 
of the scale chosen. If the perturbative series is truncated at a finite order, however, the 
prediction for A is both scale and scheme dependent. For instance, a redefinition of coupling 
or a change of scale can, in principle, change the next-to-leading order coefficient to any 
desired value. Therefore, in order to set the scale appropriate to a particular process so that 
Ci, C2, ... be meaningful, we need some sort of information on the higher order contributions 
of the perturbative expansion. There is no rigorous way to gain this information. It is possible, 
however, to make an educated guess. 

Let us concentrate again on the contribution to Ap of the (t-b) doublet. As I have men- 
tioned above, there is a very recent result of a three loop calculation [103] in the limit of very 
heavy top that modifies Eq.(4.33) into 



where I have not specified the scale, and I have used the on-shell top mass m^. As the only 
mass scale in the process is rrit, it seems reasonable to set = as{mt). It turns out that this 
is not the case, and the fact that we know only the first two coefficients of the series implies 
that the error that we make by choosing the scale in Eq.(4.79) at rrit is possibly relevant. 

For this particular case, we have two diff"erent methods that give us very similar answers. 
The first method consists in considering some particular type of potentially large eff"ects, like 
the gluon vacuum polarizations, and in trying to take them into account at higher orders. 
A heuristic procedure has been proposed for example by Brodsky, Mackenzie, and Lepage 
(BLM) [105]. An alternative procedure is based on an analysis of the scale dependence: the 
best choice for the scale is the one for which the scale dependence is minimal (Principle of 
Minimal Sensitivity, PMS) [106]. One can also use the Fastest Apparent Convergence (FAC) 
method, which consists in absorbing the whole QCD correction in a scale redefinition of the 
first term of the series. BLM, PMS and FAC applied to Eq.(4.79) give us somewhat similar 
answers [51], suggesting a smaller scale than rrit, and therefore a larger QCD correction. They 
also suggest large and increasing higher order coefficients for the series Eq.(4.79). 

The second method consists in the use of an eff"ective lagrangian approach, where the top 
is assumed to be very heavy, and integrated out. This approach is illustrated, for instance. 




(4.78) 




(4.79) 
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in Ref. [104], and the result is that the scale appropriate for Eq.(4.79) is rrit, provided that 
Eq.(4.79) is expressed in terms of the MS top mass calculated at m^. As shown in Ref. [51], 
the two approaches numerically give very similar answers. A very detailed discussion of all the 
points raised here can be found in [51], where it is also shown that a very good convergence 
pattern is obtained when the ratio between the pole and the MS mass 77x4(77x4) is first optimized 
according to one of the methods (BLM or PMS) mentioned above, and then used in Ap. The 
result of Ref. [51] can be cast in a very convenient form for numerical calculations: 



Ap 



top 



+ 3)^). (4.80) 



where ^ = 0.260^0 056- Numerically, for rrit = ISOGeV, this implies a QCD correction of 
11.95 ± 0.51%. The relative error on the QCD correction is 4.3%, and has a negligible impact 
on the electroweak observables. For instance, it induces a shift in the determination of Myy of 
^ 3MeV, and of 2 X 10"^ in sin^ ^^^(M^). 

It is clear that large higher order corrections are linked in some way to the sensitivity of 
the process under consideration to infrared gluons. For example, it is easy to realize, and can 
be explicitly checked [94], that the diagrams of Fig. 4.1b are not sensitive to low momentum 
gluons, while the one of Fig. 4.2a, having the resonance of the top propagator exactly at the 
momentum transfer, is certainly very sensitive to the gluon infrared region. This is why, 
through the on-shell counterterm (diagram Fig. 4.2a evaluated at = m^), the pole mass 
introduces in the perturbative series of Ap large terms of 0{{n — l)!/3o~^a"), where /3o is the 
first coefficient of the (3 function, connected to the appearance of singularities in the Borel 
plane, the renormalons. 

For what concerns the electroweak observables of interest, the VPF's have to be eval- 
uated at equal to M| or M^, and the arguments used for Ap cannot be applied. An 
effective field approach [104] again suggests to use 0^(77x4) for the leading quadratic contribu- 
tion expressed in terms of rht(rnt), while the logarithm In-r^ has to be shifted: In 

ln^-4//3oln^^^. 

4.6 QCD corrections to basic electroweak observables 

In this section I consider the QCD corrections entering the determination of M^i^ and sin^ dw{Mz). 
I remind that QCD corrections for the form factor fc(M|) that connects sin^ ^^^^ and sin^ 6w{Mz) 
have been considered already in Sec. 2.4. Unlike the case of processes involving external 
hadrons, like e"*" e~ hadrons, which receive pure QCD radiative corrections, the tree level 
predictions for and sin"^ 9 w{Mz) are affected only at 0{aas), by gluonic corrections to 
quark loop insertions on vector boson propagators. Despite the large value of the QCD cou- 
pling constant, we can therefore expect this kind of effects to become relevant only when the 
one-loop correction is already sizable. This is the case of the large top contribution to the 
radiative corrections Ar and Ar. 
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The incorporation of perturbative QCD effects in the basic electroweak corrections Ar and 
Ar has been considered in detail in Ref. [36]. Here I give just a summary of their results, in 
order to make the implementation of 0{aas) effects possible, using the results of this chapter 
and of App. B. Normalizing the VPF's as I have done throughout this chapter, and extracting 
from Azzi^q^) and Ayy^yy^i^q^) a factor ^ at one-loop level and at two-loop, we can write 
the contributions of a quark isodoublet (u, d) to the Z and W self-energies as 



i=u,d 



(4.81a) 



and 



(4.81b) 



where Qi and Csi = ±1 are the electric charge and the isospin eigenvalue of the quark flavor, 
respectively. Including the decoupling of the top quark according to the Marciano-Rosner 
convention [38], and neglecting terms that contribute in O(d^) but not 0{aas), we find [36] 



Ar' 



doubl 



Re 



'Azz{Ml) Ay^y^iOy 



Ml 



M2 



i=u,d 



+D 



MS 



(4.82) 



where the subscript MS indicates that the poles in e have been removed and that /x has been 
set to Mz- -D is a contribution deriving from the use of the decoupling as described in section 
2.3, and can be written as 



^ a 2P 
^ = a-2-2 



3* 



^ In^ 



15 
4 X 



(4.83) 



Similarly, we obtain for Ar in the on-shell scheme 



Ar' 



doubl 



W i=v,A 



^ n'(0, m„,md) 



=u,d 

16sW^^ E [(l-452C3iQi)2n^(M|,mi,mi) + n^(M|,mi,mi)] 

^ i=u,d 



Ss^Ml 



Re [n^(M^, m„, ma) + n^(M^, m„, ma)] 



(4.84) 



The expression for Ar in the MS scheme differs from Eq.(4.84) by subleading terms of order 
O(a^), that I do not consider here. We see that for rriu > Mz,Mw, /^j,<ioubl reduces to 
c^/s^A/9. One can use for this leading term the result of the three loop calculation of Ref. [103], 
or Eq.(4.80). From the argument of the preceding section, it seems convenient to use MS 
masses for the quarks. To this end, Eqs. (4.16) and (4.17) have to be employed. 
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Now all the ingredients for a detailed evaluation of the hadronic contribution to Ar and Ar 
have been given. Using the results of this chapter and of App. B, together with the estimate of 
higher order corrections of Ref. [51], a very accurate determination of the QCD effects in the 
prediction of and sin^ 6w(Mz) should be possible, including the light quark mass effects. 
Of course, different options are possible (the use of Eq.(4.80) or (4.79), of on-shell quark masses 
or MS, etc.) and they all contribute to define a theoretical error for the QCD corrections. 
This can be expected to be beyond the present and future experimental resolution. 



Chapter 5 

Higher order electroweak 
corrections 

Oh, il difficile non sta nel drizzare I'uovo di Colombo, ma nel covarlo. 

E. Flaiano 

Potentially large two-loop contributions originated by electroweak interactions are the subject 
of this chapter. After a brief introduction, I will present the main features of a calculation of the 
p parameter in the SM in the limit of heavy top mass [107-109]. Leading and next-to-leading 
terms in rrit are kept up to 0{G'^mfM^). Although the result cannot be applied directly to the 
case of Mw and sin^ ^e//' extrapolation indicates that next-to-leading 0{G'^mfMl) terms 
can contribute significantly to the theoretical prediction of these observables, and provides an 
estimate of the current theoretical error in the determination of high-energy observables. 

5.1 Higher order effects of electroweak origin 

In Sec. 2.5 I have anticipated that large masses, like rrit and M^, potentially induce large 
radiative corrections to precision observables. This can be regarded as a consequence of the 
violation of the decoupling theorem occurring in the SM. The investigation of the actual impact 
of the large masses on the high-energy phenomenology at the two-loop level is therefore of 
great importance. Part of this program has already been completed with the study of the 
leading dependence on rrit and of all the precision observables at two-loop level. 

5.1.1 Large Higgs mass corrections 

An elegant eff"ective Lagrangian approach has been used in Ref. [110] to investigate the leading 
two-loop eff"ect of a large Higgs mass on the various precision observables in electroweak 
physics. Barbieri et al. have shown that the physical corrections can all be reduced to 
the corresponding eff"ects in the two parameters ei and €3 [112] (or in a diff"erent language T 
and S [113]), that are often used to parametrize the eff"ect of new physics on the precision 
observables. Self-energy diagrams that contribute to 0{g'^M^) are present, but cancel in 
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all relevant physical amplitudes. The contribution to the isospin breaking parameter ei is 
the same as the one to the traditional Ap, and the actual calculation yields Ap ~ 5.03 X 
10~^ (Mjj/TeV)^, when the one-loop corrections are expressed in terms of Gfj,, a, Mz, and 
the on-shell mass M^- This result confirms the pioneer work of Ref. [111]. For a Higgs mass 
below ~ 2TeV, the correction is extremely small when compared to the current accuracy. 
Since Mjj^lTeV is unlikely to be consistent with a perturbative treatment of the SM in 
general, this means that the two-loop corrections to the p parameter and to the theoretical 
determination of and sin^ dw{Mz) are practically irrelevant. As an illustration, keeping 
such a contribution in Ar and Ar would shift by 6Mw ~ 3MeV, and sin^ dw{Mz) by a 
mere 1.7x10"^, for ~ ITeV. A similar pattern is found for the parameter €3, which enters 
the decay width of the Z boson, Aes ~ 3.16 X 10"^ (Mjj/TeV)^. Consequently, the effect on 
the Z width, Tz, is also tiny, ATz/Tz ~ 0.53 X 10"^ {MH/TeYf. 

5.1.2 Large top mass corrections 

The leading two-loop corrections in the limit of a heavy quark doublet have been first calcu- 
lated some time ago [114] in the approximation of a light Higgs mass, <C rrit, and found 
quite small. As the experimental bound on became higher, the approximation used in 
Ref. [114] turned out to be inadequate. This motivated a calculation [115] of the leading G'^mf 
corrections in the case of a heavy Higgs boson whose mass can be comparable to rrit (with 
m,}) = 0). Rather unexpectedly the result, a function of the ratio Mu/rrit, was found to give 
significant contributions to the precision observables, of the order of the future experimental 
accuracy, for values of comparable to the top mass. 

To study this kind of effects it has proven useful, both from a conceptual and a practical 
point of view, to think of the electroweak sector of the SM as a weak perturbation around an 
underlying pure Yukawa theory, obtained from the SM by turning off the gauge interactions 
[118]. Such a theory describes the ideal case of top and Higgs masses much larger than the 
gauge vector boson and light fermion masses, with gauge interactions becoming negligible 
compared to those associated with the top-scalar sector. In this scenario all relevant physical 
observables of present interest are evaluated at a momentum scale much lower than the top 
or Higgs masses and that can be set to zero, at least in first approximation. This makes 
it possible to relate physical quantities of interest, e.g. the ratio between Myy and Mz, 
to Green functions of the Yukawa theory. Indeed, in the spontaneously broken phase of 
the SM, a set of Ward identities relates vector boson Green functions at vanishing external 
momenta to appropriate Green functions of the corresponding Goldstone bosons. For instance 
the relevant self-energies entering into the radiative corrections to the p parameter can be 
expressed, neglecting subleading contributions in the heavy top, as (see Sec. 5.3) 



(5.1a) 
(5.1b) 
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where the A(0)'s are the transverse part of the vector boson self-energies at zero momentum 
transfer and the II's are the two-point functions associated with the unphysical counterparts. 
Using this correspondence one can, in practice, derive the leading contribution to the radiative 
corrections working entirely in the framework of the Yukawa model without reference to the 
original theory, this framework allowing a considerable simplification in the evaluation of 
physical eff"ects at higher orders. Along these lines, the leading 0{G'^mf) two-loop terms have 
been evaluated for both the p parameter and the Z ^ bb partial width, for an arbitrary ratio 
of the top and Higgs masses [115-117,107]. These corrections represent eff"ects proportional 
to gf, that depend only upon the two parameters of the Yukawa Lagrangian gt and A, the 
top Yukawa coupling and the quartic self-interaction of the scalars, respectively. Ap is the 
result of interest here, because it is proportional to the leading top contributions to Ar and 
Ar, the relevant radiative corrections for the determination of and sin^ dw{Mz). It can 
be written as 



Ap' 



(5.2a) 



with 



where R[x) is a function that reaches a minimum of —11.8 for x ~ 5.7, and Xt = 3 ^^'^^ . The 
eff"ect of such a contribution can be estimated using Eqs. (2.15) and (2.21). When the function 
R reaches its minimum, and for rut = ISOGeV, SM^ ^ 23MeV and ^ sin^ 0[j^^ ^ 1.4 X 10"^. 

While the above picture is theoretically clean and very useful for practical computations, 
up to now its limitations have not been fully clarified. Concerning the Higgs, it is still possible 
that this particle is much heavier than all the remaining spectrum. Instead, in the case of 
the top, we have seen that there are strong indications [1,2] that its mass is only moderately 
larger than the vector boson masses, rrit ~ 2Mz. This raises the question of the corrections to 
the asymptotic regime here referred to as the Yukawa limit. Compared to the leading order 
results, these corrections are of order (M^/m^)^, and could be sizable, as can be verified in 
explicit one-loop examples. Consider for instance the function f[mf/M^) appearing in the 
evaluation of the one-loop box diagram for B° — B° oscillation: 

1 9 3 3 z^lnz 

^^^^ " 4 + 4(1 - z) " 2(1 - z)2 + 2 (z - 1)3- ^^-^^ 

The asymptotic value /(oo) = 1/4 is rather far from the realistic value /(5) ~ 0.54, and in 
this case the subleading contributions amount to a 100% correction. 

Coming to two-loop eff"ects, it is important to establish whether the present estimates 
of the purely electroweak top contribution, based on a calculation in the Yukawa limit, are 
realistic or can be infiuenced by subleading eff"ects that cannot be neglected in view of the 
future experimental precision. The first step in this direction seems to be an investigation of 
the two-loop next-to-leading top corrections, namely contributions that scale as g^g^ where g 
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is the gauge coupling. To keep the computation as simple as possible I will focus on neutrino 
scattering on a leptonic target. For this process I present a calculation [108, 109] of the 
electroweak corrections of 0(G^mjM|) to the p parameter, defined as the ratio of neutral-to- 
charged current amplitudes, at zero momentum transfer. There are two important advantages 
in considering p instead of other observables: first, the natural relation pQ = 1 (cf. Eq.(2.4)) 
of the bare Lagrangian ensures that only one-loop renormalization is needed^ ; secondly, the 
irreducible two-point functions to be considered are evaluated at zero momentum transfer, 
which greatly simplifies the calculation. The computation can be regarded as an attempt 
to check the validity of the expansion in powers of the gauge coupling constant around the 
Yukawa limit, until the full two-loop results will be available. At the same time we will be 
able to provide an estimate of the error associated with the two-loop electroweak eff"ects. For 
the processes under examination, we find large subleading corrections, of the same sign and 
of about the same magnitude as the leading one. 



5.2 The example of A/>: a complete calculation at (9(G^m^M^) 

In this section I outline the computation of the electroweak corrections of 0(G^mjM|) to 
the p parameter. Specifically, I identify p with the cofactor of the Jz Jz interaction in neutral 
current amplitudes, expressed in units of Gfj,. As well known, radiative eff"ects lead also to 
a modification of the mixing angle, described by a form factor usually called n^q"^). These 
eff"ects will not be discussed here. 

We begin by writing the relation between the /x-decay constant and the charged current 
amplitude expressed in terms of bare quantities. At the two-loop level, neglecting terms that 
do not contribute at 0(G^mjM|), we have 

where go and Myy^ are the bare SU(2)l coupling and W mass, respectively, Ay^rw = ^ww(O) 
is the transverse part of the W self-energy at zero momentum transfer, and the quantities Vw 
and Byy represent the relevant vertex and box corrections. At the bare level, using the fact 
that M^^Cq = M^^, where cq = cos^^i^^, the p parameter can be written as: 

Azz AzzVz 1 



+ + MwoBw + -J^^ 



(5.5) 



where the self-energies are evaluated at = 0, and Yz and Bz are the corresponding ver- 
tex and box contribution in neutral current amplitude. To the order we are interested in. 



natural relation is a constraint between coupling constants and masses that results from the symmetry 
structure of the Lagrangian before spontaneous symmetry breaking. The counterterms needed to renormalize 
the parameters of the broken theory possess the symmetry of the unbroken Lagrangian, and, consequently, not 
all of them are independent. As a result, the radiative corrections to a natural relation are finite at one-loop. 
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remembering that at one-loop only self-energies give rise to terms, Eq.(5.5) reduces to: 

\ Wo Zo / 



We proceed by separating the self-energies into one-loop and two-loop contributions: 

A _ ^(1) , aW. a - aM _l aW 



(5.7) 



with the understanding that the one-loop term is still expressed in terms of bare parameters. 
The one-loop part can be further decomposed into pure bosonic (6) and fermionic (/) terms: 

■^zz — ^ZZ \ ^ZZ 1 ^ww — ^WW \ ^WWl V"""/ 

and the one-loop fermionic contribution to the p parameter, assuming a vanishing bottom 
mass, can be expressed as follows: 

( Af \^^^ 
s ^ 7T^ ">.' « r(€) (^] ' , (5.9b) 



where e is related to the dimension d of the space-time by e = (4 — d)/2 and /x is the 't-Hooft 
mass scale. 

As we want to express the final result in terms of the physical Z mass, we perform the 
shift Mj^ = M| - Re Azz{Ml) ~ M| - Azz- Using the decompositions given in Eqs. (5.7) 
and (5.8), recalling that one-loop vertices and boxes do not involve the top quark, and keeping 
only terms up to 0(G^mjM|), we obtain after simple algebra: 

A, 



( _ \ (A A \ (2) 



'■ZZ 

Ml ) ' \ Ml Ml 



{Vz - Vy,) + MlcliBz - B^) - Xd{Vy, + 2 MlB^) 

(4f -w) +«^-^-) 



(5.10) 



where Xd is the quantity introduced in Eq. (5.9), now expressed in terms of renormalized 
parameters. 

We observe that the Eq.(5.10) further simplifies if we express the one-loop fermionic con- 
tribution in terms of the Fermi constant G^j,. Indeed, as can be seen from Eq. (5.4), the first 
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line of Eq. (5.10) reproduces the effective coupling in the charged sector 



Ml 



+ MlB 



w-'-'w 



_jL_ 



Ml 



Vy, + MlB^^ f{ml,e) 



|/K,e) 



(5.11) 



Up to now, apart from the use of the physical Z mass, we have not specified any partic- 
ular renormalization condition. In order to simplify the structure of the counterterms, we 
have found it convenient to perform the calculation using the MS parameter sin^^^(M2) 
(henceforth abbreviated as i^). Indeed, while in the on-shell (OS) scheme the counterterm 
associated to the quantity = 1 — Ml /M| contains terms proportional to and gives rise 
to 0{GlmfMl) contribut ions to p, the counterterm related to does not exhibit any 
dependence and this greatly simplifies our task (cf. Sec. 2.4)^. Therefore, to the order we are 
interested in, we can directly replace Cg with in Eq.(5.10) (c^ = 1 — i^). It will always 
be possible to recover the result in the pure OS scheme, by appropriately shifting P in the 
one-loop expression for p. 

We now notice that the one-loop contribution is still written in terms of bare quantities. 
To put p in its final form, we split it into the usual 0{a) result, Ap^^\ plus the counterterm 
part, Apc namely 



(1) 



Ml 



{v,-v^)'^^HmIc\b,-b^)'^^ 



with 



Apf('^ 
NcXt = 



Ap'('\ 



Ap 



6(1) 



a 



4X5^ 



8x2^2' 

— - in c 

4i2 4 



2 cz 



(5.12) 

(5.13a) 
(5.13b) 

(5.13c) 



where Nc is the color factor, and d is the MS coupling defined in Sec. 2.4 (the decoupling of 
the top quark is irrelevant at the order we are working). In Eq.(5.13c) 



cz = j(5 - 3/3) - 3 - + 73 Q'j 
where 73 and Q are the isospin and electric charge of the target (73 



4 52 



lnc2-ln^ 1 In^ 



(5.14a) 

-1 for electrons), and 
(5.14b) 



^The mass counterterm for the unphysical scalars is related to the tadpole counterterm and therefore deter- 
mined from the tadpole diagrams [126]. 
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Figure 5.1: Vertex diagram contributing to 0{G'^mfMl). 

with ^ = M^/M^. Using Eqs. (5.10), (5.11), and (5.12), we can express p as follows: 

p = l + A/>(^' + N.xtAp^'^ + Ap^'\ (5.15) 
where the previous relation defines the two-loop contribution, Ap^^\ as: 

- (F;^ + 2 M^5;^). (5.16) 

Eq.(5.15) suggests that a possible way to take into account higher order eff"ects is to write p 
as 

P = (Tri^m) (1 + ^p'^'' + V^'), (5.17) 

where the resummation of A/?-^*^' can be theoretically justified on the basis of 1/Nc expansion 
arguments [119]. 

From Eq.(5.16) it is immediate to see that, unlike the case of the leading 0{G'^mf) 
contribution, the calculation of the self-energies is not sufficient to compute the 0{G'^mfM^) 
corrections. Indeed vertex and box diagrams do contribute to the order we are interested in. At 
the two-loop level the only box graphs that can give rise to terms are those obtained from a 
one-loop box diagram either by inserting a fermionic self-energy in a vector boson propagator 
line, or by expanding a bare coupling. The former can be taken into account, to the order we 
are working at, by replacing in the corresponding one-loop diagram the bare propagator mass 
by the physical one^. The latter do not give any contribution in the MS framework, because 
of the structure of the counterterms. For what concerns the vertex graphs, besides those with 

^In fact, the insertion of the on-shell mass counterterm cancels all mj terms coming from this kind of 
diagrams. 
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fermionic insertions in the one-loop diagrams as discussed above for boxes, one has to take 
into account the mixing between unphysical scalars and vector bosons in an internal line (see 
Fig. 5.1). 

The calculation of the vacuum polarizations is performed following the lines of Ref. [107]. 
In that paper two-loop self-energy diagrams were evaluated employing techniques used in 
the Current Algebra formulation of radiative corrections [44]. This framework provides an 
efficient way to enforce the relevant Ward identities while discussing at the same time several 
Feynman diagrams. As soon as subleading terms are taken into account, the relevant Ward 
identities do not have any more the simple form of Eq.(5.1) in a general gauge, and complicate 
considerably^. In the next section, as an illustration of the method, I will show how the Ward 
identities concerning the leading top contributions (Yukawa limit) can be derived using Current 
Algebra. The two-loop self-energy integrals are computed using the a Taylor expansion method 
outlined in Ref. [130,131], and described in App. A.2. As the transfer momentum is zero, the 
two-loop vertex integrals of the diagrams of Fig. 5.1 are equivalent to self-energy integrals, 
and can be solved accordingly. The calculation has been performed using an anticommuting 
75; we have checked explicitly that we never encounter traces involving 75 and more than two 
other gamma matrices or the antisymmetric tensor, so the inconsistency of this definition has 
no way to manifest in the actual calculation. The advantage is not only computational, as 
alternative definitions of 75 break the Ward identities of the theory. 

We define the renormalized top mass as the pole mass (see Sec. 4.1.1), using the countert- 
erm (in the 't Hooft-Feynman gauge^), Sm = m — mo, 

^ 16x2 M2, I eTm? j I 8e [ V 9 



with 




^A+zt[ """^^ I +0{zt2)\, (5.18a) 



zt , , , ht ht , , , , , ht - A r— ^ 

A = wtlnwt- A + —(1 +lnzi) + — + —(6 - ht)\a.ht Vhtg(ht). 

2 2 4 4 



I have used ht = (MH/mt) , zt = (Mz/rrit) , wt = (Mw/mt) , and 



(5.18b) 




(w - 2 arcsin y/xjtj < a; < 4 



9i^) = { , , , (5.19) 



Eq.(5.18) complements Eq.(4.12), valid for strong interactions, and it is obtained from dia 
grams like the one in Fig. 4.2a, where the gluon is replaced by scalars and vector bosons. 



^The calculation described here was performed in the 't Hooft-Feynman gauge; the Ward identities in App. C 
are valid in that particular gauge. In the framework of the Background Field Method it has been verified [120] 
that the Ward identities maintain the form of Eq.(5.1) up to 0(G^m( Af|). 

^Unless tadpoles are included, the pole mass counterterm is gauge dependent, see also Ref. [63]. 
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Explicitly, we find for A/)*"', in units Nc{d/{16TrPc'^)m^/Miy ~ NcX^: 
(-4 + ht) Vhig{ht) 



A/9('' = 25 - 4 M + 



-6-6 ht + — \ hiht 



"(^ ~ ~^ ~ ~^ ^"^^^ ~ - ht) + (9ht - 15 



3ht^\ ^.ht. 



zt 



25 



10 c' 



277 5^ 4 5^ 



2 ht i2 9 

3 (5-65^) In + 6 /5 + 1^2 - - 8 5 
0(1 2 1 565^ 



9+|r-i-6JM Inc^ 



2 2852\ 



ht 



2P P 



■K 



3 3 ht"^ ht 27 3 ht"^ ht J 

2 2 (-2 + 3/if) c^^ 



12 (-4 + 4 

/if ^ /ir 



2 /if 



/if^ 



1^2(1 - ht) 



for Mh > Mz, while in the region Mh < M^, 



A/)*'' = 19 - 2 - 4 X + ht 



27 
T 



zf 



11 3 319 5^ 



2 - 6 In /if - 5 In + 3 In zf 
^ 7 5652\ 



9 



6 /, c' + 



27 



(^7 + ^ - ^ - 45^^ Inc^ + (21 - 16 5^) Inzf 



In Eqs. (5.20) 



A(-l + -) 



-^g{x) + ^,/4/F^ 0<X<4: 



Li2{x) = - f 
Jo 



dt 



ln(l - t) 



and 



cPiz) 



< z < 1 



4^iZ7C/2(2arcsinyi) 

i [-4ii2(^) + 21n2(l^) - ln2(4z) + xV3] z > 1 
where Cl2{x) = Iin.Li2{e^^) is the Clausen function with 



A 



(5.20a) 



(5.20b) 



(5.21a) 



(5.21b) 



(5.21c) 



(5.21d) 
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The first two lines of Eq. (5.20a) represent the leading 0{G'^mf) result [115,116], which is 
completely independent of the gauge sector of the theory. The rest of Eq. (5.20a) is propor- 
tional to zt = M|/mj, and represents the first correction to the Yukawa limit. To completely 
identify what we call the 0(G^mjM|) contribution, we also have to provide a hierarchy in 
the couplings. Assuming gt ^ g,g\ where g' is the ^7(1) coupling constant, we have the two 
possibilities X> g (Eq.(5.20a)), or A < (Eq.(5.20b)). 

Eqs.(5.20) exhibit a process-dependent contribution, i.e. Qztl^c^, that comes from Bz^'^\ 
This reflects the fact that, already at one-loop, the box diagrams in neutral current depend on 
the process under consideration [37] (cf. Eq.(5.14a)). Because we choose to write our result 
in terms of the physical Z mass and c, we obtain a process-dependent term by re-expressing 
My^r, in the one-loop WW box diagram in neutral current, in favor of Mz using Eq.(2.19b), 
which at leading order reads 

Ml = Mlc^p c Ml (1 + ^) . (5.22) 

5.3 Ward identities by Current Algebra 

The local gauge symmetry of the SM Lagrangian induces relations between different Green 
functions, that are usually called Slavnov- Taylor identities [121], and generalize the Ward- 
Takahashi identities which derive from global symmetries. The use of the Becchi-Rouet-Stora 
(BRS) invariance of the SM Lagrangian makes their derivation very simple (see for ex. [122]). 
In particular, in a renormalizable gauge parametrized by ^ (the 't Hooft-Feynman gauge 
corresponds to ^ = 1), we find for the longitudinal parts of the unrenormalized gauge field 
propagators 

{d^A''{x)d^A^{0)) = -i^6\x), 
{d^A^{x)d,Z''{0)) + ^Mzd^A^{x)^2m = 0, 
{d^Z^'{x)d,Z''{0)) + eM^(9^Z^(a;)$2(0)) 

HMz{^2{x)d,Z''{0)) + eMl{^2{x)^2m = i^S\x), (5.23) 

where A^, Z^, and $2 are the photon, Z, and neutral Goldstone boson fields, and I have 
used the short-hand notation (...) for the vacuum expectation value (v.e.v.) of the covariant 
time-ordered product (0|T*...|0). The case of W is completely analogous. The identities 
Eqs. (5.23) make manifest the connection between the longitudinal degrees of freedom of the 
gauge bosons and their unphysical counterparts. As we are interested in vacuum polarization 
functions, we can extract identities for the self-energies of the corresponding fields at an 
arbitrary g^. At one-loop level, the result is essentially the one given in Sec. 4.3. At two-loop, 
however, the identities complicate considerably, and involve products of one-loop self-energies 
as well as irreducible VPF's, as can be seen, for instance, from Eq.(7.2) of Ref. [123]. Of 
course, at leading order in rrit, they reduce to Eqs. (5.1), but they are not of practical use in 
the case at hand, where the next-to-leading contributions need to be retained. An alternative 
is represented by the use of Current Algebra. 
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P; 

\ I 
\ I 
\ I 
\ / 

(e) 

Figure 5.2: Two-loop diagrams contributing to the vector boson self-energies at 0{G'^mf). 

In this section I present a Current Algebra derivation of Ward identities at = for the 
subset of diagrams contributing to leading 0{G'^mf). The same method can be applied to 
the remaining subleading contributions, and yields the identities reported in App. C. Just by 
power-counting inspection it is easy to realize that the only diagrams which can contribute to 
the leading mf term are those containing, besides top and bottom, the physical and unphysical 
scalars. A suitable choice of the tadpole counterterm allows us to neglect the tadpole diagrams 
in the calculation [126]. Consequently, the only topologies that contribute at leading order in 
rrit are those depicted in Fig. 5.2. In the figure, wavy lines represent vector bosons, dashed 
lines physical or unphysical scalars, while solid lines are fermions. 

In order to fix the notation, I write the part of the SM Lagrangian density that describes 
the interaction of the W, Z and scalars with fermions as 

£ir,t = -^{W^J!^ + h.C.) - ^Z^J^ - + $2^2 + V2 ($t5 + h.c.)] , 

(5.24) 

where and are the fermionic currents coupled to Z and W respectively, is the 
field that creates a W'^ meson, $i is the physical Higgs boson, $2 and $ the unphysical 
counterparts associated with the Z and W. 



51 =ip mP ip, 

52 =2dfj,Jz = -iil)'m° Csfsil), 



(5.25a) 
(5.25b) 
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id^J^ = 4>T^, (5.25c) 



are the hadronic currents coupled to the scalar fields. In Eqs. (5.25), ^l) represents the column 
vector tp = [t, 6)^, mP, C3 and T are the 2x2 matrices 

Cs = ( I M (5.26b) 







(5.26c) 



a_|_ = ^-y^, and the superscript on rrit refers to the bare mass. It is evident from Eqs. (5.26a) 
and (5.26c) that we consider only the third generation, and take the bottom quark as massless. 
We begin by studying Azz{0)- Using current correlation functions, and working in n dimen- 
sion, we can combine the amplitudes of Figs. 5.2a,b, where the continuous line represents a 
top and the dashed one a Higgs or $2? into the expression 



3 



X [ e-'-^y [ e'^-^^ f e-^'=-^^(0|T*[J^(t/)J^(0)5,-(a;i)5,-(a;2)] |0), (5.27) 

with the short-hand notation = J c?"fc/((2x)"/x"~^) and = / d"-x, where /x is the 't Hooft 
mass scale, T* is the covariant time-ordered product and mi = and m2 = M^. In the 
case of the unphysical scalar, Eq.(5.27) is valid in the 't Hooft-Feynman gauge. 

In order to trigger Ward identities we contract n^^j^^(g^) with . Contraction of a 
current operator with its four-momentum gives rise to a term involving the divergence of 
a current plus an equal time commutator that reduces the number of operators inside the 
time-ordered product by one unit. Noticing that 

^ {m^^W)}^,^^ = AzziO), (5.28) 



we obtain 



+ i [ e"' ^S2{y) Si{x) 

J X 

16 c-'Ml jk\\ nk^-mlj Jo, {k^ - m\f 



— ml (k — qY - 




(5.29) 



CHAPTER 5. HIGHER ORDER ELECTROWEAK CORRECTIONS 



76 



In Eq.(5.29) the notation (1 2) represents a term obtained by the previous one inside the 
curly bracket by the substitution 1 2. 

We now examine the contributions that are described by the topology shown in Fig. 5.2c. 
Let us consider the case when one dashed line represents a $i and the other one a $2- We 
write 



9 



Contraction with gives 



J k Jy J X 



[k^ - ml][{k - qy - ml] " 



(5.30) 



{2k - qf 



ScmiJk [k^ - ml][{k - qy - ml] 

X 1/ e-'-i-y [ e"'-^S2{y) Sr{x) + i [ [e^M-^(52(a;) 52(0)) 

-e^'=-(5i(a;)5i(0))]}. (5.31) 



Recalling that 



(5.32) 



it follows that the contribution of Fig. 5.2c to 1122(0) is obtained by differentiating Eq.(5.31) 
with respect to q^ and then setting = 0. Consider now the first term in Eq.(5.31), 



Jk 



{2k - qf 



[k-^ - ml][{k - qf - ml 



{ - } 



(5.33) 



where |-"| represents the three-point correlation function times the appropriate constants. 
By Lorentz invariance, T^ should have the form T^ = a{q'^) q^ , and therefore 



The contribution to A{q'^) is then given by Q;(g^). It is easy to show that 

d 



«(0) = ^2^1^Tn 



(5.34) 



(5.35) 



Using Eq.(5.35), we can write for the transverse part of 11' 



^(c)(0) 

Jk 



d 



8 c'^Ml dq^ 
1 



1/ e-'-i-y [ e"'-^S2{y) Sr{x)S2{0)) 



m\ — ml 



{k — qY — ml k"^ — m\ {k"^ — m\)[{k — qY 



mt 



(5.36) 



n 



k^ 



— / 

ml I Jx 



(52(^^)52(0)) 

(P-m2)(P-m2) 



(1-2) 
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The contribution of $i and $2 to Figs. 5.2d,e can be similarly computed 

,4 . r / 4 ^2 



Ad,e(0) 



9 



16 c^Ml Jk 



1 



(P -m2)2 



(5.37) 



The diagrams involving the charged scalar $ can be considered along the same lines, and give 
^(«-e)(0) = J^Ml^ [Ik k^-Ml 

X [ e-'-^y [ e'^-^^ j e-'^-^^{S2{y)S2{^)S\xr)S{x2))\ . (5.38) 

Jy Jxi Jx2 ' q^=0 

Finally, summing Eqs. (5.29), and (5.37)-(5.38), we find that the subset of diagrams defined 
above (see Fig. 5.2) contributes the transverse part of the Z self-energy at = a factor 



lead 



(0) d 



Ml 



ik-X2 



{S2{y) 82(0) S^{x^) S{x2)l ^ (My) 52(0) Sj{x^) Sj{x2)) 



4(P _ m}) 

-i(M' -M') I e-'-i-y I e^^- (52(t/) S^jx) 82(0)) 1 



) f e-^-i-y f 

Jy J X 

16M^ Jk (fc2 - M2)2(fc2 - M|)2 



1 _ 1 



ft TTt^ I Jx 



lj"'-^Sr{x) Sr{0)) + {1 ^ 2) 



(5.39a) 



The discussion of "leading" diagrams of the W self-energy can be performed on the same 
footing as in the Z case. Here I present only the result: 

Ay^J'-'^O) _ d 



Ml 



SM^ dq^ 



J k \.Jy J xi J X2 



ik-X2 



{S\y)S{0) {s\x^)S{x2) + h.c.) ) 2 ^^t(j/) 5(0) Sj{x^) Sj{x2)) 



J y J X 



(P _ Ml) ' ^§ 2(P _ 

((5t(t/)5(0) + h.c.) S^{x)) 



i{Ml - Ml) 
{Ml - Ml) 



(P - M2 )((fc - g)2 - Ml) 
((5t(t/)5(0)-h.c.) S2{x)) 



{k^-Ml){{k-q)^-Ml) 



/ 

Jk 



{ml - Ml)' 



16M4, A 1 (P-m2)2(p-M2,)2 



nfc2 _ M2, , 



/ e"'-^{Sr{x)Sr{0)) 

J X 
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(a) 



(b) 




(d) 



Figure 5.3: Two-loop diagrams contributing to 11$$ and to n$2$2 in a Yukawa theory. 



2 1 



nkP' — rrii 



\ I Jx 



(1-2) 



(5.39b) 



In order to understand the connection between Eqs. (5.39) and the r.h.s. of Eqs. (5.1) 
consider for example the four-point correlation functions in Eq.(5.39a). In the limit g, g' ^ 
or Myy, Mz 0, their contribution to Ap is proportional to gf. In fact each S operator 
contains an rrit = gtv factor, where v is the v.e.v., while the coefficient in front is proportional to 
l/w^. In this limit the terms involving four-point functions represent exactly the contribution 
to ^2n$2$2(?^) 2 of the diagrams shown in Figs. 5.3a and 5.3b. A similar connection 
can be established between the three- and two-point correlation functions in Eq.(5.39a) and 
the diagrams 5.3c and 5.3d respectively. We notice that Figs. 5.3c and 5.3d involve trilinear 
scalar couplings proportional to M^. Although no such coupling is present in the diagrams of 
Fig. 5.2, we recover these terms using Eq.(5.35) and simple algebraic manipulations. Equation 
(5.39b) shows an analogous connection between its various contributions and ^2n$$(g^) ^ 
computed in the Yukawa theory. 

We have explicitly verified the Ward identities Eqs. (5.39), which represent an important 
check of the calculation. It might be worth emphasizing that Eqs. (5.39) establish a correspon- 
dence between a whole class of graphs evaluated at = 0, and no additional approximation 
has been applied. They are valid for the first term of the heavy top mass expansion (Eq.(5.1)), 
as well as for the next terms of the expansion; their physical content is therefore diff"erent from 
Eq.(5.1), which they generalize. Obviously, there are other IPI diagrams that can contribute 
to subleading order. Ward identities dealing with them are considered in App. C. 
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Figure 5.4: A/?*^' for Ui^^e scattering, in units NcX^ as function of rrit for few values of Mjj- 
including only the 0{Glmf) contribution (y), and with both the 0{Glmf) and 0{GlmfMl) 
terms (g). 



5.4 Scheme dependence and numerical results 

In Sec. 5.2, the calculation of Ap was performed in the MS scheme, expressing the one-loop 
contribution in terms of a, Mz, and (with the exception of the top one-loop contribu- 
tion, which is conveniently expressed in terms of Gfj,, as we have seen). I present now the 
corresponding expression for A/?*^' in terms of a, Gfj,, Mz and the on-shell (OS) parameter 
= Ml /Ml. The relations [11] 

a _ G,Ml __ G,My 
4xJ2 2V2x2 1 + (^)^ 2^2x2 ^ • ' 

c2 = c\l - Yj^) ~ c\l - N^xt) (5.40b) 

allow us to replace the MS quantities d, P by Gfj, and c^. Furthermore, Eq.(5.40b) creates 
additional contributions to Ap^^K The one-loop result is then given by Eqs. (5.13) with the 
substitutions a/(4xs^) (G^M|c^)/(2v^x^), and s^, s^, c^, while for the two-loop 

contribution we have 



(5.41a) 
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Ncxlzt 



s " 



where 



2 



(c2 - c2 - ^ ' c2 1 - ^ 



(5.41b) 



In Eq.(5.41a) A/9*^'(i2^ c2 —i s"^, c"^) represents a term obtained from Eqs. (5.20) applying the 
same substitutions as in the one-loop case. 

Fig. 5.4 shows A/?*^' (Eqs. (5.20)) as a function of rrit for different values of M^- As a 
comparison, we also show the values obtained including only the 0{G'^mf) contribution. The 
process under consideration is i/fj^ e scattering. From the figure it is evident that the inclusion 
of corrections suppressed by a factor M|/m2 with respect to the leading term has a quite 
substantial eff"ect. For v — v scattering the eff"ect is even more pronounced. To have a better 
understanding of the size of these eff"ects and the diff"erence induced by the choice of the 
renormalization scheme, I present in Table 5.1 the values of A/?*^' and A/Oq^ for zt = 0, and 
for zt = 0.2,0.3, corresponding to rrit =~ 200, 170GeV, as a function of r = Mu/rrit. In 
preparing the table we matched the values coming from (5.20a) and (5.20b) where they meet. 
We see that in the region of light Higgs the 0{G'^mfM^) corrections are much larger than 
the mf term, that is actually suppressed by accidental cancellations, while for large Higgs 
mass, in the TeV region, their contribution is still 50% of the leading part. It is intriguing to 
note that the numbers shown in Table 5.1 are very close to the corresponding ones obtained 
in Ref. [107] in the case of a model with SU{2) symmetry, where g' has been set to zero. 

From Eq.(5.41a) we see that the process-dependence is more pronounced in the OS frame- 
work. This is easily understood by noticing that the expansion of the bare couplings in the 
one-loop box diagrams gives rise, unlike the MS case, to contributions. It is worth mention- 
ing that in the MS framework the process-dependence can be confined entirely in the one-loop 
contribution. This can be achieved by expressing the one-loop vertices and boxes using MS 
couplings and the physical W and Z masses whenever they appear in the propagators. Such 
a procedure is frequently used in one-loop calculations [11]. A numerical investigation shows 
that this procedure minimizes A/?*^', for any value of M^jfrit. I present the values of A/?*^' 
for this particular MS prescription in Table 5.2 in the case rrit = 180GeV, and I will use them 
for the estimates of the next section. Table 5.1, instead, has been obtained following Sec. 5.2. 
In Fig. 5.5 I plot A/?*^' for i/fj^ e scattering as a function of the ratio between and rrit, ex- 
pressing again the one-loop result in terms of MS couplings and physical vector boson masses. 
For a comparison, the upper curve shows only the leading 0{G'^mf) term. It is interesting to 
note the very large diff"erence in the light Higgs regime. 

Our investigation of the two-loop next-to-leading top corrections to the p parameter reveals 
several interesting features, i) The 0{G'^mtMl) corrections are not universal but depend upon 
the physical process under consideration. Indeed such terms can originate from self-energies 
as well as from vertex and box diagrams, which are process dependent, ii) The asymptotic 
result, obtained in a Yukawa theory, does not seem to be a realistic approximation for values 
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MS 


OS 


TTlt 


zt = 


zt = 0.2 


zt = 0.3 


zt = 0.2 


zt = 0.3 


0.1 


- 1.8 


-12.6 


-15.8 


-12.7 


-16.0 


0.2 


- 2.7 


-13.3 


-16.5 


-13.5 


-16.8 


0.3 


- 3.5 


-13.9 


-17.0 


-14.2 


-17.4 


0.4 


- 4.1 


-14.5 


-17.6 


-14.9 


-18.1 


0.5 


- 4.7 


-15.2 


-18.3 


-15.7 


-18.9 


0.6 


- 5.2 


-15.8 


-18.9 


-16.6 


-19.7 


0.7 


- 5.7 


-16.2 


-19.8 


-16.9 


-20.9 


0.8 


- 6.2 


-16.4 


-20.1 


-17.1 


-21.0 


0.9 


- 6.6 


-16.5 


-20.1 


-17.4 


-21.2 


1.0 


- 6.9 


-16.6 


-20.1 


-17.6 


-21.3 


1.1 


- 7.3 


-16.8 


-20.2 


-17.8 


-21.4 


1.2 


- 7.6 


-16.9 


-20.2 


-18.0 


-21.6 


1.3 


- 7.9 


-17.0 


-20.2 


-18.2 


-21.7 


1.4 


- 8.2 


-17.2 


-20.3 


-18.4 


-21.9 


1.5 


- 8.4 


-17.3 


-20.3 


-18.6 


-22.0 


1.6 


- 8.7 


-17.4 


-20 4 


-18.7 


-22.1 


1.7 


- 8.9 


-17.5 


-20.5 


-18.9 


-22.3 


1.8 


- 9.1 


-17.6 


-20.5 


-19.1 


-22.4 


1.9 


- 9.3 


-17.7 


-20.6 


-19.2 


-22.6 


2.0 


- 9.5 


-17.8 


-20.6 


-19.4 


-22.7 


2.5 


-10.2 


-18.2 


-20.9 


-20.0 


-23.3 


3.0 


-10.8 


-18.4 


-20.8 


-20.4 


-23.5 


3.5 


-11.2 


-18.3 


-20.6 


-20.6 


-23.6 


4.0 


-11.4 


-18.3 


-20.4 


-20.6 


-23.5 


4.5 


-11.6 


-18.2 


-20.1 


-20.6 


-23.4 


5.0 


-11.7 


-18.0 


-19.8 


-20.5 


-23.3 


5.5 


-11.8 


-17.8 


-19.4 


-20.4 


-23.1 


6.0 


-11.8 


-17.5 


-19.0 


-20.3 


-22.9 



Table 5.1: A/)*"' (MS) and Ap^ol {OS) relevant to p^e scattering for zt = Ml/m^ = 0.2,0.3, in 
units NcXf as a function of r = Mu/rrit. The column zt = \s the result of the Yukawa theory. 
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0.1 


-11.3 


1.5 


-15.8 


0.2 


-11.9 


1.6 


-15.9 


0.3 


-12.5 


1.7 


-16.0 


0.4 


-13.0 


1.8 


-16.1 


0.5 


-13.6 


1.9 


-16.1 


0.6 


-UA 


2.0 


-16.2 


0.7 


-15.2 


2.5 


-16.6 


0.8 


-15.3 


3.0 


-16.6 


0.9 


-15.3 


3.5 


-16.5 


1.0 


-15.4 


4.0 


-16.3 


1.1 


-15.5 


4.5 


-16.1 


1.2 


-15.6 


5.0 


-15.8 


1.3 


-15.7 


5.5 


-15.5 


1.4 


-15.7 


6.0 


-15.2 



Table 5.2: A/?*^' (MS) relevant to p^e scattering for zt = Ml/m^ = 0.2,0.3, in units NcX^ as a 
function of r = Mu/rrit. 

of the top mass compatible with the experimental constraints. In this range, the first order 
gauge corrections are numerically comparable to the 0{G'^mf) contribution. The values of 
for which A/?*^' is well approximated by the 0{G'^mf) contribution are very large. Typically, 
A/o*"' starts to be within 10% of the leading mf value for rut ^ 800GeV. iii) The 0{Glm^Ml) 
contribution has the same sign of the leading term and enhances its effect. 

5.5 Extrapolation to and sin^ ^'^^ 

The result illustrated in the previous section, obtained at = 0, cannot be directly applied to 
LEP physics. We have seen that there are terms 0(G^mjM|) that depend on the specific pro- 
cess under consideration, and we know that Ar and Ar involve two-loop vacuum polarization 
functions evaluated at = M| or M^. They differ from the self-energies at = by terms 
of 0(G^mjM|). However, the size of the correction that we have found, relative to the leading 
0{G'^mf) term, suggests to attempt a quantitative estimate of the theoretical error that we 
can make when we include the leading corrections, but neglect the ones of 0(G^mjM|). A 
complete calculation of this kind of effects in Ar and Ar, the relevant radiative corrections in 
the prediction of M^i^ and sin^ ^e//' imder way. Meanwhile, we can first notice that for rrit in 
the current range the result of the previous section is not a complete surprise. Contributions of 
next-to-leading order are expected to be roughly (d/xs^)^ mj/M| ln(mj/M|) ~ 6 X 10~^ for 
rrit = 180GeV, of the order of magnitude of the uncertainty due to the hadronic contribution 
to Aa. 
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Figure 5.5: A/?*^' for Uij^e scattering, in units NcX^ for rrit = 180GeV as a function of Mu/mt: 
including only the 0{G'^mf) contribution (upper curve), and with both the 0{G'^mf) and 
0{G'^mfMl) terms (lower curve). The one-loop bosonic contribution is expressed in terms of 
MS couplings and physical vector boson masses. 



More specifically, we can assume that the ratio between the ©(G^m^Mj) and the ©(G^m^) 
contributions in A/?*^' be representative of the unknown two-loop top effects in Ar and Ar. 
We can then use this ratio to estimate the additional contributions to Ar and Ar simply 
multiplying it by the known 0{G'^mf) terms of these quantities. The shifts in the W mass 
and the MS mixing angle, sin^ 9w{Mz) (and consequently in sin^ ^e//^' these additional 

contributions can be estimated using Eqs. (2.15) and (2.21) once more. In Table 5.5 we show, 
for a few values of rrit and M^, the effect of our estimate of the unknown top contributions 
on the W mass and sin^ ^e//' have assumed that Eq.(2.27) is not changed by effects 
of 0(G^mjM|), which we can decouple from the definition of sin^ dw{Mz), and that the 
implementation of the decoupling does not give relevant 0{G'^mfM^) effects in Ar. The ratio 
of subleading and leading terms in A/?*^' has been computed using expressions slightly different 
from Eqs. (5.20). In fact, as mentioned in the previous section, we decided to minimize the 
two-loop corrections of our MS calculation by writing the one-loop contribution in terms of 
the MS couplings and of the physical masses of W and Z. Such a procedure [11] has the 
further advantage of eliminating the process-dependent terms. From the third column, it can 
immediately be seen that, for a fixed value of the top mass, the effect is more pronounced for 
a light Higgs boson. This is not surprising, bearing in mind the fact that the 0{G'^mf) term 
is a monotonically increasing (in modulus) function of M^- 

It is important to stress that the numbers presented in Table 5.5, more than a definite 
estimate of the shifts in and sin^ ^eT/' should be taken as an indication that subleading 
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mt 
(GeV) 


M« 
(GeV) 


R 

% 


(MeV) 


(10-4) 


150 


65 


247 


-10 


0.6 


250 


100 


-8 


0.5 


800 


35 


-4 


0.2 


175 


65 


234 


-16 


0.9 


250 


94 


-14 


0.8 


800 


38 


-8 


0.5 


200 


65 


221 


-23 


1.4 


250 


88 


-20 


1.2 


800 


38 


-13 


0.7 



Table 5.3: Calculated ratio (R) of the 0{Glm^Ml) and the 0{Glmf) contributions in Ap^^\ for 

a few values of rrit and Mjj. The corresponding estimate of the shifts in the W mass and sin^ ^^^^ 
are also presented (see text). 



two-loop mt effects could be larger than what is naively expected. Their size is probably 
comparable to, and possibly larger than, the theoretical uncertainty due to the hadronic 
contribution to the photonic self-energy. As we have seen, the latter amounts to 13MeV and 
2.5 X 10-4 in Mw and sin^ ^e//' respectively. 

In conclusion, the calculation of A/?*^' presented here shows that, at least in the specific 
case of the p parameter for neutrino-lepton scattering, two-loop electroweak top contribu- 
tions are not well approximated by the 0{G'^mf) term. It seems reasonable to expect a 
similar behavior in the more relevant cases of high-energy precision observables like Myy and 
sin^ ^e// ■ -Although at present this kind of effects are well beyond the experimental resolution, 
the precision foreseen in the next decade (see Table 2.3) seems to justify a theoretical effort 
to investigate the reliability of the current analyses. In other words, a complete two-loop 
calculation for the most important precision observables may be worthwhile. 



Chapter 6 

Conclusions 



In the preceding chapters I have focussed on the theoretical determination of the two precision 
observables and sin^ ^e//' ^^^^ appear to be most promising for a significant improvement 
of the experimental precision in the next decade. I have analyzed the major sources of uncer- 
tainty in these determinations, and tried to estimate their impact, in relation to the expected 
experimental accuracy. In doing so, I have also presented original calculations that, as in the 
case of QCD corrections, find applications in other areas of precision physics. A few points 
can now summarize the results of the investigations presented in the various chapters. 

• The whole set of one-loop radiative corrections is needed to describe present data. 
Current experiments are beginning to probe radiative corrections beyond the leading 
fermionic contributions, and to become sensitive to the symmetry breaking sector of the 
theory. Even if this may not be refiected in a high sensitivity to the mass of the Higgs 
boson, the comparison with diff"erent models (supersymmetric SM, diff"erent symmetry 
breaking mechanism, etc.) will certainly be extremely valuable. 

• At the moment, the uncertainty on the hadronic contribution to the running of a is one 
of the major obstacles to an improvement of the theoretical error. New low-energy data 
can help to solve this problem, but a clarification of diverging analyses [13,22,23] is also 
needed. 

• For what concerns the QCD corrections entering the prediction of and sin^ ^e//' 
we can safely argue that they are well under control. I have provided a set of explicit 
formulae that allow the computation of the universal QCD corrections to a number 
of processes, involving vector and scalar bosons, for any value of the transfer momen- 
tum and masses. The leading term proportional to is known up to O(a^) and the 
theoretical error can be estimated, for instance, according to Ref. [51]. 

• As shown in the specific case of the p parameter at low-energy, next-to-leading 0{G'^mfMl) 
corrections to electroweak observables can be relatively sizable, and the leading 0{G'^mf) 
term does not provide a good approximation of the two-loop electroweak corrections. 
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An extrapolation to high-energy observables indicates the need for an explicit calcula- 
tion of 0{G'^m^Ml) effects on the theoretical determination of and sin^ ^e//' 
the moment, it seems appropriate to consider this extrapolation as an estimate of the 
theoretical error associated with higher order electroweak contributions. 

• A very recent work [137] confirms this point of view: a comparison between five different 
computer programs for the evaluation of radiative corrections, based on very different 
theoretical approaches, indicates discrepancies up to ~ 25MeV and ~ 2 X 10~^ for 
Mw and sin^ ^e//' respectively, in the current rrit range. These uncertainties or scheme 
ambiguities are of the same order of magnitude of the next-to-leading contributions I 
have discussed, and may be interpreted as indications of sizable higher order effects. 

• Among the higher order contributions that I have considered, it appears that only 
the QCD corrections to the one-loop quadratic top contribution may be relevant at 
the present level of precision. However, if the performance of the experiments already 
planned will be at the level of expectations, a meaningful theoretical interpretation will 
require a complete study of two-loop effects for the most relevant precision observables. 
After all, it should be kept in mind that the accuracy of the Z-resonance data at the 
end of the LEP program has abundantly surpassed all the expectations of a few years 
ago. 



Appendix A 



Two-loop integrals 



A.l Scalar integrals for 0{aas) calculations 

In this appendix I provide the expressions of the one and two-loop scalar integrals which are 
needed in the evaluation of the two-loop vacuum polarization functions discussed in Ch. 4. 

It has been shown in Ref. [123] that it is possible to express any two-loop VPF's exclu- 
sively in terms of one and two-loop integrals involving only products of scalar propagators. 
In other words, the numerator of two-loop integrals involved in the calculation of electroweak 
self-energies can always be cancelled against appropriate combinations of the factors (p? — 'ml) 
appearing in the denominators. This is not apparent from the two loop integrals reported here, 
some of which do involve integration momenta in the numerator, but can be easily understood. 
In fact, any numerator that cannot be cancelled against the corresponding denominator cor- 
responds to a vector or tensor sub-integral in one of the two variables of integration. These 
one-loop sub-integrals can be further reduced to one-loop scalar sub-integrals according to a 
well-known procedure (see for ex. Ref. [28]), reducing the whole integral to a sum of two-loop 
scalar integrals. This is the case, for instance, of the integral 



This general decomposition is not possible, however, in the case of two-loop vertices and boxes, 
for which the tensor structure is more involved. 

A. 1.1 Relevant one-loop integrals at 0{e) 

In two-loop calculations one-loop integrals are needed at 0(e). Products of two disconnected 
one-loop integrals can be found as a result of the above mentioned decomposition in scalar 
integrals, and of the insertion of counterterms in the first order diagrams. Following the 
notation previously set (n = 4 — 2e, a = —m'^/q'^, etc.), we have: 




(Al) 



1 



e'^ r(e - 1) 



^n/2 Jg2 _ ^2 



= a 
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-a 



— hi— log a + e j — log^ a — log a + 1 
6 \ 2 



(A2) 



K 



' J X"/2 



'/2 [k^ - ml][{k - qy - ml] 



2 + (a - /3 + A2 ) log a;,! - log a 



[a ^13] + 0(e). 



(A3) 



The 0(e) part of if can be found in Ref. [128]. Note that in the case of the calculations 
0{aas) discussed in Ch. 4 the 0(e) term of K is not needed: its total contribution, via terms 
like K^, KDa, KdK/da, • • • , vanishes in the final result. The derivatives oi K and Da with 
respect to a enter sometimes in the expressions of two-loop integrals and are given by 



dDj 
da 
dK 



2/1, 2 

e log a + e — log a -\ 

^ I 2 ^ 12 , 



(A4) 



= ^ [(l + a-/3 + A2)loga;„ + (l + a-/3-A2)loga;J + 0(e) 
da 9XJ L J 



2A 

The following relations are useful for numerical evaluation of the VPFs; in the equal mass 
case, for D = —1 — 4a > 0: 

1 

2" 



.£,1/2 



arctan 



X)-l/2 



where x 



where 



4ce 



(l+VT+4^^)2 



^Vl + 4a In X 
In the arbitrary mass case: 

Re [A2(lna;„ + In a;;,)] = -40(a,/3) 



(A5) 



(A6) 



0(a,/3) 



arccos 



l+a+/3 



4 In 



l+a+/3+A2 
1 

l+a+/3-A2 



A < 0, 



A > 0. 



(A7) 



A. 1.2 Scalar two-loop integrals. 



All the two-loop integrals needed for the calculations described in Ch. 4 can be expressed in 
a compact form in terms of polylogarithms of algebraic functions of the masses and external 
momentum. This is far from being a general feature of two-loop scalar integrals, and it has 
been noted [67] that this is the case only if all three particle cuts of the diagram satisfy the 
condition 



m^m^ml Il± — (rrzi ± m2 ± ^3)^) = 0, 



(A8) 



where are the three masses involved, p the total momentum going across the cut, and 
the product is over all the combinations of signs. As the gluon is massless, all the two-loop 
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integrals used in the calculations of Ch. 4 satisfy this condition. To simplify the notation, I 
use for the definition of the two-loop integrals the following abbreviations: 



^ ' J J 



(A9) 



and 



iri,2 = -ml, Ko = {h - k^f , (3i,2 = (fci,2 - qf - ml (AlO) 



With the help of the variables defined in Eqs. (2.9)-(2.10), the two-loop integrals are 



Wa 
Ta 



Ra 
Na 



< 



q^K^K2QiQ2 



>= (1 - a + I3)KDa + (1 + a - I3)KDb 



-^{Da-Db)'-\\K'; 



1 1 



-Da 



da ' 



Ua = < 
Ma 



KlK2Qi 



>-- 



1 

^ < 



q 

Va = < 



2 ^ KiK2Ko 
1 



-Da^; 
oa 

1 1 - 



al - 2e 
1 1 - e 



n2 • 



KIK2K0 



>-- 



1 

1 (fci - 2k2) ■ q ^ 



-Da (da - 2a^^^ : 
2a2 1 - 2e V da J ' 

>= Ma + {1- a + I3)Va; 



q 

Ja = < 



KiK2KoQi 
1 

> 



\ 1 ^OgXa' 




1 ajfeloga;;," 






1 - iCb . 



a;alog^ Xg _ a;;, log^ Xb 

{l-XaY ~ {l-XbY 



G{xa) + G{xb) 



< 

1 
2 
1 

^ < 



q'KrQ2Ko 



> 



-(2 + e)ir2 + + iV^ - (1 - a + /3)J^ 



>= ^[(l-a + /3)J^-i + M^]+iV^; 



KiK2Q2Ko 



Sa = < 



KIK2Q1K0 
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/ ^ < ^ >= F{1) + F{x^x,) - F{x^) - F{x,). (All) 

The functions F[x) and G{x) are given in Eq.(4.22). Some of the preceding expressions can be 
obtained from one-loop integrals and the "master" integral I by differentiation with respect to 
the masses, and exploiting symmetries in the shift and interchange of the integration variables 
ki and k2, as well as of the masses nia and m;,. Others can be found in Ref. [134,55]. All 
the calculations described in Ch. 4 can be performed on the basis of this (non-minimal) set of 
integrals, using again symmetries and differentiation with respect to the masses. The master 
integrals I has been calculated by Broadhurst in two different ways in Ref. [133]. 



A. 2 Small momentum expansion of VPF's 

The two-loop self-energy integrals necessary for the calculation of Sec. 5.2 are evaluated at = 
0; the use of Ward identities sometimes introduces derivatives of two-loop integrals at = 0, 
and they all involve many different masses. In the approximation that the bottom is massless, 
there are at least four different scales entering that calculation: M^, Myy, M^, and m^. The 
number of relevant Feynman diagrams is also quite large. It is therefore useful to develop 
an algorithm that can automatically compute the Taylor expansion of two-loop integrals in 
the squared momentum transfer [129-131]. For the case at hand we can restrict to the 
first two terms of such expansion. The coefficients of the momentum expansion correspond to 
two-loop vacuum integrals and can be computed exactly, for any value of the masses involved, 
and expressed in terms of polylogarithmic functions. 

Formally, we can write the Taylor expansion around g^ = as 

where Dg = d"^ /dqadq"', and [a)j = T[a + j)/T[a). In practice this is equivalent to expanding 
each denominator containing the external momentum according to 

1 1 / 2k-q-q^ 4(k- qf ^\ , . , 

' 1 ' ^ ^ ' ^ ^' ' O(g^) . (A13) 



(fc — g)^ — k^ — m? \ k^ — m? (fc^ — rn^Y 

A rigorous theory of mass and momentum expansions of Feynman integral can be found in 
Ref. [132]. 

As mentioned above, the coefficients of the momentum expansion are vacuum integrals, 
i.e. integrals with zero external momentum. They can always be algebraically reduced to the 
form 

1 

'(fc2 - )"(p2 - m|)^((fc - pY - mlY 



where a, 6, c are integers. Analytic expressions for any value of {a, 6, c} can be obtained 
[130] in terms of Appell's hypergeometric functions of two variables, but a more convenient 
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representation in terms of simple polylogarithms follows from the realization that any integral 
with a given set of {a, b, c} can be further reduced to a combination of integrals with a,b,c < 
1. Integrals with all exponents equal to one can be expressed for any mi^2,3 in terms of 
dilogarithms or Clausen functions [130]. Integrals with one of the exponents equal to zero 
correspond to the product of two disconnected one-loop integrals and can be trivially obtained 
from Eq.(A2). This reduction can be achieved equivalently by the use of recurrence relations 
based on integration-by-parts technique [130], or by differentiation of the basic a = b = c = 1 
integral with respect to the masses. We have checked that the two approaches give the same 
result. Note also that vertex diagrams at g2 = can be treated by the same method, and 
effectively correspond to the same kind of integrals. 



Appendix B 

QCD corrections: useful formulae 



B.l Asymptotic expressions 

The expressions obtained in Ch. 4 for the VPF's of electroweak vector bosons and scalars 
are as general as possible, but in many applications it is sufficient to know some asymptotic 
behavior of the general formulae. As the asymptotic expansions are usually not elementary, 
in this Appendix I give a list of useful expressions for the most relevant limiting cases: equal 
quark masses, one massless quark, and high momentum transfer. The expressions for vanishing 
external momentum in the case of vector boson self-energies have been discussed in Sec. 4.1.3. 
The notation will follow Ch. 4. 



B.1.1 Vector boson self-energies 

From the one-loop VPF of Eq. (4.6), it is easy to obtain for m;, = nia 



4a 



--/9„+--4a+(l + 4a)^/2(l-2a)ln— , 



n^(.) = 

n^(5) = 25a 



(1 - (1 + 6a) + ^ + 8a + (1 + .aff^ In J^^^^=y, 
- - Pa + 2 + {l + 4aY/^ln 



(1 + vrT4^)2 



and of course in this case 11^ = 0. For m;, = one has 



2 + 3a --yoj + — + 2a-a2+ 2 + 3a-a^ In— — 
Ve / 3 1 + a 



L [s) = -a 



^ - Pa + 2 + a + {1 + aflnj^ 



a 



(Bl) 



(B2) 



At the two-loop level, with all masses defined on-shell, the asymptotic expressions allow 
several checks of the general formulae of Eq. (4.20), as the VPF's have been known for some 
time in particular cases. Note that here the color factor iVc = 3 is included, unlike in the 
one-loop expressions. 
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Case rrib = rUa 

The real parts are given by 

r 1 55 26 2 

(5) = 5 I — - + — - y a + Vl + 4a(l - 6a) In a; - -a(4 + a) In^ x 

+ ^{4a'-l)[F{l) + F{x')-2F{x)\ 

- ^(1 - 2a)7rT4^ [g{x') - G{x)] } , 

f 6a 1 55 19 

n^(5) = 5 |-— + (1 + 24ap - 22a)— - (1 + 12a/> - 22a)p + " y« 

2 

+ 4a^ + (1 + 12a + 4a^)Vl + 4a In a; + -a(5 + 11a + 6a^)ln^ x 

3 

- ax^ - -(1 + 2a)(l + 4a) [f(1) + ^(a;^) - 2F(a;)] 

- ^{l + 4af/' [g{x') - G{x)]Y (B3) 

{3 / 1 1 \ 1 3 
— ^+ 6/9 - + --6a + 3(3- 2a) Vl + 4a In x 

+ (11 - Qp)p + (3 + 4a - 6a2) In^ x 

- 2(1 + 2a) [f(1) + F(a;2) - 2F(a;)] - 4^1 + 4a [^(a;^) - G{x)\ } , 



with X = 4a/(l + + 4a)^, and the functions F(x) and G{x) have been introduced in Chap. 
4. The vector part of the longitudinal component vanishes in this case: this is expected to 
occur as a consequence of a QED-like Ward identity. 

In the equal mass case, 11^(5) coincides, up to a color factor, with the irreducible part 
of the photon self-energy in QED which has been calculated in the fifties by the pioneers of 
Ref. [58]. IIt'^{s) has also been derived using a dispersive approach in [56], and more recently 
in [57]. In fact, we have expanded the expressions of the previous sections around m;, = nia 
and m,}) = retaining terms up to order (m;, — rUaY j s and rn^l s, respectively. The rather 
lengthy expressions [136] will not be given here. 

Case rri}) = 

For one massless quark, the coefficients of n^j,(s) vanish, so that n^j,(s) = n^j,(s) = 
s n^j,(s), and one obtains for the real parts 

- — a+(l + a) h +-a- -a^ j Ina; + -(1 - 3a)(l + a) In^ a; 

- ^a' + 1(2 + a)(a - l)G(a.) + i(a - 2)(1 + a)' [F{1) - F{x)]^ , 
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2e^ 



11 



13 7 11 



e 8 2 



^(1 + af Ina; + (3 + 2a)(l + af^^ + (1 + a)G{x) 



{1 + af [F{1) - F{x)]} , 



(B4) 



with X = a/(l + a). The expression of Il^'^(s) in this special case has been obtained in 
Ref. [57] by means of a dispersion integration; the connection with Kniehl's result is discussed 
in Sec. 4.4. 



Case s > ml f^ 

From Eqs. (4.20), in the limit |s| ^ j^, one obtains for 11^,'^ (-s) 

1 , -5 55 



oo) 



2e 



In- 



12 



4C(3) 



_3_ 
2^ 
9 



^ ( Y - 3/>„ - 3/>, 



11, , 11 

— [Pa + Pb)- Y 



-{Pa + PhY - T(lna + ln/3) - - In a In /3 + 6C(3) + — 



+ (m^-m,^) In -f 



-^ + ^(/'a + /'6)-^(lna + ln/3) + ^ 



3(m2 +mfe2)(lna + ln/3), 



(B5) 



n 



V,Aj 



^ / n2 r 3 1 /II „ „ \ 11, ,3 

oo) = [ma T "ib) + ^ (^Y " - ^/'fej - Yl/'a + /'&)- g 

+ \{Pa + Z'b)' - ^(Ina + ln/3) - ^lnaln/3 + 6C(3) + ^ 



mt 



_3_ 
2^ 



^(/'a + /'6)-^(lna + ln/3)-5 



(B6) 



with the upper (lower) sign corresponding to vector (axial) current, and ^(3) = 1.20206 . . . 
As expected, only the transverse part of the vacuum polarization function is quadratically 
divergent for |g| oo. This divergent term, the expression of which is in agreement with 
the one obtained in Ref. [73], is the same for axial and axial- vector currents as expected from 
chiral symmetry. Moreover, as required by the Kinoshita-Lee-Nauenberg theorem [61,62], this 
term does not contain any mass singularity as nia^b tend to zero. As the vector part of the 
longitudinal component should vanish for nia = m,}), it must be proportional to (rrZa — m;,) or 
(m^ — ml), and since the axial- vector component can be obtained by changing the sign of one 
of the two masses, it must be proportional (rrZa + m;,) or (m^ — ml) (\nm'^/ml alone would 
have introduced mass singularities). This behavior is exhibited in Eq.(B6). 
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B.1.2 Scalar self-energies 

Referring to the results of Sec. 4. 2, I begin giving the one-loop scalar VPF in the limit of 
transfer momentum much larger or much smaller than the quark masses; in this case the 
self-energies read 



4 

— 
2 ' 



n-(»«<.) = -i+ln=^-l + l4±4ln4, (B7) 



sll+{s » ml,) = 1 - lln-J + 1. (B8) 
We can now consider the two-loop case, listing some limiting cases of the result of Eqs. (4.50). 
Equal masses 

In analogy with the one-loop case illustrated in Sec. 4.2.1, we can derive from Eqs. (4.20) the 
expressions of the self-energies of neutral scalar and pseudoscalar Higgs bosons, setting the 
quark masses to rUq, and using the appropriate couplings. Using x = 4a/(l + -^/l + 4a)^ with 
" = -iTT-l/s, 

n^,'^{s) = ^^^^shl{S„A,) (B9) 

with Sq = n+(s) - an-(s) and Aq = 11+ (s) + an-(s) given by^ 

3 , ,1/11 \ 11 9 

^^-^(1 + 12a) - 7 ( Y - S/fa + 42a - 36a/>J + -jPa " ^pl + 84a/>„ 

3 3 1 

- 36apl + - - 73a (1 + 12a) + -(1 + 4a)2(14a + 3)lna; 

8 4 2 

3 3 

+ (- + 14a + 29a2) In^ a; - 2(1 + 4a)t [^(a;^) - G{x)] 

- (1 + 2a)(l + 4a) [f{1) + F{x^) - 2F{x)] ; 

3 , , 1 /II \ 11 9 

A^- — {1 + 4a) - - ( — - 3pa + Ua - Uapaj + —pa - 3pl + 28a/>„ 

3 3 1 

- I2apl + - - 33a (1 + 4a) + -(1 + 4a)2(3 - 2a)lna; 

8 4 2 

+ (^ + 2a - 3a2) In^ a; - 2(1 + 4a)t [^(a;^) _ ^(a;)] 

- (l + 2a) [F(l) + F(a;2)-2F(a;)]. (BIO) 



^The expression of recently derived in Ref. [84] in the case of the SM Higgs boson is in agreement with 
our result. 
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One massless quark 

In the limit where one of the quarks is massless, = 0, the coefficients of 11^ and 11" vanish 
while takes the simpler form (with x = a/[l + a)) 



11 

— + 14a - 3pa - I2apa 



3 53 9 
- 3pl{l + 4a) + ---a- -(1 + 4a) + -(1 + aflnx (Bll) 

+ ^(1 + a)2(3 + 2a) In^ a; + (1 + af[F{x) - F{1)] + (1 + a)G{x) 
As expected, in this limit the expression of n+(s) is free of mass singularities. 
Case 5 > m^^ 

When the masses are very small compared to the momentum transfer, the coefficient of 11" 
vanishes and 11^ reads 

, 3 1 / ll\ 3 

+ 101„4-31n'4-|ln^ + |ln^!=i (B12) 

Case 5 < m^^ 

Finally, in the limit where the momentum squared is much smaller than the internal quark 
masses squared, the components 11+ and 11" read 

,61, , 3 mi 2 

SUA = -zrriA. H — (14m_|_ — 3m_|_/9_ — 3m_/9_ — 6m_|_/9_|_) oP- 

^ e 4 m_ 

3mi , / , 9 , \ 

+ i:^P- + ^+ [^P+ - 14/^+ - 7p- + ip+P- + -^pl + 30 j 

+ xV+ + m_ {^p'i_ - lp_ + 3/9+/9_^ ; (B13) 
n- = 4 + -(14 - 6/>+) - I'^pl + 30 - 14/>+ + 2p\ + + IpI, 

with p± = In ^ ± In ^ and m± = ± m^. 



B.2 Imaginary parts and partial widths 



In all the expressions for the VPF's the momentum transfer has been defined to be in the 
space-like region, < 0. When continued to the physical region above the threshold for the 
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production of two fermions, s > (rrZa + mbY, the self-energies acquire imaginary parts, related 
to the decay widths of the bosons into fermions, at this order in perturbation theory, by 



Imll{q^) = M T{q^) 



(B14) 



where M is the mass of the boson^. 

Adding a small imaginary part —ie to the fermion masses squared, the analytical continu- 
ation is consistently defined. From the general expressions of Eqs. (4.6,4.20), the imaginary 
parts can be straightforwardly obtained using 



ImlnXa b = 



Im\v? Xab = 2xln \Xa b\ 



In particular, for some of the integrals present in the two-loop VPF, we obtain^ 



J = -Iml = -2 [4Li2{xaXb) - 2Li2(a;a) - 2Li2(a;i,) + 21n \xaXb\ 



■K 



X ln(l — XaXb) — In |a;a| ln(l — Xa) — In \xb\ ln(l — Xb) 



(B15) 



J' = -Iml' = ln\xaXb\ i- + a + /3--X2] + 4A2 ln(l - a;„a;6) (B16) 

+ In I - (/3 - a + A2 ) ln(l - x^) - {a - (3 + A2 ) ln(l - Xb). 

2 p 

B.2.1 Vector boson self-energies 

From the complete one-loop expressions of Eqs. (4.6) we have the following imaginary parts 



Imll+{s) = 



(1 + « + /?)- -A 



X , 1 



2:mn+(s) = -A2 [A-l-a-/3], 



Iml[j,[s) = xA2 
Tmn^(s) = xA2 



(B17) 



or, equivalently. 



Imn.l'^{s) = ^sA2 



l + a + /3 --A ± 2^^ 

6 S 



X-l-a-p ± 2 



rriamb 



(B18) 



From the complete two-loop expressions of Eqs. (4.20) we obtain for TmHj, 



-Imll+ = -(11 + 19a + 19(3 + 12a(3 - 5A)A2 + -(a - (3) In 



(l-Xa) 



^At higher orders in perturbation theory, and for = M'^ , one should take into account the wave function 
renormalization of the decaying particle, including a factor 1 — li'(M^) in the denominator of the r.h.s. of 
Eq.(B14) 

^There are small differences with the formulae of Ref. [59] and [74] due to trivial algebraic simplification, 
and to the correction of a couple of misprints. 
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i^^-^^ [3(1 + a + + a + /3 + 2a/3 - A) + 4a + 4/3 + 26a/3] 
2 l + a + /3 i\ a 



A2^ (a-/3)ln- 



3 2 y ' /3 

+ 1 [(a + /3 - 2)A - 12a/3] J - ^ [3(1 + a + /3) - A] J'; 
-Iml[~ = (9 + a + /3)A2 + [4(1 + a + /3) + 2a/3 - A] In \xaXb\ 

TT 

+ (a - /3)Ai In ^ - 2(1 + a + /3) J - 4 J'; 
a 

-Imllt = -(9A - 9 - 9a - 9/3 - 12a/3)A2 - [(a + /3)A - 4a/3] J 
X 2 

+ ^ [(2A - 2 - a - /3)A - a/3(6a + 6/3 + 16a/3)] In |a;„a;6| 

+ ^(2A + 1 + a + (3)X2{(3 - a) In ^ + 2[1 + a + /3 - A] J'; 
z a 

-TmH" = 3(3 - a - /3)A2 + [4(1 + a + /3) - 6a/3 - A] In |a;„a;6| 

+ (a-/3)A2ln^-2(l + a + /3)J-4j'. (B19) 
a 

Equal masses 

2x f 3 1 

2:mn^(5) = — 5 <^ -A2(l - 6a) - 2a(4 + a) In |a;| - (1 - 4:a^)J 

- 2(1 - 2a) J'}, 
Stt f 3 1 

Imll-^{s) = — 5 <^ -A2 (1 + 12a + 4a^) + 2a(5 + 11a + 6a2)ln |a;| 
3 1^ 2 

- (1 + 2a)AJ - 2AJ'} , (B20) 
Iml[f{s) = 4x5a |^A2 (3 - 2a) + (3 + 4a - 6a2)ln |a;| - (1 + 2a)J - 2J'| , 

where 

J = -4 [2Li2(a;2) - 2Li2(a;) + 21n|a;|ln(l - a;^) - ln|a;|ln(l - a;)] , 
J' = (3 + 4a - 3A2)ln|a;| +4A2 ln(l - a;^) - 2A2 ln(l - a;), (B21) 
with A2 = a/1 + 4a. 

One massless quark 

Imlll'^{s) = X5 |l + ^a + ^a^ - ^a^ - 1(1 + a)2(4 - 5a) ln(l + a) 

(5a^ + 4a - 5)ln(-a) 

3 
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'-{2-a){l + aY 



2Lio 



a 



a 



ln(l + a)ln- 



-a 



a 



f r 9 

Iml[l'^{s) = -X5a(l + a)Ul + a) -- + (5 + 2a) ln(l + a) 
- (3 + 7a + 2a2)ln(-a) 



2(1 + a) 



2Li, 



a 



a 



ln(l + a)ln- 



-a 



a 



(B22) 



IJ ^ ' 1 

Close forms for Imll^'^^(s) in the general case rUa ^ rrib ^ have been also derived in the past 
by a number of authors [56,68,69] (in the first reference only the transverse part is given) by 
directly calculating the QCD corrections to the flavor changing decay of a vector boson. The 
results that we obtain here using a completely difl"erent method agree with those of Ref. [56,68] 
and also with Ref. [69] once some obvious mistakes in the integrals of their Appendix (Ji and 
J2) are corrected; see also Ref. [70]. 



B.2.2 Scalar self-energies 

We now consider the partial decay widths of the various Higgs bosons into quark-antiquark 
pairs. At 0{Gfj,as), the partial decay width of a charged Higgs boson into ud pairs is given 
by {s = M]j^) 

G^asMH+ 



T{H+-^ ud) 



2^2x3 



hllmll+{s) + hjlmll^{s) + 2KhJ^^^^^^ImU.- {s) 



The imaginary parts of 11^^ and 11 can be derived along the same lines as discussed previously, 
and one obtains for IniU.^ in the case of arbitrary masses 
1 



^(1 + a + /3)A2 + ("^ + a + /3 ) A + 5a/3 



(In \XaXh\ 



(1 + a + /3)(a - /3)A2 + In 



9 



(1 + a + - (1 + a + pfj - 2(1 + a + /3) j'; 



TmllKs) = 2:mn+(s) [ m„ ^ rrzd ]; 



(B23) 



-Imnr{s) = -2[A+2(l + a + /3)]ln|a;„a;6| 

- I2A2 - 2(a - /3)A2 In I + 2(1 + a + 13) J + 4 J'. 

From this formulae one can derive again the expressions of the hadronic decay widths in 
the previous special situations of physical relevance. In the limit where one of the quark is 
massless, 0, one has for Iml[^[s) 

= -2:mn+(5) = ^(1 + a)2 + (1 + a)(3 + 7a + 2a2) In — ^ - 2(1 + a)2 
TT 2 1 + a 



ln(l + a) 



2Li2 



a 



a 



a 



ln(l + a)ln- 



-a 



a 



(B24) 
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In the case of scalar and pseudoscalar neutral Higgs bosons, the partial decay widths T[S, A 
qq) are given by 



T[{S,A)^qq] 



2^2x2 X « 



where Im[Sq, Aq) = Iml['^[s, (3 = a) ^ almH. {s,l3 = a) are given by 



(B25) 



X 

1 



ImSq = -(1 + 4a)2(14a + 3) + {58a^ + 28a + 3)ln|a;| - 4(1 + 4a)/: 



" 2 
3 



where 



-ImAq = -(1 + 4a)2(3 - 2a) + (3 + 4a - 6a^)ln|a;| 
X 2 



£ = Ij' + 7(1 + 2a)J. 
2 4 



4:£ 



(B26) 



(B27) 



Finally, we note that in the limit where the quark masses are much smaller than the Higgs 
boson masses, the QCD corrections to the decay widths will exhibit the well known logarithmic 
behavior [78, 79] which, because of chiral symmetry, is the same for the scalar, pseudoscalar 
and charged Higgs boson 



-Imlln 



X 



(B28) 



The large logarithm is responsible, for example, of the very large QCD corrections to the 
decay width of the SM Higgs boson into bb, which are of the order of 35% for a Higgs mass of 
only ~100GeV. However, it is possible to resum the large logarithmic terms: this is equivalent 
to replacing the on-shell quark masses by running masses defined at M| in the first order 
amplitude [79]. We can see that explicitly, in the case the MS mass mq(M$) is used. The 
relation between on-shell and MS mass is shown in Eq.(4.13), and we can write the partial 
width of the scalar $ in the limit of vanishing quark masses as 



with 



2 a, /9 

3^(2+^^^ Ml 



3G'^M$ ,2 



4V2x 



(B29a) 



(B29b) 



As hq is proportional to niq, shifting the mass according to Eq.(4.13) cancels the logarithm in 
the r.h.s. of Eq.(B29a). 



Appendix C 



Ward identities for electroweak 
VPF's 



In this appendix I report some Ward identities for the two-loop self-energy diagrams that 
contribute to 0(G^mjM|). They were derived by the method described in Sec. 5.3 and have 
been used in the calculation of Sec. 5.2. We can write the W and Z self-energies at = as 
the sum of different contributions: 



where the A'^"'^(0)'s are given in Eqs. (5.39), and represent the set of diagrams of Fig. 5.2, 
i.e. diagrams contributing to 0{G'^mf). A" and A"'°" represent all the other contributions 
of Feynman diagrams involving the top quark. Some of the graphs contributing to them 
can be obtained from Fig. 5.2 replacing the scalar lines with vector boson lines, the others 
are displayed in Fig. C.l. The diagrams involving fermionic insertions on a scalar propa- 
gator (Fig. C. lb) formally contribute to leading 0{G'^mf). However, when the appropriate 
mass counterterms of the physical and unphysical scalars are included, they are effectively of 



I first consider only the irreducible two-point diagrams that contain the Higgs boson and the 
top quark, but do not contribute to 0{G'^mf). They are depicted in Fig. C.l, where the 
dashed line represents a Higgs boson. The short-hand notation (...) stands for (0|T* . . . |0), 



A{0) = A^''"^{0) + A"{0) + A"°"(0) 



(CI) 



0{Glm- 




C.l Subleading with Higgs 




and = J d^x. We have explicitly checked these Ward identities as 



^f.(O) = 
n 2c4 dq' 




{P - Ml)[{k - qY - Ml] 



,ik-x—iq-y 



{ S2{y)Jn^)Si{0) ) 
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9 
"4c4 



// 

Jk Jx 

Jk Jx 



oik-i 



{P-Ml){P-Ml) 

ik: {j^{x)j,c.m 



n c' 



-{Ml-M: 



i d 



{P-Mlf{P-Ml) 
1 9' f !j^-^S2{x)S2m-i{Sm 



n 2c4 Jk 



{P-Mlf{P 



(C2) 



n 4 dq' 
-II 

4 JkJx 



/// 

J k V X V y 



-iq-y-\-ik-x 



{[S\y)j^{x)- h.c.)5i(0)) 



- Ml) [{k - qf 
'^■^Sr{x)Sr{Q)) g^Ml 



II 

Jk Jx 



^'=-(J«(a;)jt?(0)) 



n 



{Ml - M, 



{P-Ml){P-Mlf 2n JkJx{P-Mlf{P 
/,e'^-(5t(a.)5(0))-i(5i(0)) 



Ml) 



Jk 



{P 



Ml)\P 



Ml)' 



(C3) 



C.2 Subleading without Higgs 



Here I consider the two-loop diagrams involving the top quark but not the Higgs boson that 
contribute to 0(G^mjM|). Some infrared divergences (that cancel in the sum) have been 
regulated by a photon mass A. Note that the last two lines of Eq.(C5), corresponding to 
the diagrams in Fig. C.le and to some residual terms coming from the four- and three-point 
functions, do not contribute to 0(G^mjM|), as can be easily verified. 



AnoH 



(0) 



9 



4 d 



8c2 9g2 



// / / 

V k V x-\ J Xo V y 



,iqy-\-ik{xi -X2) ^ 



X {S2{y)S2{0) 
d 



Jw{xi)Jw\{x2) , Jz{xi)Jz\{x2) , S^J^{X1)J^X{X2) 



- Ml 



k^ 



9 

2nc2 



— // 

Jk Jx,y 



^iqy—ikx 



{S2{y) {S^{x)j^{0)+h.c.)) 



+ Mlg^- 



I 

Jk. 



n 

Akx 



// 

Jk Jx 



Ml) 



{k-' - Ml){{k - qf 

Akx 



(Jg(a;)Je.^(0)) 

k,x {k-'-MiY 



{S^{x)S{(^)) 



{k^ 



Ml)' 



{c' + 2s') [ 
Jk 



{Si{0)) 

k {k^-MD^- 



(C4) 



A7w{0) = 
X {S{y)S\0) 



II I I 

J k J xi J X2 J y 



,iqy+ik{xi -X2) 



4 9g2 

Jw{xi)Jw\{x2) , Jz{xi)Jz\{x2) , S'J^{X1)J^X{X2) 



Ml 



c2(F - M|) 



fc2 
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— 

(a) 



(b) 



(c) 



(d) 




(e) 



Figure C.l: Some of the self-energy diagrams contributing to 0(G^mjM|). Solid lines represent 
fermions, dashed lines scalars, and wavy lines vector bosons (see text). 



^4 ^ » ^-iqy+ikx 

Td^ jk,.,y {P - Ml)[{k - qf - M, 
.9 



4 Q . e-i-iy+i^'^ q^{(s\y)S{{))-h.c) 
Jk 



g„ ((5t(t/)J«(0) + h.c.) S2{x)) 



2 9g2 Jk, 
.g' d 



[{k - qf - M, 

Ak-x—iq-y 



c^{k^ - Ml) 



k^ 



) 



2 dq' JJJ^[(k-qy-Ml] 



{{S\y)j^{0)+ h.c.) 



,2 A/f2^ ''zolK-I') p Jia{-1-) 



2^ )/..(o))/^ 



2 /(4 - n)fc2 + nA^ 2(n - 5) 



6 — n 



\ (fc2-A2f k^{k'-Ml) 
AMl{Ml - Mlf 6 - n A-n 



{k'-Ml)\ 



Xk' - Mlfik' - Mlf {k'-Mlf {k'-Mlf 



2(n- 5) 



_ M|)(fc2 _ Ml) (^2 _ ^2)(^2 _ Ml)' 
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n JkJx {P -Mlf{P -Ml) 



{Ml - Ml)\nMl + (4 - n)k') ^ ^ - P 



2c^ 



ik' - Ml)' 
9' f Le''--{S2{x)S2m-i{Sr{0)) 



Ac' 



2n Jk {k' - Ml){k' - Ml) 
^ f J^e''^-{S\x)Sm-l {S^m 
n Jk 



{k' - Ml) 
4 4 



(n-6)c2 + 2i2 
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